Skip Navigation


Cerebral Cortex Advance Access originally published online on December 15, 2005
Cerebral Cortex 2006 16(10):1474-1486; doi:10.1093/cercor/bhj085
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrowOA All Versions of this Article:
16/10/1474    most recent
bhj085v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (2)
Google Scholar
Right arrow Articles by Schüz, A.
Right arrow Articles by Dortenman, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Schüz, A.
Right arrow Articles by Dortenman, M.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2005. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oxfordjournals.org The online version of this article has been published under an open access model. Users are entitled to use, reproduce, disseminate, or display the open access version of this article for non-commercial purposes provided that: the original authorship is properly and fully attributed; the Journal and Oxford University Press are attributed as the original place of publication with the correct citation details given; if an article is subsequently reproduced or disseminated not in its entirety but only in part or as a derivative work this must be clearly indicated. For commercial re-use, please contact journals.permissions@oxfordjournals.org.

Quantitative Aspects of Corticocortical Connections: A Tracer Study in the Mouse

Almut Schüz, Denis Chaimow, Daniel Liewald and Monika Dortenman

Max-Planck-Institut für biologische Kybernetik, Spemannstrasse 38, 72076 Tübingen, Germany

Address correspondence to Almut Schüz, Max-Planck-Institut für biologische Kybernetik, PO Box 21 69, 72012 Tübingen, Germany. Email: almut.schuez{at}tuebingen.mpg.de.


    Abstract
 Top
 Notes
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Appendix
 References
 
This study provides neuroanatomical data relevant to models and simulations of the propagation of activity over the cortex. We administered small injections of the anterograde tracer biotinylated dextran amine to various regions of the mouse cortex (1 per animal). Two-dimensional reconstructions of the cortical surface were made, showing the distribution, size, and density of the terminal fields. Within the injected hemisphere, the largest part of the terminal field always surrounded the injection site and extended over neighboring areas. On average, axons from injection sites of ≤0.1 mm2 (containing several thousand neurons) diverged onto a region about 180 times larger than the injection site. The density of stained fibers in distant terminal fields could reach about 25 m/mm3. More than half of the total terminal field from an individual injection site consisted of weak projections with densities of 3 or 4 m/mm3. The number of main axons entering an individual distant terminal field ranged between 14 and about 890. By indirect arguments we estimate that the density of stained fibers close to the injection site is 3–6 times that in the most densely labeled distant terminal fields. In addition to symmetric projections to the opposite hemisphere, nonhomotopic callosal projections were found.

Key Words: anterograde • cortex • dextran amine • divergence • fiber density


    Introduction
 Top
 Notes
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Appendix
 References
 
One of the most outstanding properties of the cerebral cortex is its high degree of connectivity within itself. The great majority of the synapses in the cortex are synapses between cortical neurons (Cragg 1967Go; Braitenberg 1974Go; 1978aGo). Broadly speaking, there are 2 kinds of connections: local connections, which knit the cortical neurons together into one huge continuous network, and distant connections, which connect distant places directly to each other in a rather selective way (e.g., Jones and others 1978Go). All types of neurons seem to participate in the local connections via axon collaterals ramifying within the gray matter. Distant connections are made almost exclusively by pyramidal cells (in which we include the spiny stellate cells) either via the white matter or via long fibers within the gray matter. (For occasional reports on distant nonpyramidal connections, see Peters and others 1990Go; McDonald and Burkhalter 1993Go.) We will see that—at least in the mouse—the 2 systems of local and distant connections are not neatly separable.

The present study seeks to contribute to the understanding of the global connectivity of the cortex. Considerable research has already been made in this field, in days gone by, with the degeneration methods or with "physiological neuronography" (e.g., Bailey and others 1950Go) and, more recently, with the development of tracer methods. We therefore have a fairly good idea as to how the cortical areas are connected to each other. In the present study we shall not deal with the connectivity between particular areas but rather address the modelers who are concerned with the mechanisms of propagation of activity over the cortex (e.g., Abeles 1991Go; Kötter and Sommer 2000Go; Gerstein and Kirkland 2001Go; Hehl and others 2001Go; Knoblauch and Palm 2002aGo, 2002bGo; Mehring and others 2003Go; Robinson and others 2003Go; Izhikevich and others 2004Go). Therefore we focus as far as possible on quantitative aspects.

With the aid of small injections of an anterograde tracer, we seek to answer the following questions. 1) How large a part of the cortex can be reached by the neurons located in a small cortical "column"? (See Size of Terminal Fields in the Injected Hemisphere.) 2) What is the density of fibers in a distant terminal field? (See Fiber Densities.) 3) What proportion of the neurons in such a column projects to an individual distant projection site? (See Numbers of Neurons Projecting to a Distant Spot and Number of Neurons at the Injection Site.)

We used the anterograde tracer biotinylated dextran amine (BDA; Veenman and others 1992Go), and our experimental animal is the mouse. We chose this animal for 2 reasons. First, its small brain makes it easy to investigate the whole cortex for each tracer injection, and second, the results can be combined with other quantitative data gained in previous studies (summarized in Braitenberg and Schüz 1991Go).

Preliminary results have been published in Schüz and Liewald (2001)Go and in Chaimow and Schüz (2004)Go. A similar study using a retrograde tracer was made previously in our laboratory (Greilich 1984Go; Palm 1986Go).


    Materials and Methods
 Top
 Notes
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Appendix
 References
 
Experimental Animals and Tracer Injections

Adult white female laboratory mice (21 animals, at least 4 months old) were used for this study. Each animal received one iontophoretic injection of 10% BDA (BDA-10,000; Molecular Probes, Mo Bi Tec, Göttingen, Germany) in 0.01 M phosphate buffer, each in a different position of the neocortex, such that samples from anterior, posterior, dorsal, and lateral regions were obtained. The exact position of the injection was later determined on the histological sections and localized with reference to the map of Caviness (1975)Go. The injections were placed approximately at middle depth of the cortex that, depending on the particular area, corresponds to lower layer III to upper layer V.

For the injection procedure and the histology, we largely followed the protocols by Veenman and others (1992)Go and Dolleman-Van der Weel and others (1994)Go. (The deviations described below were probably not crucial for the success of the procedure.) The injections were performed under deep anesthesia, using a combination of ketamine hydrochloride and Rompun (xylazin hydrochlorid), intraperitoneally. The glass electrodes had outer tip-diameters of 20–30 µm. The tracer was applied in pulses of 5-µA positive current (7 s on, 7 s off) for 15 min. After a survival time of 7 days, the animals received a lethal intraperitoneal injection of 33% chloral hydrate and were first perfused with 0.9% NaCl, containing 0.05% heparin, and then with the fixative containing 4% paraformaldehyde in 0.1 M phosphate buffer (pH 7.4).

Histology

The brain was removed from the skull the next day and kept overnight in a 30% solution of sucrose in phosphate buffer on a rotator in the refrigerator. The brain was then cut with a freezing microtome into coronal sections of 50-µm thickness. Alternating sections were collected for BDA and Nissl stain in phosphate-buffered saline (PBS; the phosphate buffer was 0.01 M). In some cases where a second tracer (rhodamine-conjugated dextran amine) had been injected (not evaluated in this study), 3 alternating series were collected.

The Nissl stain was carried out with cresyl violet on mounted sections. The sections for the BDA stain were kept in PBS in the refrigerator for up to 4 days. They were then rinsed several times in cold PBS (3 x 10 min) and kept for 1 h in PBS containing 0.5% Triton X-100. For this and all the following steps, the sections were stored in the refrigerator on the rotator. BDA was visualized using the avidin–biotin complex method (ABC-kit PK-4000 Vectastain, Linaris, Wertheim-Bettingen, Germany) and diaminobenzidine (DAB). The sections were kept overnight in the ABC solution in darkness. They were then rinsed several times (2 x 10 min) in PBS and for 10 min in 0.1 M phosphate buffer. The DAB reaction was carried out according to the instructions of the substrate kit for peroxidase (Vector, SK-4100, Linaris, Wertheim-Bettingen, Germany) for about 15 min. The sections were then rinsed several times in phosphate buffer, mounted, air-dried, dehydrated, and covered in DePeX.

Shrinkage

We did not make systematic measurements on shrinkage of the tissue. However, to relate our measurements to data gained on different material, a short reflection on the shrinkage of the material used in the present study is in order.

There was no difference in size between adjacent BDA- and Nissl-stained sections. The shrinkage of Nissl-stained frozen sections was determined in an earlier article (Schüz and Palm 1989Go). There, we found a linear shrinkage of the sections to 88% in the 2 dimensions parallel to the plane of section. However, in contrast to our assumption in the previous study, this is not true of the third dimension, which a later study (Schüz A, unpublished data) showed to suffer a linear shrinkage to only 95%. The linear shrinkage factors relevant to the present study are therefore 0.88 in the 2 dimensions of the plane of section and 0.95 in the direction perpendicular to it. The volume shrinkage is then 0.74.

Two-Dimensional Surface Representations of the Terminal Fields

Two-dimensional maps of the cortical surface were reconstructed from 10 brains as follows: in a first step, the outlines of every second BDA-stained coronal section were drawn under a table projector at a magnification of x27. Terminal fields in the cortex stained densely enough to be recognized at this magnification, together with blood vessels and other landmarks, were marked on the drawings. In a second step, the sections were scrutinized under the microscope at higher magnification for less densely stained terminal fields. These were then also inserted on the drawings with the aid of the landmarks. In a third step, the length of the outline of the cortex as seen on the coronal section was measured between the corpus callosum (or otherwise defined medial landmarks rostral and caudal to the corpus callosum) and the ventrolateral margin of the cortex. This length was then plotted as a straight line, and the positions of the terminal fields were marked on it. The surface lines from the various sections were aligned side by side to produce the 2-dimensional maps showing the location and extent of the terminal fields (Figs. 4 and 5). Coronal sections through the very tip of the frontal and occipital pole (which contain only cortical gray matter) were not included. Thus, our maps are about 1 mm shorter in their rostrocaudal length than the real cortex is from pole to pole.


Figure 4
View larger version (32K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 4. Two-dimensional surface maps of the cortex (neocortex and part of the allocortex) of 6 mice. They are arranged according to the macroscopic location of the injection: central region of the dorsolateral convexity (ac), lateral (d), anterior (e), posterior (f). Right and left maps correspond to right and left hemispheres, respectively, rostral is at the top. The straight vertical line in each map corresponds to the border between dorsal and medial cortex, the elongated gray field medially corresponds to the corpus callosum (cc), the elongated gray field laterally is the olfactory tract (ot). The pyriform and entorhinal cortex are located laterally of the dotted line (the dotted line marks the dorsal border of the dark layer II of these allocortical areas, except for the most rostral part where it marks the location of the rhinal fissure). The medial piece caudally of the corpus callosum is the transition zone to the subiculum and hippocampus, both of which are not included on the maps. The white spot within the black field indicates the region of uptake. In (a) and (b), the injection was placed into the left hemisphere, in the other cases into the right hemisphere. The black region indicates the region of highest fiber density, and the regions of dense and wide hatching correspond to middle and low fiber densities, respectively.

 

Figure 5
View larger version (35K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 5. Two further maps, showing terminal fields as in Figure 4, but containing, in addition, measurements of 1) the number of pyramidal cells that have taken up the tracer (long arrows), 2) the number of stained fibers entering some of the terminal fields, and 3) maximum densities of stained fibers in these terminal fields.

 
Regions with fibers of passage showing no or only very few axonal boutons were not included on the maps.

Areal Measurements

On 8 of these maps (Figs. 4 and 5), the size of the terminal fields within the injected hemisphere and the total area of the neocortex of one hemisphere (but also including the cingulate cortex and the presubiculum) were measured. The measurements were carried out with a digital planimeter.

The Problem of Retrograde Staining

BDA is mainly, but not exclusively, an anterograde tracer. Individual retrogradely stained neurons are sometimes found in the sections, particularly in the vicinity of the injection site (e.g., Fig. 1). Some of them have a beautiful Golgi-like appearance, with extensive staining of axon collaterals. In most cases, these collaterals could be distinguished from antrogradely stained fibers emanating from the injection site, and so they did not affect the assessment of direct anterograde staining.


Figure 1
View larger version (86K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 1. Coronal section, showing the fiber stain close to an injection site. Numbers and arrowheads on the surface of the cortex illustrate the 3 different categories of fiber density used in this study. Long arrows indicate the width of the region of descending axons. The width of this region was taken as the diameter of the region of uptake (Table 1). A limited number of retrogradely stained neurons (upper left) can occur beyond this region. Bar: 100 µm.

 

View this table:
[in this window]
[in a new window]

 
Table 1 Areal measurements of terminal fields

 
However, we cannot exclude that some of these neurons had been stained via an axon collateral and that the tracer then moved on into the distant axonal tree. Such a retrograde/anterograde staining of a neuron located outside the region of uptake (via cell bodies and dendrites) may occasionally simulate an increased divergence if its distant axonal tree does not coincide with those of neurons located within the region of uptake. To minimize this problem, we took the following precautions: 1) the maps on Figures 4 and 5 were made from preparations containing only a few retrogradely stained neurons, 2) we excluded from the maps terminal fields consisting of very few fibers, and 3) The borders of weak projection zones were cut off in Figures 4 and 5 at a corresponding threshold.

Size of the Region of Uptake

The size of the injection, that is, the size of the region within which are located the cell bodies of the neurons that have taken up the tracer (via the cell bodies or the dendrites), was defined as the width of the band of descending axons (see Fig. 1).

Quantification of Fibers in the Terminal Fields

This was first performed in a semiquantitative fashion. The result is represented in Figures 4 and 5. The regions of highest density (category 1) were visible with the naked eye and are drawn black on the figures. Dense hatching (category 2) represents the regions that were visible under the low-power table projector. Those regions indicated by wide hatching (category 3) contained fibers that could only be recognized under the microscope. The arrowheads in Figures 1 and 3 indicate borders between these 3 categories. These categories (strong, middle, and weak) reflect the amount of stained fibers under a given surface area of cortex, in that dense hatching can either stand for a high density of fibers in one layer or a lower density over several layers. In some cases (Fig. 5b), it seemed appropriate to make 4 categories, where the 2 middle ones can be considered as subdivisions of category 2 in the other cases.


Figure 3
View larger version (97K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 3. Example for the appearance of the fiber stain on a coronal section through the telencephalon of the mouse, close to an injection site. Apart from the large terminal field close to the injection site in the right hemisphere, a distant terminal field can be seen in both hemispheres, as well as a terminal field in the thalamus (Th). The numbers and arrowheads on the surface of the cortex indicate the 2 different categories of fiber density that can be recognized at low magnification. Magnification of figure is x7.5.

 
In 2 of the brains, we made a more precise quantitative assessment of the terminal fields by actually measuring the fiber densities in them. Fiber densities can be measured as length densities. Length density is the total length of a filiform component in a given volume divided by that volume.

We used standard stereological methods to obtain length density estimates (e.g., Hennig 1963aGo, 1963bGo; summarized in Weibel 1979Go; Østerby and Gundersen 1988Go; Mayhew 1991Go; Russ and Dehoff 2000Go) and adapted them to our specific task. Circular test lines were superimposed onto the microscopic image (Fig. 2), and their intersections with the stained axons were counted. Length density estimates can be obtained from the number of intersections by the following formula:

Formula 1(1)
where LV is a length density estimate, N the number of intersections, l the length of the test line, and t the thickness of the section, in this instance 50 µm. A derivation of this formula can be found in the Appendix. The circular shape of the test lines guarantees that the sampling does not have any directional preference in the plane of section. For the above formula to hold, it had to be assumed that axonal orientations in any plane parallel to the cortical surface, that is, perpendicular to our vertical sections, are uniformly distributed. (For assumption-free approaches, see Cruiz-Orive and Howard 1990Go; Gokhale 1990Go; Larsen and others 1996Go; Mouton and others 2002Go.)


Figure 2
View larger version (162K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 2. Micrograph showing BDA-stained fibers in a terminal field and the circular test lines superimposed onto them.

 
The measurements were made with the aid of a computerized system (Eutectic Neuron Tracing System), connected to a microscope, using an objective lens x40/1. By focusing through the section, we could ensure that all intersections were recorded, regardless of the depth within the section. Four symmetrically arranged circular test line segments, each a quarter of a circle, were placed onto the densest spot of a region of termination (Fig. 2). Using the above formula, an estimate for the mean density was calculated.


    Results
 Top
 Notes
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Appendix
 References
 
Figure 3 shows a coronal section close to the site of injection in the right hemisphere. Axon collaterals radiate from the injection site within the gray matter, and distant terminal fields can be seen in the cortex of both hemispheres and in the thalamus. The axons are completely filled up to—and probably also within—their terminal fields.

Figures 4 and 5 show patterns of anterograde staining on maps of the flattened cortex as described in Methods. For each animal, both hemispheres are shown. In 6 animals, the injection was administered into the right, and in 2 animals into the left hemisphere. The black regions correspond to the highest density of staining, dense and wide hatching to progressively lower densities. The white spot in each figure shows the location of the injection site. Its size indicates the extent of the region of uptake as defined in Methods.

Injected Hemisphere

In spite of the large variation in projection patterns, some properties are common to the various injections: a large field of stained fibers around the injection and patches or bands at more distant places, numbering between 1 and 8 in the injected hemisphere. In the large field around the injection site, the overall density of stained fibers decreases toward the periphery, although denser patches sometimes appear within regions of lower density. The core of highest density (black on Figs. 4 and 5) extends about 500 µm (sometimes more) from the injection site and is due to the bulk of axon collaterals of the neurons affected by the injection. The adjacent hatched regions are fed by fibers running horizontally over long distances and by other fibers that take a bow-shaped course through layer VI, some of them scratching the white matter.

In distant terminal fields, individual axons can ramify such that they contribute fibers to several spots located closely to each other. This is, for example, the case in the 3 densely hatched lateral spots in the right hemisphere in Figure 5a.

Comparison with the Caviness Map

It is, of course desirable to have some idea of the cortical areas involved in Figures 4 and 5. Figure 6a shows the classical map in which Caviness (1975)Go has adapted the scheme of Brodmann (1909)Go to the mouse neocortex. On this map, we have marked the locations of the injection sites of the brains in Figures 4 and 5 (Fig. 6b). The localization was made by identifying the cortical areas on the Nissl sections according to the cytoarchitectonic criteria used by Caviness (1975)Go. The shape of the map differs from that of the neocortex in our maps by its larger rostrocaudal extent. This is due to the different treatment of the frontal and occipital pole regions: flattening by distortion in the Caviness map, cutting the tip of the poles, and keeping the rest metrically correct as far as possible in our maps. By means of the landmarks indicated in Figure 6b, it is possible to localize approximately some of the terminal fields.


Figure 6
View larger version (21K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 6. (a) Two-dimensional map of one hemisphere of the neocortex of the mouse (reprinted from Caviness 1975Go; with permission). It shows the relative position of the architectonic areas. (b) Location of the injections of Figures 4 and 5 on this map (large black dots), determined with the aid of coronal Nissl sections. The numbering corresponds to that of the corresponding maps in Figures 4 and 5. The Caviness map is longer than ours in its rostrocaudal extent. This is due to the 2-dimensional representation of the frontal and occipital poles in this map and corresponding distortions of surrounding regions. The following landmarks can be used for comparison with the maps in Figures 4 and 5: the lateral border of area 13 coincides with the dotted line in our maps. The broken line on the Caviness map corresponds to the straight vertical line in our maps. The dotted lines on the Caviness map serve to roughly indicate what parts are not contained on our maps. CC, corpus callosum; asterisks, assumed positions of the tip of the frontal and occipital pole in the Caviness map.

 
Some observations shall be mentioned. 1) Although the dark field around the injection site can be restricted to the corresponding cortical area, the large terminal field in continuity with it always extends over adjacent areas. 2) There is a striking similarity in the projection patterns in Figure 4b,c. This is obviously due to the fact that these injections were located in the same cytoarchitectonic area. Both injection sites project to a rostromedial region (area 6/8) and laterally to an elongated field comprising the border between area 35 and the adjacent entorhinal cortex and to a shorter elongated field of higher density, which according to the cytoarchitectonic criteria described by Caviness is at least partly contained in area 14. Both injections are located just caudally to the barrel field. 3) The 3 lateral spots in the right hemisphere in Figure 5a are located in different areas: areas 40 and 14. 4) No fiber stain could be found in the pyriform cortex, and in the entorhinal cortex fiber stain was restricted to its lateral border (e.g., Fig. 4b,c). Fiber stain was occasionally found in area 27, the presubiculum (medially posteriorly in Fig. 4f). 5) The small hatched field close to the olfactory tract in the right hemisphere of Figure 5b is located in the orbital region of the frontal pole (area 11).

Size of Terminal Fields in the Injected Hemisphere

Measurements on the size of the terminal fields in the injected hemisphere were made on the maps shown in Figures 4 and 5.

The field in continuity with the injection site was always larger than the sum of the remaining terminal fields (Table 1). In some cases, nearly all projections merged into one large field (Fig. 4e and 5b).

Table 1 shows the results of the measurements (not corrected for shrinkage). In column 1, the sum of the areas of all terminal fields in the injected hemisphere is reported for each of the 8 brains. Columns 2 and 3 show separate values for the areal size of the large terminal field in continuity with the injection site and the sum of the distant terminal fields. Columns 4–6 show separate values for the sums of the strong (i.e., black), middle, and weakly stained regions. The table also shows the approximate diameter of the region of uptake, as well as the cytoarchitectonic area in which the injection site was located.

The total area of the terminal fields ranges between about 7 and 22 mm2, with an average of 12 mm2. The difference in size can largely be explained by the inevitable variation in the strength of the injections and the corresponding variation in the size of the regions of uptake, with their diameters ranging between 200 and 400 µm (200 µm is the diameter of the dendritic extension of a large pyramidal cell in the mouse). The region of uptake ranged between about 0.03 and 0.1 mm2 in area.

As can be seen from Table 1, in all injections the main part of the terminal field consisted of weak projections (widely hatched regions in Figs. 4 and 5). If we exclude these, the size of the terminal fields (medium plus strong) ranges between 2.5 and 8.8 mm2, with an average of 4.7 mm2.

The average size of the surface area of the entire neocortex of one hemisphere was 64 mm2. Thus, the neurons from a middle-sized column of about 300 µm in diameter can reach about 7% of the total neocortex by projections of high plus middle strength, and even about 19% if the weak projections are included.

If we translate the measurements of Table 1 into degrees of divergence by dividing the area of the terminal fields by the area of the region of uptake, we get an average factor of 73 (with individual values between 45 and 122) for the strong plus middle projections and an average factor of 182 (with individual values between 122 and 298) if we include the weak connections. (The large range in the individual values may be partly due to the fact that the values for the size of the region of uptake are approximations only.)

Fiber Densities

In 2 of the brains (Fig. 5) we quantified the densities of stained fibers in distant projections. The contralateral projection in Figure 5b was one of the densest projections in our material (see micrograph in Fig. 7), whereas the rostral spot in the ipsilateral hemisphere was one of the weakest (see Fig. 8). We therefore took these 2 spots as instances of an upper and a lower limit. The densities of the 2 projections measured in Figure 5a were in between.


Figure 7
View larger version (175K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 7. Micrograph of the terminal field contralateral to the injection site in Figure 5b, as an illustration of a terminal field of high density.

 

Figure 8
View larger version (128K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 8. Micrograph of the terminal field of low density measured in the ipsilateral hemisphere in Figure 5b.

 
The results are shown in Figure 5 and Table 2. For comparison with other data (see Discussion), we express these densities in fiber length per cubic millimeters. Table 2 shows the raw data. The data in Figure 5 indicate the maximum densities measured in each spot because the maximum value may be relevant for the question of how strongly distant places in the cortex can influence each other functionally. These maximum values range between 4.6 and 21 m/mm3. The higher value may be an underestimate by some meters per cubic millimeters because in dense regions some intersections with the test line may escape the observer. Thus, the maximum density in the contralateral spot in Figure 5b may be rounded up to 25 m/mm3.


View this table:
[in this window]
[in a new window]

 
Table 2 Raw data of the measurements on the density of stained axons in 4 different terminal fields

 
Numbers of Neurons Projecting to a Distant Spot

We also counted the number of axons entering a terminal field from below. Assuming that each of these axons comes from a different pyramidal cell, this number is a measure of the number of neurons projecting from the injection site to that terminal field.

The results of these counts are also reported in Figures 5a,b. In Figure 5b, the number of axons entering a terminal field ranged between 14 for the weak spot to about 890 fibers running toward the dense and middle dense region of the callosal spot. The number for the contralateral spot in Figure 5a was about 230 and for one of the ipsilateral spots 540.

Number of Neurons at the Injection Site

We also made an estimate of the number of neurons that had taken up the tracer in the brains in Figure 5.

To a first approximation, one can calculate the number of neurons within the region of uptake. This calculation is made on the basis of the diameter of the region of uptake given in Table 1 and the number of neurons under 1 mm2 of cortical surface area, which is about 105 in the mouse (Rockel and others 1980Go; Schüz and Palm 1989Go).

However, not all the neurons within the region of uptake defined by the bundles of stained descending axons have necessarily taken up the tracer. In the core region, all neurons appeared to have taken up the tracer because the density of stained cell bodies there was as high as the known cell density in the mouse (Schüz and Palm 1989Go). Toward the periphery, the density of stained cells decreases relatively sharply. According to our estimates, an equivalent cylinder, containing as many neurons as stained in the core region and its halo of lesser staining, has a diameter about 60 µm lesser than the region of uptake defined by the stained descending axons.

With this reduced region of uptake, we end up with a number of about 9074 stained neurons in the brain in Figure 5b. Because we are concerned with the pyramidal cells that project over larger distances, we subtract from this number 15% nonpyramidal cells (Peters and others 1985Go; Braak H and Braak E 1986Go), which leads to an estimate of about 7700 pyramidal cells.

The injection in Figure 5a was smaller, and some of the cortical gray substance was mechanically destroyed by the injection. The number of stained pyramidal cells was about 2900 in this case.

Contralateral Projections

In all cases (with the exception of the brain in Fig. 4f) fiber stain could also be found in the contralateral hemisphere in a position symmetrical to the place of injection.

In some of the injections, additional stain could be found in other (i.e., nonhomotopic) positions on the contralateral side. When such projections were present in the contralateral hemisphere, terminal fields were always found in the corresponding position in the injected hemisphere. This raised the question as to whether these nonhomotopic callosal projections come directly from the injection site or rather are mediated by neurons stained retrogradely through one of their axonal branches and reaching the opposite side through another.

In some cases, marked by an asterisk in Figures 4 and 5, the latter possibility could not be excluded entirely because in these cases some retrogradely stained neurons were found on the side of the injection. Nor can we fully exclude that some of the weakly hatched regions around the homotopic contralateral spots are slightly enlarged due to fibers coming from retrogradely stained neurons in the vicinity of the injection site.

However, in other cases it was clear that nonhomotopic terminal fields were fed directly from the injection site. For example, in the brain in Figure 5a, some 230 fibers entered the contralateral terminal field symmetrical to the injection site, whereas at least twice as many crossed the corpus callosum. The rest of the stained fibers bypassed the homotopic site to reach the terminal field located more laterally.

Occasionally, an individual axon could be seen to ramify in the other hemisphere to reach 2 closely located contralateral terminal fields (e.g., brain in Fig. 4e) or to supply both the symmetrical contralateral field and its weakly stained surroundings (e.g., brain in Fig. 4b).


    Discussion
 Top
 Notes
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Appendix
 References
 
Methods

Completeness of the Fiber Stain

Wherever axons run in the plane of sectioning, they could be followed in their entirety. They showed rich ramifications upon reentering the cortex. This makes it probable that the tracer filled the axonal trees completely.

Aberrant Staining

When employing anterograde tracers one should be aware of the possibility of a certain amount of aberrant staining: transneuronal staining, retrograde staining, or staining of fibers passing the injection site. This raises the question of the extent to which the pictures in Figures 4 and 5 are partly the result of such aberrant staining.

We scrutinized our whole collection of BDA tracer preparations, which included deep cortical and subcortical injections, but found no evidence for transneuronal staining. However, evidence could be found for staining of fibers passing an injection site in both retrograde and anterograde directions. In our material, this phenomenon seemed to have been restricted to fibers injured by the electrode. The existence of a limited number of Golgi-like retrogradely stained neurons has already been mentioned, a phenomenon that may be activity dependent (Jiang and others 1993Go).

Both the fiber stain due to retrogradely stained neurons and the staining of injured fibers of passage are, however, responsible for a small percentage of stained fibers only. We can exclude that these artifacts had any effects on the size of the middle and densely stained terminal fields. Nor do we believe that they increased the size of the weak projection fields more than slightly because we had taken the precautions described in The Problem of Retrograde Staining.

Depth of the Injection

The location of the tip of the electrode varied between layers III/IV and V. This had some influence on the width of the region of uptake in the various layers, but in each case, neurons in all layers had taken up the tracer, as could be gathered from the distribution of darkly stained cell bodies.

Contralateral Projections

The existence of nonhomotopic callosal projections in rodents has been debated. For a review see Innocenti (1986)Go. The presence of such projections is evident in our material.

The complete lack of contralateral fiber stain in Figure 4f requires some explanation. This map shows an injection into area 17, the main part of which is known to be "acallosal" (also in the mouse; Yorke and Caviness 1975Go). In addition, as can be seen on the figure, this brain had an exceptionally small corpus callosum, possibly lacking some of the normal callosal fiber components.

Local versus Distant Connections

The axonal tree of a typical pyramidal cell is divided into local axon collaterals radiating in all directions within the gray matter and a main axon running through the white matter to a selected distant location. In the small brain of the mouse, this neat distinction into gray matter axon collaterals and white matter main axons can only be made with respect to the callosal and subcortical projections. In the large terminal field continuous with the injection site, the 2 systems cannot be separated because the main axons do not enter the white matter and can therefore often not be distinguished from long axon collaterals. Even in more distant projections within the same hemisphere, the main axons tend to be reluctant to enter the white matter.

In the following discussion, we shall therefore use the terms "local connections" for those fibers constituting the black field around the injection site in Figures 4 and 5 and "distant connections" for all those constituting the hatched regions. The term "distant spots" refers exclusively to the spots of terminal fields disconnected from the large field in continuity with the injection site.

Quantitative Aspects

Let us assume that a nest of neurons has been activated somewhere in the cortex, for example due to thalamic input. The activity will be relatively strong in the immediate vicinity of the activated neurons because all of them contribute axon collaterals to the local field. A particular distant spot will, however, only be reached by some of the activated neurons because only some of the pyramidal cells at a given place will project to the same distant spot. Most pyramidal cells probably have only one main axon projecting to a distant spot in the cortex. This assumption is supported by double-labeling experiments that show that callosally and intrahemispherically projecting neurons are largely separate populations (Andersen and others 1982Go; Schwartz and Goldman-Rakic 1982Go; however, for evidence of ramifying main axons in subcortically projecting pyramidal cells see Cowan and Wilson 1994Go; Deschênes and others 1994Go; Levesque and others 1996Go).

What then are the chances that the activity at a given place A transmitted to a particular distant spot B can activate neurons in B? We leave the final answer to the neurophysiologist or the theoretician modeling cortical activity but will deal with this question below as far as possible from an anatomical point of view.

Percentage of Pyramidal Cells Projecting from A to B

The percentage of pyramidal cells that project from a given injection site to some distant spot may be taken as an indication of the relative strength of distant connections. In Figure 5b the number of pyramidal cells projecting to the contralateral spot (ignoring the weak halo) is about 890. This is 12% of the about 7700 pyramidal cells that have taken up the tracer. In the case of Figure 5a, 540 out of 2900, that is, 19% of the stained pyramidal cells, project to the distant ipsilateral spot investigated. This means that 1/8 to 1/5 of the neurons from a given injection site can contribute fibers to a distant spot. Let us now investigate how this translates in fiber densities.

Density of Stained Fibers in Distant Spots

The maximum density of stained fibers measured in a distant spot in the brain in Figure 5b was about 25 m/mm3. This is very little compared with the total axonal length in 1 mm3, which is about 4 km in the mouse (Braitenberg and Schüz 1991Go; for comparable values in the cat see Foh and others 1973Go), that is, only about 0.6%. (If we take into account a volume shrinkage to 74% in the BDA preparations [see Methods], the maximum density reduces to 19 m/mm3, corresponding to 0.5% of the total axonal population.) In view of the overall arrangement of synapses along axons (Braitenberg 1978bGo; Amir and others 1993Go; Hellwig and others 1994Go; Houzel and others 1994Go; Anderson and others 2002Go), we take fiber density as a direct indication of synaptic density. Functionally speaking, this implies that of the 8000 postsynaptic sites on a dendritic tree of a pyramidal cell in the mouse (Braitenberg and Schüz 1991Go), about 48 synapses can be activated from a distant site of activation (in the event of a strong projection from that site).

Densities of Stained Fibers Close to the Injection Site

How do these values compare with the situation in the local terminal field in the vicinity of the injection site? This question cannot be answered by direct measurement because the density of stained fibers in the black regions in Figures 4 and 5 is too high to be measured by stereological methods.

However, the density of fibers in the black region can be estimated as follows. In former work we came to the conclusion that individual pyramidal cells in the mouse have a total axonal length between 20 and 40 mm within the gray matter, with an average closer to the upper than to the lower value (see Chapter 17 in Braitenberg and Schüz 1991Go). This enables us to calculate the total length of stained axons in the brains in Figure 5, where we have an estimate of the number of pyramidal cells that have taken up the tracer. In the case of Figure 5b, the estimated number of labeled pyramidal cells was 7700. Taking 35 mm as an average axonal length per neuron, we end up with a total length of stained fibers of 270 m. (In view of the uncertainties inherent in this calculation, we ignore the fact that some of these 7700 pyramidal cells will be output neurons to subcortical structures and therefore contribute less fibers to the cortical gray matter than the average pyramidal cell.)

We can now estimate how much of these 270 m is contained in the hatched subregions of both hemispheres. The remainder must then be contained in the black region. For this estimate, we must assign a fiber density and a volume to each of the hatched subregions. For the middle and widely hatched regions we used average values from Table 2 (12 m/mm3 and 3 m/mm3, respectively) and for the densely hatched regions we used the maximum value of the densest spot in Table 2 (21 m/mm3) instead of the average value (17 m/mm3) because—as already mentioned—the measurements in high-density regions tend to be underestimated by a few meters per cubic millimeter. To assign a volume to each subregion, we took the areal size of the region and the cortical depth within which fibers could be found in a given subregion. This depth varied between 100 and 600 µm in the various subregions.

For the brain in Figure 5b, we ended up with a total value of 81 m in all the hatched regions together. This leaves 270 – 81 = 189 m for the black region close to the injection site, that is, 70% of the total length of stained axons. The volume of the black region is 1.34 mm3 (the cortical depth at the injection site was 0.6 mm). The fiber density is then 189 m/1.34 mm3 = 141 m/mm3. This is 5.6 times higher than the maximum density of 25 m/mm3 in a distant field and corresponds to 3.5% of the total axonal density of 4 km/mm3.

We now perform the same kind of estimate with the brain in Figure 5a. There, about 2900 pyramidal cells had taken up the tracer, making a total fiber length of 2900 x 0.035 = 102 m. In the hatched regions, we end up with a sum of 46 m. This leaves 102 – 46 m = 56 m for the dark region around the injection site (including the disconnected black spot). Thus, 55% of the total length of stained fibers is located within the dark field. The volume of the black region is 1.47 mm3 (the cortical thickness is 1 mm here) and the fiber density is then 56 m/1.47 mm3 = 38 m/mm3. This is 2.7 times more than the highest distant fiber density measured in this brain (see Table 2).

Axonal Length per Neuron

The above calculations on fiber density can, in turn, also be used to provide a lower bound for the possible average axonal length per neuron within the range of 20–40 mm mentioned above. If we had made the above calculation with an axonal length of 23 mm instead of 35 mm, we would have ended up with a fiber density in the local field of 14 m/mm3 in the brain in Figure 5a, that is, equal to that in the densest distant field in that brain. Because the stain is always strongest close to the injection site, an average axonal length of 23 mm is definitely too low. (For comparison with axonal length of individual neurons in the cat see Tettoni and others 1998Go; Binzegger and others 2005Go.)

If we go to the upper extreme of the range given above, we get an estimate of the upper bound of the percentage of stained axons in the black region. With an average axonal length per neuron of 40 mm, we get a total length of stained axons of 308 m in the brain in Figure 5b and 116 m in the brain in Figure 5a. The percentage of stained pyramidal cell axons in the black region is then 74% and 60%, respectively.

To Summarize These Considerations on Fiber Densities

The density of stained fibers in the near field surrounding the injection is always higher—usually much higher—than in a distant spot. However, strong distant projections can reach densities that are one-sixth or even one-third of the local density. High density of fiber stain in the close vicinity of an injection site as compared with more distant patches has also been found in other species (Malach and others 1993Go; Bosking and others 1997Go).

The finding that 55–70% of the stained axons were located in the near field provides an answer to the longstanding question of the relative importance of local versus distant projections of the axonal tree of a pyramidal cell. The data suggest that the local axonal tree is at least as large and probably somewhat larger than the parts projecting to the hatched regions.

Global Aspects of Connectivity

The high degree of divergence of fibers from a small column supports the view of the cortex as a large mixing machine (Braitenberg 1986Go). The axons originating from a piece of cortex 0.1 mm2 in size or smaller were distributed in a region measuring about 70 times that area, and even about 180 times that area, if weak projections are included. Injections of comparable size in the cat show an even larger degree of divergence (Henry and others 1990Go). (For further comparison with larger species, see Burkhalter and Charles 1990Go; Amir and others 1993Go; Levitt and others 1993Go; Lund and others 1993Go; Pucak and others 1996Go; Porter 1997Go; Barone and others 2000Go; Levitt and Lund 2002Go; Kondo and others 2003Go.)

The low density of projections in more than half of the total terminal field from an injection site needs some discussion. One wonders what contribution they can make to the firing of neurons at a distant place. However, weak projections make sense in connection with the theory of cell assemblies (Hebb 1949Go), the theory of synfire chains (Abeles 1982Go), and other similar ones that rely on cooperative effects in large numbers of neurons. It should also be mentioned that under certain conditions the membrane potential of cortical neurons can be poised just below threshold for firing (Metherate and Ashe 1993Go; Cowan and Wilson 1994Go; Stern and others 1997Go; Timofeev and others 2001; Destexhe and others 2003Go; Léger and others 2005Go). In this case, weak connections could have a decisive effect.

The large size of the terminal field in continuity with the injection site, extending over neighboring areas, suggests that most of the connectivity between areas is devoted to closely located areas. This seems to be true for small and large brains alike because it is in accordance with tracer studies in cats and monkeys (summarized in Young and others 1995Go) and with the distribution of fiber lengths in the cortical white matter in humans (Schüz and Braitenberg 2002Go).


    Appendix
 Top
 Notes
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Appendix
 References
 
Here we will give a derivation of the length density formula (formula (1) in Methods) based on our specific task. We want to estimate length density, that is, the total length of a filiform structure in a given volume divided by that volume. The number of intersections between line segments with random position and orientation and a superimposed test line is proportional to the length density.

In our microscopic image, we can see only the 2-dimensional projection of these line segments. The first thing we therefore need to know is how the projected line length is related to the actual line length in space. As mentioned in Methods, we assume that all axonal directions in a plane parallel to the cortical surface are uniformly distributed. To our knowledge, this is a reasonable assumption.

The projected lengths p are shortened relative to the real lengths r by the factor cos {alpha}, where {alpha} is the angle between the fiber and the plane of section:

Formula 2(2)
Given a uniform distribution of 3-dimensional orientations, we are more likely to find small values for {alpha} than large ones (see Fig. 9). The probability density for finding the angle {alpha} is equal to its cosine.

Formula 3(3)
The expected shortening is the integral of p/r over all possible angles {alpha} weighted with P({alpha}).

Formula 4(4)
With this result, we can now reduce our original problem to the estimation of the length density in 2 dimensions.


Figure 9
View larger version (11K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 9. A 3-dimensional orientation can be represented by a vector v on a sphere that is parametrized by an angle of elevation {alpha} and an azimuthal angle {theta}. For a uniform distribution of 3-dimensional orientations, the probability density for finding a vector v is equal at any place on the sphere. Given a specific elevation {alpha}, v is constrained to lie on the circle c, which is the horizontal section of the sphere at the elevation {alpha}. The probability density for any vector to lie on that circle is proportional to its circumference, which is proportional to cos {alpha}. Therefore the probability density of finding a vector v with an elevation {alpha} in a uniform distribution of orientations is proportional to cos {alpha}.

 
Estimation of length density in 2 dimensions is a classic problem in stereology (Buffon 1777Go; Smith and Guttman 1953Go). The main idea is that the number of intersections between a linear structure such as the projected axons in our case and a test line in a given area depends on how much total length is distributed in that area. The formula for length density in 2 dimensions is

Formula 5(5)
N is the number of intersections and l is the length of the test line. This formula assumes isotropic and uniform random sampling. One easy way of acquiring isotropicity is by using circular test lines.

We will now give an intuitive explanation of the above formula. Imagine a set of parallel lines of length b that is crossed by a perpendicular test line of length lstraight (Fig. 10a). In this case, the number of intersections N with that test line tells us exactly how much line length L is contained in the shaded area A.

Formula 6(6)
However, we require a formula that works for arbitrary orientations. This is achieved by replacing the straight test line by a circular one (Fig. 10b). The number of intersections remains identical, whereas the length of the line is now {pi}/2 times longer. This leads to the following formula:

Formula 7(7)
Let us now combine formulae (4) and (7). Formula (7) gives us the total length of axonal segments per area, ignoring their shortening due to their projection onto the plane of section. According to formula (4), the average shortening factor is {pi}/4. If we incorporate the section thickness t = 50 µm, we arrive at formula (1) of Methods:

Formula 8(8)


Figure 10
View larger version (24K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 10. A set of parallel lines running through a rectangular area (gray) has a length density that equals the number of lines times the length b of a single line segment contained in the shaded area divided by that area. (A) A perpendicular test line of length lstraight is used to estimate the length density based on the number of intersections. (B) A circular test line is unbiased to any orientation. The number of intersections is the same as in A, but the length of the test line has changed.

 

    Notes
 Top
 Notes
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Appendix
 References
 
Funding to pay the Open Access publication charges for this article was provided by the Max-Planck-Institute for Biological Cybernetics, Tübingen, Germany.


    Acknowledgments
 
Valentino Braitenberg participated throughout in the discussions regarding this work. We are grateful for his precious advice. We also thank Günther Palm, Fritz Sommer, and Robert Miller for discussions and suggestions. Special thanks to Günther Palm for preventing us from neglecting the point made by formula (3). We also thank Michaela Schweizer for technical advice and Lea Arndt for carrying out some of the measurements and drawings. In a later phase, discussions with Ad Aertsen, Stefan Rotter, and Nicole Voges strongly influenced this work. We are grateful also to Shirley Würth, who improved our English.


    References
 Top
 Notes
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Appendix
 References
 
Abeles M. 1982. Local cortical circuits. Berlin: Springer-Verlag.

Abeles M. 1991. Corticonics. Neural circuits of the cerebral cortex. New York: Cambridge University Press.

Amir Y, Harel M, Malach R. 1993. Cortical hierarchy reflected in the organization of intrinsic connections in macaque monkey visual cortex. J Comp Neurol 334:19–46.[CrossRef][ISI][Medline]

Andersen RA, Asanuma C, Cowan WM. 1982. Observations on the callosal and associational cortico-cortical connections of area 7a of the macaque monkey. Soc Neurosci Abstr 8:210.

Anderson JC, Binzegger T, Douglas RJ, Martin KAC. 2002. Chance or design? Some specific considerations concerning synaptic boutons in cat visual cortex. J Neurocytol 31:211–229.[CrossRef][ISI][Medline]

Bailey P, Bonin GV, McColluch WS. 1950. The isocortex of the chimpanzee. Urbana, IL: University of Illinois Press.

Barone P, Batardiere A, Knobloch K, Kennedy H. 2000. Laminar distribution of neurons in extrastriate areas projecting to visual areas V1 and V4 correlates with the hierarchical rank and indicates the operation of a distance rule. J Neurosci 20(9):3263–3281.[Abstract/Free Full Text]

Binzegger T, Douglas RJ, Martin KAC. 2005. Axons in the cat visual cortex are topologically self-similar. Cereb Cortex 15:152–165.[Abstract/Free Full Text]

Bosking WH, Zhang Y, Schofield B, Fitzpatrick D. 1997. Orientation selectivity and the arrangement of horizontal connections in tree shrew striate cortex. J Neurosci 17(6):2112–2127.[Abstract/Free Full Text]

Braak H, Braak E. 1986. Ratio of pyramidal cells versus non-pyramidal cells in the human frontal isocortex and changes in ratio with ageing and Alzheimer's disease. In: Swaab DF, Fliers E, Mirmiran M, Van Gool WA, Van Haaren F, editors. Progress in brain research. Amsterdam: Elsevier. p 185–212.

Braitenberg V. 1974. Thoughts on the cerebral cortex. J Theor Biol 46:421–447.

Braitenberg V. 1978a. Cortical architectonics: general and areal. In: Brazier MAB, Petsche H, editors. Architectonics of the cerebral cortex. New York: Raven Press. p 443–465.

Braitenberg V. 1978b. Cell assemblies in the cerebral cortex. In: Heim R and Palm G, editors. Lecture notes in biomathematics. Vol. 21. Berlin: Springer. p 171–188.

Braitenberg V. 1986. Two views of the cerebral cortex. In: Palm G, Aertsen A, editors. Brain theory. Berlin, Heidelberg: Springer-Verlag, p 81–96.

Braitenberg V, Schüz A. 1991. Anatomy of the cortex—statistics and geometry. Berlin: Springer-Verlag. [Revised edition (1998) under the title: Cortex: statistics and geometry of neuronal connectivity].

Brodmann K. 1909. Vergleichende Lokalisationslehre der Großhirnrinde. Leipzig: Verlag Johann Ambrosius Barth.

Buffon GLL de, 1777. Essai d'arithméthique morale. Suppl. à l'histoire Naturelle. Vol. 4. Paris: Imprimerie Royale. p 47–148.

Burkhalter A, Charles V. 1990. Organization of local axon collaterals of efferent projection neurons in rat visual cortex. J Comp Neurol 302:920–934.[CrossRef][ISI][Medline]

Caviness VS. 1975. Architectonic map of neocortex of the normal mouse. J Comp Neurol 164:247–264.[CrossRef][ISI][Medline]

Chaimow D, Schüz A. 2004. Quantification of anterogradely stained axons in the cerebral cortex. In: Bülthoff HH, Mallot HA, Ulrich R, Wichmann FA, editors. Proceedings of the 7th Tübingen Perception Conference. Kirchentellinsfurt: Knirsch-Verlag. 100 p.

Cowan RL, Wilson CJ. 1994. Spontaneous firing patterns and axonal projections of single corticostriatal neurons in the rat medial agranular cortex. J Neurophysiol 71:17–32.[Abstract/Free Full Text]

Cragg BG. 1967. The density of synapses and neurones in the motor and visual areas of the cerebral cortex. J Anat 101(4):639–654.[ISI][Medline]

Cruiz-Orive LM, Howard CV. 1990. Estimating the length of a bounded curve in three dimensions using total vertical projections. J Microsc 163(Pt 1):101–113.