The completion method is based on the calculation of the neuronal mass densities inside the slices and extrapolating these densities to the area outside the slices. For this extrapolation it is assumed that the mass densities of the neuronal arborizations are axial symmetric (Figure 2A). This is a plausible assumption for pyramidal neurons, which have the apical main shaft as symmetry axis. Using an axial-radial coordinate system, the neuronal mass is determined for a given height and a given radius, thus summing all the mass in a ring as is illustrated in Figure 2A. For a fully intact neuron, the radius of the integration ring ranges from zero up to the maximal radial extension of the arborizations. The total mass distribution of a neuron becomes a function of height and radius. Evidently, in the case of pyramidal neurons, these distributions can also be obtained separately for axons and (apical and basal) dendrites. For a sliced neuron, its radial extension in the XY direction of the slice can become as large as in the intact neuron, but in the Z-direction it is limited by the boundary planes of the slice. Also the integration ring will be complete for small radii but becomes incomplete at one or two sides for larger radii. In the last case the ring is reduced to two separate ring segments (Figure 2B). The estimation of the mass densities is then restricted to the remaining parts of the integration rings.

With the assumption of axial symmetry, the densities in the fractional parts of the integration rings inside the slice can be extrapolated to their complementary parts outside the slice. The volume of the integration ring inside and outside the slice is fully determined by the position of its center (i.e., the symmetry axis), the radius of the integration ring, and the thickness of the slice. Three situations can be distinguished. The ring can be fully contained in the slice, it can extend the slice on one side, or it can extend the slice on both sides (Figures 2C–E). The neuronal masses can only be obtained for the volume fraction of the ring inside the slice. For an integration ring with radius r (pointing to the center of the ring), the volume fraction F(r) of the ring area (or volume) within the slice can be expressed in terms of the thickness of the slice (T), the position of the ring center within the slice (H) and the radius of the ring (r).

An estimate of the mass in a full ring M_ring(r, h, δr, δh) with radius r and at height h and with radial thickness δr and height δh can be obtained by dividing the experimentally observed neuronal mass M_obs(r, h, δr, δh) in the part of the ring inside the slice by the ring volume fraction. Thus,

For the mass density in the ring D_ring (r, h) we now obtain Dring​(r, h)=Mring(r, h, δr, δh)Vring=Mring(r, h, δr, δh)2πrδrδh with V_ring the volume of the full ring. In our analysis, integration rings with a radial thickness of δr = 1 μm and height δh = 1 μm are used.

Neuronal arborizations may extend beyond the boundaries of a slice but more distal parts may bend back into the slice (Figure 2F). Such branches have lost their parent branches and are called orphan branches. Orphan branches are no longer connected with the proximal parts of the arborization inside the slice. As the reconstruction procedure quantifies only the contiguous parts of the arborization inside the slice, the orphan branches are lost. In the calculation of the mass densities from the contiguous part inside the slice, this then results in an underestimation of the neuronal mass. As the completion procedure extrapolates from this neuronal mass, it does not solve the problem of lost orphan branches.

For the completion procedure it is important to know the thickness of the slice and the position of the cell body with respect to both boundary planes of the slice (Z-axis). Although the thickness of the slices is provided in most reconstructions, this is not the case for the cell’s position in the slice, so that we had to estimate the cell’s position from the reconstructions themselves. For this estimation we used the fact that all cut terminal tips have approximately the same Z-coordinate. To this end we calculated the Z-coordinates of all terminal tips of both axons and dendrites in the reconstruction and obtained a frequency distribution along the Z-axis. When the distribution showed a sudden increase in frequency (number of tips) at one or both of its ends, this was interpreted as the effect of cutting. When a frequency increase was observed at the low-Z end of the distribution, the cut part of the cell’s arborization was aligned at the low-Z end of the slice. When a frequency increase was observed at the high-Z end of the distribution, the cut part of the cell’s arborization was aligned at the high-Z end of the slice. Finally, when a frequency increase was observed at both sides of the distribution, it indicated the actual thickness of the slice. The thickness information and the estimated cell positions were used in the completion procedure.

The mass density completion procedure was applied to the arborizations of individual neurons because the procedure depends on the position of the individual neuron inside the slice, which varies from neuron to neuron. To obtain a population average of the estimated mass densities, the cells were aligned by their somata and all the cells were rotated in such a way that their apical main stem was pointing into the Y-direction.

For the alignment of the apical dendrite, the orientation of the apical main stem was needed. As the data formats in the reconstruction files generally do not distinguish between apical main stem, and apical tuft and oblique branches, this information needed to be derived from the supplied data. To this end, an iterative procedure was applied of pruning terminal line pieces off the apical dendrite, until it was reduced to a single segment with one terminal line piece. The tip coordinate of this terminal line piece together with the exit coordinate from the soma of the apical dendrite provided the alignment line piece with which to rotate the whole cell in the XY plane (i.e., around the Z-axis) in such a way that the projection of the alignment line piece onto the XY plane was pointing into the Y-direction. Alternatively, a rotation in the XYZ space could have been applied so that the alignment line pieces were all pointing into the Y-direction. However, such a rotation would also have changed the Z-coordinates of the terminal tips in the arborizations and would have hampered the use of tip coordinates for estimating the cell’s position in the slice. As a consequence, the orientation of the alignment segment could maintain a slightly tilted angle with respect to the Z-axis, while the Z-axis itself was taken as symmetry axis for the calculation of the mass distribution. Given the wide spread of arborization mass, the effect of a possibly slightly tilted orientation on the final estimated mass distribution was assumed to be negligible.

Crucial for the completion procedure is its validation, i.e., whether the masses of the completed neurons are equal to those of the original non-sliced neurons. Such a validation is not possible for experimental data sets of sliced neurons, but can be done for a set of model-generated neurons. For this validation, we used a set of 50 neurons, generated with our NETMORPH simulator (Koene et al., 2009). These model neurons were subsequently sliced according to several slice thicknesses, followed by the density field completion procedure. The masses of the original and of the completed neurons were subsequently compared. Because in this case the mass of orphan branches was also known, the comparisons were made from completions with and without inclusion of orphan branches.

For one of the NETMORPH-generated neurons an example is given in Figure 3 to demonstrate the impact of slicing on the remaining morphology within the slice.

To validate the completion procedure, the 50 NETMORPH generated neurons were artificially sliced with slice thicknesses of 100, 200, 300, 400, 500, 600, 1000, or 2000 μm. First, the neurons were aligned with their apical main stem pointing into the Y-direction, and with the slicing perpendicular to the Z-axis. For the position of the neurons in the slice, two choices could be made. With the somata positioned in the center of the XY plane of the slice, the Z-coordinate was put either in the center of the slice (thus with equal distances to both cutting sides of the slice), or was uniform randomly selected within an 80% Z-range (between 10 and 90%) or within a 60% Z-range (between 20 and 80% of the Z-range). After slicing, the axial-radial mass distribution was obtained for each individual neuron followed by the completion procedure, which was applied excluding or including orphan branches. The results of the completion procedure were finally summed and averaged over all the NETMORPH neurons in the data set.

The results for a slice thickness of 200 μm, 80% range random Z-coordinates of the somata, and ignoring orphan branches are shown in Figure 4. Figure 4A displays the axial-radial mean dendritic mass distribution, in the top part for several axial positions (summed per 50 μm height steps) and in the bottom part summed over all heights. The dashed histograms present the masses of the sliced dendrites, while the solid lines present the completed dendritic mass distributions. Both distributions coincide at small radial distances from the axis, but the completed ones become increasingly higher at larger radial distances. Evidently this is caused by the larger correction factors at higher radial distances. The distribution for the non-sliced original full dendrites is also included in the bottom panel of Figure 4A as a red solid line. The total mass of 5264 μm of the dashed (sliced) distribution deviated 19.7% from the original dendritic mass of 6552 μm, while the completion procedure resulted in a mass of 6399 μm, deviating only 2.3% from the original mass (Table 1). Figure 4B displays the mean mass distribution as a function of the radial distance to the soma of the sliced dendrites (dashed), of the completed ones (black solid line), and of the original non-sliced dendrites (red solid line). This distribution is related to a 3D Sholl diagram which counts the number of intersections with a set of concentric spheres (in practice, however, the Sholl method is most frequently applied to 2D projections, see Uylings and Van Pelt (2002) for a discussion on Sholl diagrams). The tail in both distributions originates from the apical dendrites. Figure 4C shows the positions of the somata in the Z-direction of the slice (filled circles) as well as the distribution of Z-coordinates of the terminal tips per neuron. The effect of cutting is clearly seen in the accumulation of terminal tips at the boundary layers of the slice. Figures 4D–F show the findings for the mean axonal mass distribution. Note that panel 4D extends over larger axial (Y-axis) and radial distances than panel 4A for the dendritic mass distribution. The summed radial distribution shows a substantial difference between the sliced axonal distribution (dashed) and the completed one (black solid line), with a mean total axonal mass per neuron of 4239 μm for the sliced axon and of 8531 μm for the completed mass, while the original axon (red solid line) had a total length of 10655 μm. Thus, in 200 μm thick slices, 60.2% of the original axonal mass was lost, while the completion procedure reduced the axonal loss to 19.9% (Table 1). Comparison of the black and red solid lines in the bottom panel of Figure 4D clearly shows the level of axonal recovery obtained in the case of 200 μm thick slices. The difference in the level of recovery between dendrites and axons appears to be related to the fraction of conserved mass in the slice (see the results shown in Figure 5E).

Table 1 summarizes the effects of slicing on the dendrites and axons of 50 NETMORPH-generated neurons, depending on slice thickness and soma positions in the slice. Also are shown the results of the completion algorithm, taking or not taking into account the orphan branches. The original mean (apical and basal) dendritic and axonal lengths are 6552 and 10655 μm, respectively.

As expected, slicing has a significant effect on the length of axons and dendrites. For instance, centrally located neurons in 100 μm thick slices loose 34.5% of their dendritic length and 85.9% of their axonal length; in 200 μm thick slices the neurons loose 12.0% of their dendritic length and 61.1% of their axonal length; and in 300 μm thick slices they loose 4.0% of their dendritic length and 37.3% of their axonal length. A significant recovery of the lost parts was obtained by applying the completion procedure. For instance, when orphan branches were not included, dendritic loss in 100 μm thick slices was reduced from 34.5% down to 7.4%, and axonal loss from 85.9 to 61.1%. When orphan branches were included, dendritic loss was further reduced to 4.5% and axonal loss to 7.0%. In 300 μm thick slices, dendritic recovery of centrally located neurons excluding orphan branches resulted in 0% dendritic loss and 7.9% axonal loss, and including orphan branches in 0% dendritic loss and 0.8% axonal loss.

When the neurons were randomly placed in the slice, the losses were slightly higher for the dendrites, but more or less similar for the axons because the larger extent of the axons makes them less sensitive to the precise position in the slice. Also with random placement the completion procedure was able to greatly reduce the loss. For instance, in the case of a Z-range of 10–90%, the dendritic loss of 41.1% in 100 μm thick slices was reduced by completion without orphan branches to 9.3%, and with orphan branches to 3.1%. Axonal loss of 86.1% was reduced by completion without orphan branches to 61.3%, and with orphan branches to 3.5%.

The results in Table 1 show that the completion procedure that takes orphan branches into account is able to fully recover the original mass, with loss values around zero. Small negative loss values also occur, indicating an overcompensation, which can be expected when, as a result of statistical fluctuations in the spatial distribution of the arbors, the mass densities inside the slice are somewhat larger than outside the slice. The full recovery of dendritic and axonal mass when orphan branches are included can be considered as a validation of the completion procedure. In the absence of knowledge about orphan branches, which is usually the case, the completion procedure is still able to recover dendritic mass to values very close to the original mass, with a deviation of less than 2.3% in the case of 200 μm thick slices, and less than 1.6% in 300 μm thick slices.

Intracortical axons extend their branches at large distances from the slice and their loss in 100 μm and 200 μm thick slices is substantial (up to 85.9 and 61.1%, respectively). Although the completion procedure significantly reduces these losses (to 61.1% in 100 μm thick slices, to 25.6% in 200 μm thick slices, and to 7.9% in 300 μm thick slices), knowledge of orphan branches is required for a full recovery.

The loss of dendritic and axonal mass by slicing depends on the spatial extent of dendritic and axonal arbors in relation to the slice thickness. The NETMORPH neurons have a mean radial center of dendritic mass at 70 μm, and a mean radial center of axonal mass at 204 μm. The dependence of mass loss on slice thickness as well as the results from the completion procedure is depicted in Figure 5, displaying the mass loss of neurons with their somata in the center of the slice. Clearly is shown that for 300 μm thick slices, completion without orphan branches results in a full recovery of the dendrites (dashed red curve in 5A), with a loss of 0% (Table 1), and an almost full recovery of the axons (dashed green curve in 5B), with a loss of 7.9% (Table 1).

Figure 5C shows a scatter plot of the mass losses of the axonal and dendritic arborizations of each individual neuron vs. the radial center-of-mass of these arborizations after slicing. These data were obtained for neurons with their somata uniform randomly positioned in a 80% range of the slice thickness of 300 μm. Thus the closest distance of a soma to the cutting plane is 30 μm. Figure 5D shows a scatter plot of the recovery results vs. the radial center-of-mass of the individual axons and dendrites within the slice. Figures 5C,D reveal a large scatter in the individual data points, originating from the large variability in morphologies of the axonal and dendritic arborizations and the effect of cutting. The larger extent of axonal arborizations in comparison with dendritic arborizations (as reflected by the larger values for their radial center-of-masses) causes also larger values for their mass losses by slicing, as shown in Figure 5C. Nevertheless, for these 300 μm thick slices, there is almost full recovery of both axons (95.7%) and dendrites (98.3%). Note that the axonal recovery is better than shown in Figure 5A because of the different positions of the somata in the slices. The uniform random positions of somata are more realistic than a central positioning. The findings in Figure 5 thus show that the completion procedure results in (almost) full recovery of the mass distributions of axonal and dendritic arborizations when they were cut by 300 μm thick slices. The amount of scatter in the data points appears to depend also on the locations of the somata in the slices. In Supplementary Figure S5, the results are also shown for centrally located neurons, and for uniform random placements in a 60% range of the 300 μm thick slices.

Whether recovery results are related to the fraction of the original neuronal mass conserved in the slice is shown in Figure 5E using the data in Table 1 (excluding orphan branches). Both the axonal (green) and dendritic (red) data points show a clear dependency of recovery result on the fraction of conserved mass. The data points are labeled by the section thickness and the positioning scheme used. Averaged over the three positioning schemes we observe for the dendrites, that a conserved mass of 62.2% (100 μm) relates to a recovery of 91.8%, a conserved mass of 84.1% (200 μm) relates to a recovery of 98.1%, and a conserved mass of 92.5% (300 μm) relates to a recovery of 99.2%. For the axons, a conserved mass of 14.3% (in 100 μm thick slices) relates to a recovery of 39.9%, a conserved mass of 39.1% (200 μm) relates to a recovery of 76.4%, and a conserved mass of 61.6% (300 μm) relates to a recovery of 95.3%. Clearly is shown in this figure how the red dendritic data points for 100–300 μm thick slices intermingle with the green axonal data points for 300–1000 μm thick slices. Thus dendrites and axons show similar relationships obtained for different slice thicknesses, which leads us to the conclusion that recovery results relate to the fraction of conserved mass, while this relation is independent of the slice thickness.

The 11 rat parietal cortical L2/3 pyramidal neurons were reconstructed from 300 μm thick slices (Shepherd and Svoboda, 2005). Figure 6 shows a selection of four of these neurons as projections on the XY plane and the YZ plane. A full display of all 11 neurons is shown in Supplementary Figure S1. While the XY projections show the full extent of the arbors, the YZ projections clearly show the accumulation of terminal tips at the intersection of the XY and YZ planes.

The axial-radial mass distributions of the reconstructed neurons after mass completion are shown in Figures 7A,D for (apical and basal) dendrites and axons, respectively. Particularly the summed radial distributions in Figures 7A,D and the radial to soma distributions in Figures 7B,E show the extent of the corrections, which were much larger for the axons than for the dendrites. Figures 7C,F clearly show the cut endings, which appear as an accumulation of cut terminal tips at one side of the arborizations. This accumulation was used as a criterion to align the neurons within their slices. The total length of the completed mass distributions is shown in Table 2A. The mean dendritic length increased from 7865 to 9289 μm, and the mean axonal length from 4692 to 9118 μm. The radial center-of-mass of the reconstructed dendrites was equal to 76 μm and that of the reconstructed axons 195 μm. Based on the validation results with the NETMORPH neurons (Figure 5) and the thickness of the slices in the Svoboda data set (300 μm), we may expect an almost full recovery of the axonal and dendritic mass. This implies that by slicing 15% of the dendritic mass was lost, and 49% of the intracortical axonal mass (within an uncertainty range of a few percent).

The 33 Wistar rat somatosensory cortical L2/3 pyramidal neurons were reconstructed by Wang et al. (2002) from 300 μm thick slices, and are shown in Figure 8 as projections on the XY plane and the YZ plane. While the XY projections show the full extent of the arbors, the YZ projections clearly show the accumulation of terminal tips at the intersection of the XY and YZ planes.

The axial-radial mass distributions of the original neurons and after mass completion are shown in Figures 9A,D for (apical and basal) dendrites and axons, respectively. Particularly the summed radial distributions in Figures 9A,D and the radial to soma distributions in Figures 9B,E show the extent of the corrections, which were much larger for the axons than for the dendrites. Figures 9C,F clearly show the cut endings as an accumulation of cut terminal tips at one side of the arborizations, which was used as a criterion to align the neurons within their slices. The length of the completed mass distributions is shown in Table 2B. The mean dendritic length increased from 3790 to 4567 μm, and the mean axonal length from 3211 to 6177 μm. The radial center-of-mass of the reconstructed dendrites was 59 μm and that of the reconstructed axons 160 μm. Based on the validation results with the NETMORPH neurons (Figure 5) and the thickness of the slices in the Markram data set (300 μm), we may expect an almost full recovery of the axonal and dendritic mass. This implies that by slicing 17% of the dendritic mass was lost, and 48% of the intracortical axonal mass (within an uncertainty range of a few percent).

The data from Parnavelas and Uylings originate from a study of (basal) dendritic development in rat visual cortex pyramidal neurons (Uylings et al., 1994). Golgi stained dendrites were reconstructed from slices with a thickness of about 120 μm. Reconstructions were obtained from 153 layer 2/3 pyramidal cells at different ages of cortical development, i.e., at 10, 14, 18, 24, 30, and 90 days postnatal (PN). Figure 10 shows a selection of the 24 L2/3 pyramidal reconstructions at the age of 90 days PN. A full display of all the neurons is given Supplementary Figure S4.

The results of the completion procedure for the 90 days PN data set are shown in Figure 11. With a mean length of 2022 μm for the reconstructed basal dendrites and 2532 μm for the completed masses, the mass loss by slicing becomes 20%. However, the validation study (Table 1) has shown that for 100 and 200 μm thick slices dendritic recovery still leaves a deficit of 7.1 and 1.7%, respectively, for somata within a 60% range of the slice thickness. Based on these validation findings, we may expect that the outcome of the completion procedure of 2532 μm deviates less than 7% from the original mean dendritic mass, thus in the range of 2532–2723 μm. The outcomes of the completion procedure for the other age groups are listed in Table 2C.

The problem of incomplete reconstructions by slicing has recently also been studied by Hill et al. (2012). The algorithm they devised for morphological repair was to derive a statistical growth model from the intact parts of the arborizations and then to regrow the cut portions. Using Bayesian spatial distributions, cut dendrites were regrown point by point. Axons were separately repaired by pasting subtrees from the intact parts. Invariance to axial rotations was also assumed. While Hill et al. (2012) attempted to recover the individual branches outside the slice volume, our approach aimed at recovering the population mean axonal and dendritic mass density fields by extrapolating the observed mass distributions within the slices to the outside space by simple geometrical relationships. Our approach is not model-based, and the only assumption used is that of axial symmetry in the mass distributions.

For the validation of the completion procedure we used a set of 50 neurons generated with our NETMORPH simulator. The neurons were subsequently sliced with different slice thicknesses, and subjected to the completion procedure. Comparison of the sliced masses and completed masses then showed how well the completion procedure was able to recover the lost masses. It turned out that a complete recovery was indeed possible provided that all the original mass of the sliced neurons within the slices was used. This outcome validated the completion procedure. However, the original mass within a slice includes the contiguous part seen from the soma, as well as the mass of orphan branches that lost their connection with the contiguous part within the slice. Because the orphan branches are not included in experimental reconstructions, they remain a missing part of the reconstructed neurons and therefore affect the degree of mass recovery. However, the number of orphan branches depends on the slice thickness and the spatial extent of the arborizations. From the validation data it turned out that this missing part for dendrites in the case of 300 μm thick slices was negligibly small, while for intracortical axons the missing part was less than 8% for centrally located neurons (Figure 5). For randomly located neurons in 300 μm thick slices a full (i.e., better than 98%) dendritic recovery and an almost full axonal recovery (i.e., better than 95%) was obtained by the completion procedure.

The validation study revealed that the mean loss of mass in a population of neurons depends on the positions of the somata within the slice. However, a systematic trend could not be derived from Table 1, because the outcomes for the three different location options also depended on the slice thickness and the extent of arborization (dendritic or axonal). Apparently, these geometrical parameters all play a role in the final outcome.

Recovery results appear to relate to the fraction of conserved mass in the slice. Both the dendritic and the axonal data show a similar relationship. In particular the observation that dendritic data points in this relation for 100 μm thick slices coincide with axonal data points for 300 μm thick slices, and for 200 and 300 μm thick slices are intermingled with axonal data points for 400–1000 μm thick slices shows that this relationship is independent of the slice thickness itself. Thus, the fraction of conserved mass in the slice determines the level of recovery obtained. For instance, a recovery result better than about 95% requires a conserved mass in the slice larger than about 60%. This relation may be of practical value as it provides guidance for the required slice thickness which will be different for axons and dendrites. Whether a similar relationship between conserved mass and recovery result is obtained when the method is applied to experimentally fully reconstructed neurons instead of NETMORPH generated neurons is still an open question.

Based on the positive outcomes of the validation study, we have analyzed three data sets of reconstructed neurons, two of which were obtained from the NeuroMorpho.org data base and one was provided by one of the authors of this study.

The Svoboda data set (Shepherd and Svoboda, 2005) resulted in a recovery of the mean dendritic length of 9289 μm (from the actual measured value of 7865 μm, indicating a loss by slicing of 15%). Axonal length was recovered up to 9118 μm (from the actual measured value of 4692 μm, indicating a mass loss of 49%), but this outcome may still be about 5% from the true value (see discussion in the previous two paragraphs).

The Markram data set (Wang et al., 2002) resulted in a recovery of the mean dendritic length of 4567 μm (from the actual measured value of 3790 μm, indicating a loss by slicing of 17%). Axonal length was recovered up to 6177 μm (from the actual measured value of 3211 μm, indicating a mass loss of 48%), but this outcome may still be about 5% from the true value (see discussion in the previous two paragraphs).

The Parnavelas-Uylings data set (Uylings et al., 1994) resulted in a recovery of the mean basal dendritic length of the 90 days age group of 2532 μm (from the actual measured value of 2022 μm, indicating a loss by slicing of 20%). These dendritic reconstructions were, however, made from 120 μm thick slices and the recovery result may still be about 7% from the true value (see Table 1, dendritic deviation for 100 μm thick slices with random positions in 20–80% range).

The outcomes of the three data sets are not directly comparable, because of different ages of the rats used (25–36 days PN for the Svoboda data, 13–15 days PN for the Markram data, and 90 days PN for the Parnavelas-Uylings data), and the restriction in the last data set to basal dendrites only.

It is interesting to note that the Svoboda and the Markram data set show similar mass losses for axons (48–49%) and dendrites (15–17%) in 300 μm thick slices. Apparently, the mean mass loss is not so sensitive to the difference in the mean radial center-of-masses of the cut axons (195 and 160 μm) and cut dendrites (76 and 59 μm) in the Svoboda and Markram data set, respectively. This is also in line with the more or less uncorrelated scatter of the individual NETMORPH neuron data in Figure 5C. Nevertheless, the axonal and dendritic losses do differ significantly in all these cases.

The described completion/recovery method, and thus the quality of its outcomes, is independent of homogeneous tissue shrinkage. The 3D metrical properties of the arborizations, however, are affected by shrinkage. In quantifying the 3D geometry of neuronal arborizations, one needs to take tissue shrinkage into account, which occurs in histological and staining procedures. The extent of shrinkage is different for different staining techniques, and is even different for different Golgi techniques. For the Markram data set based on HRP staining, Wang et al. (2002) reported a 25% shrinkage of the slice thickness and ~10% anisotropic shrinkage along the X- and Y-axes. Only shrinkage of thickness was corrected. Shrinkage correction in the Svoboda data set (Shepherd and Svoboda, 2005) was not reported. No shrinkage correction was applied for the Parnavelas-Uylings data set (Uylings et al., 1994). Uylings et al. (1986, 1989) summarized tissue shrinkage values for different Golgi staining procedures. Linear shrinkage between 5 and 10% was reported for Golgi-Cox and rapid Golgi staining, and of 5–20% for Golgi-Kopsch staining. Thus, depending on histological and staining procedures, tissue shrinkage is an important issue in quantitative studies of 3D neuronal arborizations. Correction for shrinkage is only possible when the actual amount of anisotropic shrinkage in X, Y, and Z direction is known.

A key assumption in the completion procedure is axial (rotational) symmetry in the distribution of axonal and dendritic mass. Because such symmetry is often assumed for pyramidal neurons, we have selected cortical layer 2/3 pyramidal neuron reconstructions for this study. While this assumption may be reasonably valid for basal dendritic arborizations that locally innervate space, it may not be valid for non-pyramidal dendrites (Parnavelas and Uylings, 1980) and axons that extend their arbors not only to local but also to remote locations. Clearly, single slice axonal reconstructions visualize only the local part of axonal arbors. Because the completion procedure recovers only this local axonal part, the axial symmetry assumption may still be valid.

The paper dealt with the estimation of the axial-radial distribution of axonal and dendritic mass. The calculation of axonal and dendritic densities is a straightforward extension and proceeds by dividing the mass in an integration ring (see Figure 2) by the volume of the ring (see Materials and Methods; see also Van Pelt and Van Ooyen, 2013).

Recently, we have shown that potential synaptic connectivity between neurons in a network can be estimated from their axonal and dendritic density fields (Van Pelt and Van Ooyen, 2013; McAssey et al., 2014). Thus, for constructing neuronal networks and their inter-neuron connectivity one does not need as many neuronal reconstructions as there are neurons in the network, but can use the population mean density fields instead. The availability of density fields for a variety of neuronal cell types is thus important. While for a large variety of cell types reconstructions have become available in open-access data bases, their incompleteness hampers a full use of the data. When with our method full density fields can be recovered from incomplete single-slice reconstructions, the open-access data become even more valuable, as they now can also be used for building neuronal networks and connectivity studies.

It has to be noted that the estimation of the number of synapses not only requires the number of potential synapse locations but also the probability that synapses actually are formed at these locations. A detailed EM study of the hippocampal neuropil by Mishchenko et al. (2010) showed that this probability was variable and dependent on ultrastructural details, such as dendritic circumference and actual axo-dendritic touches. Helmstaedter (2013) emphasizes in his review on dense neural circuit reconstruction that in “mapping neuronal circuits, it is important to detect synaptic contacts between neurons, but it is in many cases even more important to be able to exclude synaptic connectivity between neurons to determine the structure of a wiring diagram,” Clearly, the estimation of potential synapse locations is only one, but still crucial, factor in estimating synaptic connectivity.

Nevertheless, the realism of network connectivity estimates based on overlapping axonal and dendritic arborizations has recently been demonstrated by Hill et al. (2012) and Van Ooyen et al. (2014). In a statistical study, Hill et al. (2012) found that random alignment of axonal and dendritic arbors provides a sufficient foundation for specific functional connectivity to emerge in local neural microcircuits. In a computational study, Van Ooyen et al. (2014) found that the synaptic connectivity emerging between neurons that grow out in the absence of any guidance cues showed a good agreement with available experimental data on spatial locations of synapses on dendrites and axons, number of synapses by which neurons are connected, connection probability between neurons, distance between connected neurons, and pattern of synaptic connectivity.

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.