Abstract
Barrelless mice (BRL) homozygous for the BRL mutation that disrupts the gene coding for adenylyl cyclase type I on chromosome 11 lack spatial segregation of layer IV cortical cells and of the thalamocortical axons (TCAs) into barrel domains. Despite these morphological perturbations, a functional topographic map has been demonstrated. We reconstructed individual biocytin-injected TCAs from thalamus to barrel cortex in NOR (normal) and BRL mice to analyze to what extent the TCA arborization pattern and bouton distribution could explain the topographic representation of the whisker follicles. In BRL, the geometry of TCA is modified within layer IV as well as in infragranular layers. However, in both strains, the spatial distribution of TCA in layer IV reflects the spatial relationship of their cell bodies in the ventrobasal nucleus of the thalamus. The morphometric analysis revealed that TCAs of both strains have the same length, branch number, and number of axonal boutons in layer IV. However, in barrelless, the boutons are distributed within a larger tangential extent. Analysis of the distribution of boutons from neighboring thalamic neurons demonstrated the existence in layer IV of domains of high bouton density that in both strains equal the size and shape of individual barrels. We propose that the domains of high bouton density are at the basis of the whisker map in barrelless mice.
Introduction
In direct sensory pathways of the mammalian brain, axons projecting from thalamus to the primary sensory cortex have a stereotypic morphology. In the thalamus, they send collaterals to the reticular nucleus and subsequently ascend through the internal capsule to terminate in cortical layers VI and IV (Jones, 1985, 1986). Subcortically, the topological organization of thalamocortical axons (TCAs) undergoes distinct transformations (Adams et al., 1997). In cortex, the spatial segregation of TCAs underlies the functional parcellation (Blasdel and Lund, 1983; Humphrey et al., 1985; Jensen and Killackey, 1987a; Antonini et al., 1998, 1999). During development, sensory activity is an important factor in determining the spatial distribution of TCAs as was demonstrated after visual deprivation in cats (Antonini and Stryker, 1993) and whisker follicle denervation in rats (Jensen and Killackey, 1987b). Detailed studies of thalamic projections in the whisker-to-barrel pathway revealed that modifying neuronal transmission in cortex influences the spatial configuration of TCA during development. Blocking the activity through the NMDA receptor in the conditional NR1 knock-out affects TCA morphology but does not modify the cytoarchitectural differentiation into barrels (Iwasato et al., 2000; Lee et al., 2005). Barrel formation is perturbed in several mouse lines: the knock-out for the metabotropic glutamate receptor 5 (mGluR5) receptor (Hannan et al., 2001), the monoamine oxidase A (MAOA)-deficient mouse (Cases et al., 1996), the GAP-43 (growth-associated protein 43) (Maier et al., 1999) and in the mouse mutant barrelless, which lacks functional adenylyl cyclase I (Welker et al., 1996; Abdel-Majid et al., 1998; Lu et al., 2003). Whereas a morphological analysis of the TCA axons is not available for the mGluR5 knock-out, in the other mouse lines lack of barrel formation is attributable to a disrupted clustering of TCA within barrel domains (Welker et al., 1996; Rebsam et al., 2002; McIlvain et al., 2003). Despite this altered arborization pattern of TCAs, a functional topographic map was demonstrated in SI of both MAOA and barrelless mice (Welker et al., 1996; Yang et al., 2001).
The aims of the current study are the following: (1) to characterize the arborization pattern of TCA in wild-type and barrelless mice to identify the level at which the mutation affected axonal morphology, and (2) to analyze whether the tangential distribution of the TCA boutons in barrelless mice can explain the establishment of a functional map of the sensory periphery. We reconstructed individual biocytin-labeled TCAs from thalamus to cortex and registered the spatial distribution of the axonal boutons that correspond to TCA synapses (White et al., 2004). We demonstrated that TCAs in barrelless are identical to those in wild type with respect to total axonal length and number of axonal boutons, but that the mutation affects the tangential extent of the axonal plexuses in layers IV and VI, confirming our previous report (Welker et al., 1996). Superposition of the distribution of axonal boutons from TCAs with cell bodies that are adjacent in the thalamus showed that within the cortex tangential domains are created with high density of axonal boutons that equal the size of individual barrels.
Materials and Methods
Animals
The barrelless mice (Welker et al., 1996) used in this study are from a breeding line that was started after the discovery of a mutation that had occurred in the normal (NOR) strain (i.e., mice bred for a standard pattern of mystacial vibrissas) (Van der Loos et al., 1984). Barrelless mice are homozygous for the BRL mutation that disrupted the gene coding for adenylyl cyclase type I on chromosome 11 (Welker et al., 1996; Abdel-Majid et al., 1998). In total, six NOR and 11 barrelless animals were used in this study. Mice of both strains were exposed to the same experimental procedures (see below) that were approved by the Office Vétérinaire Cantonal (Lausanne, Switzerland), in accordance with Swiss laws.
Tracing and reconstruction of thalamocortical axons
Adult females (age, 6–10 weeks) were anesthetized with Nembutal (60 mg of sodium pentobarbital/kg of body weight, i.p.) and placed in a stereotaxic frame. Glass pipettes filled with 2% biocytin in 1 m potassium acetate were lowered into the ventrobasal nucleus of the thalamus (VB) of the right hemisphere. Stimulation-dependent neuronal activity was recorded through the pipette to identify the representation of the mystacial whisker follicles. Subsequently, the tracer was iontophoretically injected using a 150–200 nA positive current for 10–15 min (1 s on, 1 s off) according to the protocol of Pinault (1996). In most cases, the injection was made in the representation of follicle C2. After 24 h, mice were reanesthetized with a lethal dose of Nembutal and were fixed via transcardial perfusion of 4% paraformaldehyde in phosphate buffer (0.1 m; pH 7.2). Brains were removed and postfixed for 2 h in same fixative after which they were placed in a 30% solution of sucrose overnight. Brains were cut coronally on a freezing microtome at 50 μm, and the sections were processed for the revelation of biocytin using the avidin–biotin reaction kit from Vectra stain. Sections were Nissl stained with methylene blue. For additional details on the histological procedure, see Dolleman-Van der Weel et al. (1994). Cortical areas were assigned using cytoarchitectonic criteria as described by Caviness (1975).
The number of labeled thalamic neurons was counted. A total of 35 NOR and 17 BRL axons were hand-drawn using an ordinary light microscope fitted with a drawing tube. These drawings were used to study the axonal trajectory and the principal axonal bifurcations in the somatosensory cortex. From these sets of axons, eight NOR axons (from four mice) and nine BRL axons (from three mice) were randomly assigned for three-dimensional reconstruction using a Zeiss (Oberkochen, Germany) microscope equipped with a video camera (930XP 3CCD Color Video; Sony, Tokyo, Japan), a microcomputer-controlled stage, and a reconstruction program, Neurolucida (MicroBrightField, Williston, VT). Figure 1 gives the positions of the injection sites of these reconstructed axons.
Per mouse, two to six axons were reconstructed simultaneously from VB to the cortex over 10–25 sections. In the cortex, the entire axon was reconstructed, and the branch points (nodes) and endings were plotted. Axonal branches represent the axonal segments between two consecutive nodes. The branch was “incomplete” if, for any reason, its reconstruction was not possible. Branches with a length ≤4 μm were not traced but their position was identified. The position of the axonal boutons was also plotted. A correction (2.66×) was applied for the shrinkage in the z-axis because of the histological procedure.
The NEUROLUCIDA data file was transferred to a UNIX station (Silicon Graphics, Mountain View, CA) for additional analysis using MAXSIM software (Tettoni et al., 1996).
Measurements
MAXSIM software enabled measurements of the entire axonal intracortical arborization and its components: the trunks and the tangential plexuses. Trunks were defined as parts of the axon that traversed the cortical layers and gave off branches in the two target layers: the upper layer VI and layer IV. The tangential plexus of layers I–IV are all the branches coming off by the axonal trunks in layers IV, III, II, and I. The tangential plexus of layers V–VI is formed by the branches that course tangentially at the white matter (WM)/cortical border and at the VI/V interface. The following morphological parameters were determined.
The “neurometric” parameters.
To characterize the axonal tree and the distribution of boutons, we measured the axonal length, number of axonal branches, and the position of each axonal bouton. Measurements of these parameters were made first on the entire intracortical arborization by selecting the axon at its cortical entrance (Fig. 2A, arrow). Using the Strahler centripetal ordering (Fig. 2C), the program further allows us to topologically characterize the axon. Two compartments were distinguished: a “transmission” compartment formed by the Strahler's orders 0 and 1, and a “conduction” compartment formed by the sum of the orders 2, 3, and 4 (Fig. 2B,C).
To obtain the parameters of the layer IV tangential plexus, the entrance of the trunk in this layer was identified (Fig. 2E, green crosses). For the V–VI plexus, the tangential oriented infragranular branches (Fig. 2E, blue) were selected at the emergence of their trunk. The points at which the axonal trunks entered layer IV were further used to determine the “entrance surface” of the axon in layer IV (see below).
The “surface” parameters.
These parameters are two-dimensional and represent the areal extent of the axon projection in a plane tangential to the pial surface overlying the representation of the large mystacial whisker follicles (Woolsey and Van der Loos, 1970). In this plane, we determined the entrance surface of the axon in layer IV (Fig. 2F) and analyzed the distribution of axonal boutons of the tangential plexuses in layer IV and layer VI using a grid consisting of 20 × 20 μm squares. In each square, the density of boutons was calculated. The total surface area is the sum of squares in which one or more boutons were present. For additional analysis, we used five density ranges: (0,1], (1,2], (2,4], (4,8], and (8,64].
Statistics
Statistical analyses were performed using the SAS software package (SAS Institute, Cary, NC) and according to SAS/STAT user's guide, version 6. The distribution of parameters was tested for normality, and rank transformations were used if necessary. The differences were considered to be statistically significant at p ≤ 0.05.
The morphometric parameters were compared using a multivariate ANOVA (SAS/STAT user's guide, version 6; GLM procedure). A factor analysis (FACTOR procedure) was used to obtain a global image of the interdependencies between the parameters. Only factors with an eigen value ≥1 were considered to well characterize the variability inside a tangential plexus. A regression model was used (GLM and REG procedures) to study the hypothesis of homogeneity of slopes between strains for each layer and between layers inside each strain.
A χ2 test (FREQ procedure) was used to tests the differences between the frequency distribution of the number of terminals reached after a given number of nodes in NOR and BRL.
Results
Neuronal labeling and subcortical trajectory
Subcortical axonal trajectory is common in NOR and BRL
Biocytin injections were made in the VB of the thalamus where the large mystacial whisker follicles are represented (Fig. 1A,B,D). For both NOR and BRL, injections situated more dorsally resulted in labeled cells whose axons gave arborization in SI, whereas those situated more ventrally, gave arborizations in SI and SII. The labeled axons were seen to arise directly from the soma or from a dendrite close to its emergence from the soma.
The study of the hand-drawn TCAs from NOR (n = 35) and BRL (n = 17) showed that the subcortical trajectory is very similar in both strains (Fig. 3). Inside VB, the TCA change direction several times at variable distances from the cell body before arriving at the VB border. TCAs cross the thalamic reticular nucleus (RTN) parallel to each other. In this nucleus, the axons give off collaterals (1–4 per axon; mean number, 1.9) that form a dense plexus oriented perpendicular to the trajectory of the labeled axons. Leaving the RTN, TCAs enter the internal capsule in which axons follow diverging trajectories through the striatum. Figure 3 illustrates this divergence both in a coronal and sagittal plane. It shows that axons from neurons that are close neighbors in VB do not maintain this proximity within their trajectory to the somatosensory cortex. On reaching the subcortical WM, some TCAs run for long distances at the border with layer VI; other axons enter the cortex and travel tangentially for some distance in layer VI. The cortical entrance and the eventual axonal branching at the WM/VI interface are, in all cases, within the boundaries of the somatosensory cortex.
Using the set of computer-reconstructed axons, we found that the subcortical trajectory of TCA of NOR axons has a mean length of 3.4 ± 0.4 mm and is not different from that of BRL axons (3.7 ± 0.4 mm).
Progressive convergence in NOR
In NOR, the majority (32 of 35) of labeled axons restrict their cortical entrance to the territory of three neighboring columns, including the “target” barrel. However, some axons (3 of 35) enter SI at its border with SII, with occipital (Fig. 3, NOR, sagittal view, arrow) or with motor cortex. These axons curve at 90 or 180° and travel inside layer VI to reach the column of their target barrel.
Of the 35 axons, only one enters the cortex within the appropriate barrel column. Of the 34 axons, 19 converge (with one or more of their trunks) to the appropriate column within layer VI; 10 others direct their trunks toward the appropriate barrel column within layer V or at its border with layer IV. One axon corrects the trajectory of one of its trunks at the level of layer IV, whereas two axons use the supragranular layers for turning into the appropriate barrel column. The trajectories of these 33 axons illustrate the trend to progressively converge toward the appropriate barrel column from layer VI upward. Two axons do not converge with their trunks into a single barrel column; their trunks enter two barrels, where they give raise to an extensive terminal plexus in one barrel and give off some terminal branches in the other.
In layer VI and in layer IV, NOR axons give off branches that give rise to bouton-rich plexus (Fig. 4A–C). In general, the infragranular plexus extends into several barrel columns, whereas the plexus in layer IV is highly concentrated in a single barrel. Two of the 35 analyzed axons terminated in more than one barrel; one of these is illustrated in Figure 4C. The layer VI plexus is essentially formed by branches that run tangentially in upper VI and that originate from radial trunks that ascended to layer IV (Fig. 4A,B). In two cases, we identified branches that arise from nodes at the WM/VI border climb in upper VI where they contribute to the infragranular plexus (Fig. 4B).
Progressive divergence in BRL
At the cortical level, BRL axons are highly divergent (Fig. 3) and do not enter SI within a three barrel column area as described above for NOR. Of the 17 hand-drawn BRL axons, 12 enter SI with several, highly divergent, branches (Fig. 5A) that climb up to layer IV. In addition to these long, bouton-bearing branches, we also notice a large number of very short axonal branches (<4 μm) that are formed by the axon at the level of the border between the white matter and layer VI (Fig. 1F). However, four axons enter SI with a single trunk, but bifurcate in the infragranular layers and give off branches that run ∼1.5 mm tangentially inside layer Vb, before entering layer IV (Fig. 5B). One axon approaches layer IV with a single axonal trunk that bifurcates just before entering (Fig. 5C).
The number of branches by which a TCA enters layer IV in barrelless is twice that of NOR [BRL, 5.4 ± 2.2 (SD); NOR, 3.3 ± 1.4 (SD); p < 0.005]. In addition, the surface of entrance of these branches is 31 times larger than in NOR (Table 1, p = 0.001). Several of these branches are given off at the border between white matter and layer VI or during the intracortical trajectory of the axon. The result of this branching pattern is a progressive divergence of the axonal arborization while reaching layer IV, which is in contrast with the increase convergence by which TCA approaches layer IV in NOR.
As in NOR, TCA in barrelless give rise of axonal bouton-rich plexuses in layers IV and VI. However, inside the layer IV plexus in barrelless, axonal branches seem not to encounter a spatial restriction, which adds another factor to the total divergence of the axonal arbor.
Despite the divergence, in their cortical termination, BRL axons conserve their topographic relationship as found in the thalamus (Fig. 3, BRL). An example is depicted in Figure 6. In this BRL mouse, the biocytin injection was performed more ventrally in the thalamus and TCA made branches in both SI and SII. The biocytin-labeled cell bodies were close together in VB. One axon, instead of reaching SI by traveling at the WM/VI interface, entered in the cortex at the SI/SII border and gave branches in layer VI of the SII. It subsequently traveled in layer III for ∼1500 μm to reach the same region in SI as its VB neighbor.
Quantitative analysis of individual axons
Analysis of the intracortical arborizations
Table 1 lists the quantitative parameters used to analyze the differences between NOR and BRL axons. On average, the intracortical arborizations of TCA in NOR mice have a length (L) of 17.0 mm (SD, 8.1 mm), carry 1234 (SD, 836) boutons (B), and are formed by 209 (SD, 196) branches (BR). In the granular and supragranular layers, these axons spread over a mean tangential surface (S4) of 0.083 mm2 (SD, 0.026), the intracortical arborization of BRL axons have a tendency to be longer (21.5 ± 5.6 mm) than those in NOR, but possesses the same mean number of boutons (1336 ± 590) and mean number of branches (196 ± 41). The difference in length is mainly attributable to a 40% increase of the conduction compartment in BRL. The tangential surface S4 in BRL is 0.17 mm2 (SD, 0.055 mm2), twice that in NOR (p < 0.002).
The differences between NOR and BRL axons were statistically analyzed using a one-way multivariate analysis and a factor analysis on the combination of the following five parameters: B represents the total number of boutons of the entire intracortical arbor (Fig. 2A, the yellow portion of the axon); LT represents length of the transmission compartment of the cortical arborization (all the branches of order 0 and 1 found mostly in the layers I–IV and V–VI) (Fig. 2B, the red portion of the axon); LC is the length of the conduction compartment of the intracortical arborization (the axonal trunks plus all the branches of order 2, 3, and 4 inside the tangential plexuses) (Fig. 2B, blue); SE is the entrance surface of the axonal trunks in layer IV (Fig. 2F); S4 is the tangential surface of the axon inside layers I–IV.
LC and SE parameters characterize the part of the axon by which it reaches layer IV. The parameters V, S4, and LT characterize the tangential plexuses where most of the boutons and the 0 and 1 order branches are found.
A one-way multivariate analysis on the combination of the five parameters LT, LC, B, S4, and SE demonstrated an important “strain effect” (p = 0.009). To analyze the interdependencies of the five parameters, a factor analysis was performed on NOR and BRL axons as a single sample. It demonstrated that in both strains individual TCAs are characterized by (only) two factors. Figure 7A indicates how the five parameters correlate to these two factors in NOR and in BRL. For NOR axons, a single factor (“factor 1”) was sufficient to explain 72% of the variability with an eigenvalue of 3.63. This factor comprises four of the parameters, except SE that was found to be part of second factor (“factor 2”) responsible for only 18% of the variability and which had no statistical significance (eigenvalue <1). In BRL axons, factor 1 covered only 56% of the variability and the analysis identified a significant contribution of factor 2 comprising 27% of the variability with an eigenvalue of 1.38. For these axons, factor 1 highly correlated with B, LT, and S4, and factor 2 correlated with LC and SE. This multivariate analysis therefore indicates that the strain difference between TCA of NOR and BRL is expressed by two parameters, LC and SE, whereas other parameters seem to have been unaffected by the mutation. Figure 7B is a graphical representation of the distribution of these two factors and demonstrates how factor 2 segregates the group of BRL axons from those of NOR.
Because LC represents the length of the conduction compartment, we analyzed whether the strain difference for this parameter also affects the geometrical complexity of the thalamocortical axon. For this, we determined for each terminal branch (order 0) the number of nodes (axonal bifurcations) that an action potential has to pass before reaching it from its entrance in the cortex. Subsequently, terminal branches were classed according to the number of nodes, and Figure 8 displays their frequency distribution. Statistical analysis identified that this distribution is significantly different between the two strains: terminal branches in BRL axons are reached by passing fewer nodes than in NOR (χ2 test, p < 0.001). This shows that thalamocortical axons in BRL have a longer conduction compartment (LC) but have a simpler branching pattern than those in NOR.
We further analyzed the strain difference in the geometry of the TCAs by testing the correlation between the branch number and the size of their tangential spread. As shown in Figure 9A, in both strains a significant linear correlation between the two parameters is found. However, the correlation is significantly different between the two strains (p = 0.005): compared with NOR, the same variation in branch number produces in BRL a significantly larger increase in tangential spread in layer IV.
The multivariate statistical analyses described above identifies that the BRL mutation affects the parameters characterizing the conduction compartment of the thalamocortical axon as well as the geometry of the axonal arborization. Therefore, some of the results of the univariate analysis are worth mentioning (Table 1).
The total length of the conduction compartment in BRL was 1.6 times longer as in NOR (Table 1, p = 0.02), and its surface of entrance was 31-fold greater (Table 1, p = 0.001). No statistical differences were found for the transmission compartment (Table 1). Significant strain differences were detected for the extent of the tangential surface (S4), which was twice as large in BRL as in NOR (p = 0.002). Also, the mean mediolateral extent of the axons, as determined in the coronal plane, showed a significant difference by being three times larger in BRL than in NOR (BRL, 1.3 mm; NOR, 0.4 mm; p = 0.002).
Also, the correlation function between the tangential surface of the axonal plexus and the axonal length is significantly different between the strains (Fig. 9B,C). This correlation function is, in both strains, significantly different between the plexus in layers I–IV and that in layer V–VI: for NOR, p < 0.01; and for BRL, p < 0.008.
Distribution of bouton densities areas in layers I–IV
As mentioned above, by using the same axonal length, number of boutons, and axonal branches, BRL axons made surfaces twice as large as those of the NOR axons. Therefore, the density of boutons per unit of tangential surface was greater for NOR axons than for BRL axons (p value for mean density analysis). In this final part of the analysis, we studied the distribution of axonal boutons within the plexuses in layer IV and in layer VI. For this, we calculated the density of boutons per 20 × 20 μm area of cortical surface and analyzed the distribution within a series of density ranges. Figure 10 displays the density distribution for three NOR and three BRL axons as well as the relative distribution of the density classes per strain. In NOR, the TCA distributes axonal boutons in a restricted and compact part of layer IV, with a large area of high bouton density range situated in the center of the axonal plexus. In BRL, the total cortical area receiving TCA boutons in layer IV is larger, but the area covered by high bouton density is distributed throughout the plexus, not forming a compact area.
Performing a one-way univariate statistical ANOVA on the absolute values of the density areas inside layer IV confirmed the strain difference and showed that, in BRL, the surface with the lowest bouton density (0,1] was more than twice as large as in NOR (p < 0.001) and formed 33% of the total surface in NOR and 43% in BRL (Fig. 10G–I). Surfaces of intermediate density (1,2], (2,4], and (4,8] were approximately twice as large in BRL than in NOR (p < 0.001; p = 0.006; p = 0.01), and each of them represented 14–20% of the total tangential surface in both NOR and BRL. A smaller part of the total tangential extent of the axons was occupied by surfaces with a high bouton density (8,64] (20% in NOR vs 6% in BRL). This analysis revealed a disproportionate distribution of axonal boutons in the layer IV tangential plexus of BRL axons that resulted in a twofold larger area of low and medium densities as compared with NOR. The high densities (8,64] form a compartment that in absolute value was of similar size in both strains.
Performing the same statistical analysis on the density areas formed by the boutons of the plexus in layers V–VI showed no significant differences between strains (Fig. 10J–L).
Realizing that cortical compartmentalization is not attributable to a single TCA, but several hundreds, we analyzed the distribution of boutons formed by a pair of reconstructed axons that have their cell bodies close to each other in VB. Figure 11A depicts two NOR axons that reached the same cortical barrel and that together formed a maximum density peak of 39 boutons/400 μm2. In BRL, the two axons tended to form their high-density surfaces in the same cortical region (Fig. 11B) and reached a maximal density that was similar with those formed by NOR axons. In BRL, this area of high density is surrounded by a large area of low density. Adding a third BRL axon to this analysis resulted in a even more pronounced delineation of a high-density area (Fig. 11C) that is only slightly larger in its tangential dimension than the barrel-shaped area in NOR. The addition of the bouton distribution of two additional TCAs shows the emergence of a second area of high density (Fig. 11D). From these cumulative histograms, we conclude that in barrelless the high-density domains are at the basis of a topographic representation of the whisker follicles.
Discussion
To elucidate the possible structural basis underlying the barrelless phenotype, we reconstructed TCAs in NOR and BRL mice. Comparing the population of hand-drawn axons with the morphometrical data from the three-dimensional (3D)-reconstructed axons confirms that the conclusions drawn from the quantitative analysis are representative for the total sample of labeled TCAs. These conclusions are the following: (1) The subcortical trajectory of BRL axons is the same as of NOR axons. TCAs in both strains have the same total axonal length, bear the same number of branches and axonal boutons. (2) BRL axons possess branches at the interface between white matter and layer VI, that were not observed in NOR. (3) The number of branches by which TCA axons in BRL reach layer IV is increased in comparison to NOR; also, the areal surface covered by these entering branches is larger in BRL than in NOR. (4) The comparison of the geometry between the two populations of TCA revealed that an action potential in BRL axons has to cross a smaller number of axonal bifurcations before reaching the terminal branches. (5) Within layer IV, BRL axons distribute axonal boutons over a larger cortical area, forming a terminal plexus in which the bouton density is lower when compared with the compact, high bouton dense terminal plexus invariably formed by NOR axons. However, NOR and BRL axons have the same high-density area of axonal boutons, and in both strains, axons with the cell bodies close together in the thalamus made their high-density areas in the same cortical region.
A number of the morphometric parameters (total length, number of branches, and axonal boutons) are indistinguishable between the two strains studied here and, interestingly, have an almost identical value for TCAs terminating in the primary visual cortex of the mouse (Antonini et al., 1999). TCA axons in the primary visual cortex of the tree shrew (Raczkowski and Fitzpatrick, 1990) and macaque (Blasdel and Lund, 1983) bear comparable numbers of boutons as the axons studied here (see supplemental table, available at www.jneurosci.org as supplemental material). This transmodal and interspecies comparison indicates the constant nature of TCA axons in the vertebrate sensory pathways.
Although we analyzed the spatial distribution of axonal boutons in three dimensions, we included in Table 1 the mediolateral extent of our TCA allowing comparison with the literature. NOR axons have a mediolateral extent comparable with TCA in VI of the mouse (Antonini et al., 1999), the barrel cortex of the rat (Jensen and Killackey, 1987a), the primary visual cortex of the tree shrew (Raczkowski and Fitzpatrick, 1990), and the primary visual cortex of macaque (Blasdel and Lund, 1983). Combined with the comparison of the number of boutons, this indicates that the density of the projection of individual TCAs within a tangential area of primary sensory cortex is highly constant between cortical areas and between species. This conclusion is further strengthened by the linear correlation found between axonal length and tangential extent in both strains of mice studied here. Although variability exists between individual TCAs, their contribution to the innervation of the cortex clearly follows common rules. This interstrain/interspecies comparison predicts that variations in the impact of the thalamocortical projections in a given part of the primary sensory cortex are either attributable to a modification in the number of TCAs entering the cortex or the consequence of an alteration in the geometry of TCAs.
Two factors influence the final geometry of TCAs: sensory activity and intracortical neuronal transmission. In NOR, the columnar shape as well as the topological relation between TCAs is established by progressive convergence inside the cortex. In the subcortical trajectory, TCAs with cell bodies close together in the thalamus may lose this spatial relationship. However, from the cortical entry point until their arrival at the boundary between layers IV and V, TCAs redirect to finally converge within the same barrel. This progressive convergence reflects the developmental sequence by which TCAs grow into the cortex and establish whisker-related patterns first in the deep layers before the arrival and pattern formation in layer IV (Agmon et al., 1993, 1995).
Several studies have identified the molecules involved in subcortical pathway finding of TCA (Dufour et al., 2003; Lopez-Bendito et al., 2006) and the signals that TCA follow inside the cortex to terminate in target layers IV and VI (Cang et al., 2005; Shimogori and Grove, 2005). It was further shown that TCA can correct an altered development both during their subcortical trajectory (Higashi et al., 2005) as well as their subsequent growth within cortex (for review, see Rash and Grove, 2006). Although our observations are made in adulthood, it seems that the BRL mutation has neither influenced the subcortical trajectory, nor the termination of TCAs in the appropriate cortical layers.
We showed that BRL mutation disrupts TCA morphology both below and inside layer IV. First, BRL TCAs possess supernumerary WM/VI branches that resemble those described for TCAs during development but that are normally not found in adults (Naegele et al., 1988; Senft and Woolsey, 1991a,b; Ghosh and Shatz, 1992; Catalano et al., 1996; Auladell et al., 2000). Some branches look similar to the “interstitial budding” branches described by O'Leary and Terashima, (1988). Second, the number of axonal branches by which BRL TCAs reach layer IV is increased, and together they form an area of entry in layer IV that is several times larger than that in NOR. Third, in layer IV, BRL axons extend across a cortical surface that is twice that in NOR. These observations indicate that developmental mechanisms as branch pruning/retraction (O'Leary and Terashima, 1988; Skaliora et al., 2000) may be deficient during infragranular and granular stages of TCA development in BRL. This is in agreement with the observations of Rebsam et al. (2002) and Lee et al. (2005) and congruent with the hypothesis that axonal pruning/retraction are mechanisms involved in geometrical shaping of TCA (Senft and Woolsey, 1991a,b; Skaliora et al., 2000; Portera-Cailliau et al., 2005).
These deficiencies can be explained by two, not mutually exclusive, mechanisms. It can be attributable to a dysfunction of BRL TCAs to respond to guidance molecules comparable with the maintenance of ectopic retinal projections in the superior colliculus of barrelless mice (Ravary et al., 2003). In this retinal projection, it was shown that BRL axons maintain aberrant branches because of a modified response to ephrin-5 signaling (Nicol et al. 2006). Alternatively, the extensive BRL arbor can be the consequence of synaptic dysfunction during the first postnatal week when the induction of long-term potentiation was shown to be impaired in the BRL mutant (Lu et al., 2003), a hypothesis that further underlines the role of activity-dependent synaptic plasticity in shaping TCAs (Erzurumlu and Kind, 2001).
In layer IV, a TCA in NOR creates a compact terminal plexus with the size of a barrel, 20% of which is characterized by an area of high density of thalamocortical (TC) boutons. In BRL, the layer IV plexus has an area comparable with up to 2.5 barrels and only 6% of this area is characterized by a high density of TC boutons. The high degree of divergence of TCA termination in layer IV is the basis for the lack of segregation of whisker input to individual layer IV neurons in BRL strain (Welker et al., 1996). The total areal extent of the high-density domains in BRL and NOR are similar, indicating that this aspect of the TCA plexus tends to be resistant to the genetic variation. A similar result has been described by Antonini and Stryker (1993), who demonstrated that decreasing sensory activation in the visual pathway did not alter the high-density area of geniculocortical axons either.
However, the current study reveals that, in BRL, TCAs of neighboring thalamic neurons form a compact area of high-density innervation that has the tangential size of a single barrel. We propose that this high-density area underlies the topologic parcellation of SI in barrelless mice shown previously (Welker et al., 1996). This implies that not all parts of the layer IV plexus of an individual TCA contribute in equal manner to the establishment of a functional map in SI. This notion is in harmony with the Gaussian distribution of TCA boutons in SI of the monkey cortex (Garraghty and Sur, 1990), where the peak of the curve coincides with the corresponding part of the sensory periphery. Enlargement of the functional representation during adult plasticity may therefore be the consequence of a proportional increase in the activation of the part of the TCA termination that has a low bouton density (Garraghty et al., 1989; Garraghty and Sur, 1990; Rausell and Jones, 1995; Jones and Pons, 1998).
In conclusion, we propose that developmental mechanisms that shape thalamocortical axons before their arrival in cortical layer IV are deficient in the BRL mouse, but that thalamic neighborhood relationships are maintained in the cortical projection, in which the spatial distribution of TC boutons underlie the establishment of a functional map of the mystacial whisker follicles.
Footnotes
-
This work was supported by the Swiss National Science Foundation (310000-108246). We thank Laurent Tettoni for support throughout this study; Graham Knott for helpful comments on this manuscript; Nathalie Borrajo-Trapp and Nathalie Mueller for histology; and Gilles Bronchti, Vincent Castagné, Christel Genoud, Charles Quairiaux, Fabien Pichon, Aouatef Abaza, and Jan-Harry Cabungcal for helpful discussions.
- Correspondence should be addressed to Egbert Welker, Département de Biologie Cellulaire et de Morphologie, Faculté de Médecine, Université de Lausanne, Rue du Bugnon 9, CH-1005 Lausanne, Switzerland. Egbert.Welker{at}unil.ch