Selective and constructive mechanisms contribute to neural circuit formation in the barrel cortex of the developing rat

The cellular strategy leading to formation of neuronal circuits in the rodent barrel cortex is still a matter of controversy. Both selective and constructive mechanisms have been proposed. The selective mechanism involves an overproduction of neuronal processes and synapses followed by activity dependent pruning. Conversely, a constructive mechanism would increase the number of axons, dendrites, and synapses during development to match functionality. In order to discern the contributions of these two mechanisms in establishing a neuronal circuit in the somatosensory cortex, morphometric analysis of dendritic and axonal arbor growth was performed. Also, the number of synapses was followed by electron microscopy during the first month of life. We observed that axonal and dendritic arbors retracted distal branches, and elongated proximal branches, resulting in increased arbor complexity. This neuronal remodeling was accompanied by the steady increase in the number of synapses within barrel hollows. Similarly, the content of molecular markers for dendrites, axons and synapses also increased during this period. Finally, cytochrome oxidase activity rose with age in barrels indicating that the arbors became more complex while synapse density and metabolic demands increased. Our results support the simultaneous use of both selective and constructive mechanisms in establishing the barrel cortex circuitry.


INTRODUCTION
During the ontogeny of the nervous system, axonal and dendritic branches, as well as synapses, are thought to be formed in numbers that far exceed the amount required by the mature brain to function properly.The predominant theory for neuronal circuit formation is the selective mechanism, which states that a period of overproduction of neuronal processes and synapses is followed by the elimination of a significant fraction of them.This process of elimination is driven by competition.Axons or dendrites compete for synaptic partners that may provide neurotrophic factors in sufficient amount to consolidate synapses and retain their partners [1].Axonal and dendritic branches that establish synaptic-like contacts with incorrect partners are eliminated because they do not receive enough trophic support.Even though the notion described above constitutes the dominant view about how the developing brain sculpts its circuitry, increasing anatomical evidence supports that the selective mechanism is not a universal rule.Indeed, anatomical evidence gathered during the development of different brain structures indicates that axons may grow towards and within their targets by preferentially adding branches with a high degree of precision until they attain their adult size [2].Similarly, multiple reports also support the notion that the dendritic arbors of some classes of neurons also develop preferentially by adding individual branches until they reach their adult configuration [3].The rodent primary somatosensory cortex (S1) processes tactile information and contains a representation of the body formed by modules called barrels.Each barrel contains an anatomical map that represents collections of mechanosensory receptors distributed throughout the body.
The barrel field is a great tool for understanding how circuits are constructed [4].Barrels primarily contain thalamic afferents from neurons of the ventro postero medial (VPM) nucleus and dendrites of the spiny stellate neurons located in cortical layer IV.The axons arising from neurons in the ventral posterior complex of the thalamus project to cortical neurons located primarily in layer IV of the somatic sensory cortex.Interestingly, both glutamatergic and GABAergic synapses increase their numbers during the first month of life [5][6][7].In addition, there is no evidence of synapse elimination in S1 from the first month of life throughout adulthood [5][6][7].Therefore, it is plausible that other mechanisms contribute to neuronal circuit formation in S1.A possible alternative is the constructive mechanism which states that axons, dendrites, and synapses increase in number during development.According to this model, the number of neural process and synapses should not decrease because there is not elimination of unsuccessful branches.Previous studies, aimed at unraveling the strategy used by S1 elements to establish synaptic connections, have reached conflicting conclusions.While some groups support a selective mechanism [8,9], others suggest a constructive mechanism [10][11][12].In order to get an insight into which of these two possible mechanisms is used for establishing the neuronal circuit from the VPM thalamic nuclei to the rodent primary somatosensory cortex, the dendritic and axonal morphology was analyzed, and the number of synapses was determined at different times during the first month of life.In addition, fractal dimension analysis was also performed to quantify morphological complexity.Finally, the expression of molecular markers for dendrites, axons and synapses was also followed by immunocytochemistry and Western blot.The aim of our work was to simultaneously study the development of the three major components of the neuronal arbor: the axons, dendrites and synapses.Through these analyses we were able to identify different stages of the thalamocortical pathway construction in the barrel cortex.

Animals
Wistar rats were raised at the animal facility of Instituto de Investigaciones Biomédicas, Universidad Nacional Autónoma de México.Litters were adjusted to eight pups each and were weaned at postnatal day 21.Rats were housed with free access to food (Sani-Chips, Harlan Teklad, IN) and water in rooms with controlled temperature and lighting.Animal handling and experimental procedures were conducted according to the Guide for the Care and Use of Laboratory Animals [13].Experimental protocols were also approved according to the guidelines of Lineamientos para el cuidado y uso de animales de laboratorio published by Instituto de Investigaciones Biomédicas.

Dendritic Morphology
At postnatal day 5, 10, 15, 20, 25 and 30, rats were euthanized with pentobarbital (45 mg/kg of body weight; Pet´s Pharma, Nezahualcoyotl, Estado de México, Mexico).Euthanized rats were decapitated and their brains were rapidly dissected, immersed in a Golgi-Cox solution to stain neurons, and kept in light-protected jars for three months at room temperature [14].After this period, 120 μm vibratome (Campden Instruments; WPI, Berlin, Germany) coronal slices were obtained, dehydrated, cleared and placed onto gelatin coated slides, and sealed with a coverslip using Cytoseal 60 (Richard-Allan Scientific, Kalamazoo, MI).Stained cortical layer IV spiny stellate neurons were drawn with the aid of an Optiphot brightfield microscope equipped with a camera lucida (Nikon; Tokio, Japan).Drawings were scanned and digitalized with a scanner model Lexmark X75 at 1200 dpi (Lexmark International de México, S de R.L. de C.V. Mexico City).Digital images were then used to perform a morphometric analysis by counting the total branch number and total branch number per order.Throughout the paper, dendritic branches that grow from the soma are called first order branches, the branches that emerge from the first order ones are called second order branches and so on.The total branch length and total length per branch order were also measured manually by tracing individual branches.Morphometric analyses were carried out using the software Scion Image (Beta 4.0.2Scion-Corp, MD).Values obtained for each parameter were compared statistically among different age groups by using Kruskal-Wallis ANOVA and Tukey post-hoc tests (Sigma Stat 3.0, Systat Software Inc.; San Jose, CA).Differences were considered statistically significant at p < 0.05 values.In addition, to evaluate dendritic arbor complexity, a fractal dimension analysis was performed following the box counting method with the software Harmonic and Fractal Image Analyzer (Version 5.0; HarFa, Brno, Czech Republic, downloaded from http://www.fch.vutbr.cz/lectures/imagesci/). Measuring the fractal dimension is another way in which lengthrelated measurements can distinguish between two objects that have the same linear length but are structurally different [15].

Axonal Morphology
Wistar rats were euthanized and perfused with saline (0.15 M NaCl) followed by a fixative solution containing 4% paraformaldehyde (Sigma-Aldrich, Saint Louis, MO) and 0.5% glutaraldehyde (Sigma-Aldrich) in 0.1 M phosphate buffer, pH = 7.4.After decapitation, heads were placed in amber jars containing the same fixative solution.These jars were kept protected from light at room temperature for a week.Solutions of 0.5 μl of 0.025% DiI (1,1'-dioctadecyl-3,3,3',3' tetramethylindocarbocyanine perchlorate) from Molecular Probes (Eugene, OR) dissolved in ethanol were injected stereotaxically into the ventrobasal complex of the thalamus, while avoiding crossing over into the somatosensory cortex.Heads were then placed back into the same jars.In heads from younger animals (postnatal days 5, 10 and 15), proper DiI labeling was obtained after two months.In heads from older animals, it was necessary to allow dye diffusion for up to nine months.At the end of the dye diffusion-time, brains were collected and were then transected along a plane at 45 degrees to the brain base.Sections of 200 μm from the barrel field were obtained using a vibratome (Campden Instruments, WPI; Berlin, Germany).Sections were placed onto gelatin-coated slides and mounted with glycerol (Sigma-Aldrich).Barrels of the postero-medial subfield were identified and individual axonal arbors within them were drawn with the aid of a camera lucida.Because axons at postnatal day 25 were highly complex, they were reconstructed using a LSM 5 Pascal scanning confocal microscope (Zeiss MicroImaging, Inc.; Overkochen, Germany).The image stacks were later used to generate two-dimensional reconstructions.All drawings were scanned, digitized and analyzed as described above for dendrites.The values obtained for each parameter were compared among different age groups by Kruskal-Wallis ANOVA followed by Dunn's post-hoc test (Sigma Stat 3.0, Systat Software Inc.; San Jose, CA).Differences were considered statistically significant at p < 0.05 values.

Synapse Density
Wistar rats from postnatal days 5, 15 and 30 were euthanized and perfused with saline solution followed by Karnovsky's fixative (5% glutaraldehyde EM grade, 4% formaldehyde in 0.08 M sodium phosphate buffer).Brains were obtained and cortices dissected, flattened and placed for at least 24 hours at 4˚C in fixative solution.Vibratome tangential sections (200 -250 μm) were cut and used to identify C3 and D3 barrels of the postero-medial barrel subfield by trans-illumination under a dissecting microscope (M420 Stereomicroscope.Leica Microsystems; Wetzlar, Germany).Once identified, barrels were dissected using a 2 mm diameter dissecting needle as previously reported [16].Tissue slabs were placed in Karnovsky's fixative for 24 h, washed three times with sodium carcodylate buffer 0.1 M, for 1 h at 4˚C and then in 1% osmium tetroxide (Sigma-Aldrich) for 1 h at 4˚C.After three rinses with water, tissue slabs were dehydrated with ethanol and then immersed in propylene oxide for 40 minutes.Finally tissue slabs were immersed in propylene oxide-Epon for 1 h, and then included in Epon (Electron Microscopy Sciences; Pennsylvania, USA).Tissue slabs were then cut using an ultramicrotome (Leica Ultracut EM UC6), and three ultrathin sections 65 μm apart from each other were mounted in formvar-coated grids (Electron Microscopy Sciences).This allowed sampling of each barrel at its upper, middle and lower levels.To avoid observer's bias, Epon blocks containing the samples were coded, before being sectioned by a technician who was not involved in performing final observations and analyses.The average number of synapses per 100 square microns was then estimated.

Antibody Characterization
Source and properties of antibodies used in this paper are summarized in Table 1.The antibody against microtubule associated protein-2 (MAP2; Chemicon International Inc., Temecula CA) was used to analyze the developmental pattern of dendritic arbors of neurons at cortical layer IV.This antiserum recognized two bands of approximately 300 kDa by Western blot of rat cortex lysate.Labeling of cell bodies and neuropil by the anti-MAP2 antibody was consistent with the pattern observed in rat and chicken brains [17,18].The rabbit anti-neurofilament antibody (Chemicon International) was used to analyze the development of thalamocortical afferents.The staining pattern was similar to patterns observed in previous studies of rat brain [19].This antibody recognized a band of 150 kDa by Western blot analysis, and also labeled axonal terminals that innervate Merkel cells in rat skin Mouse; monoclonal (data not shown).To analyze the synaptogenic process during S1 development, we performed immunostaining against pre-and postsynaptic markers.The synaptophysin antibody (DAKO A/S, Glostrup, Denmark) reacted with a protein of 38 kDa band by Western blot analysis.The punctate labeling pattern obtained with the anti-synaptophysin antibody was similar to that obtained in neural cultures and in central nervous system [20,21].The antibodies against the postsynaptic markers PSD 95 (Chemicon International) and Shank2 (Zymed Laboratories Inc., San Francisco, CA) also gave a punctuate pattern as described previously [22].The antibody against PSD 95 recognized two bands around 100 kDa and a single band of 71 kDa in rat brain lysates.

Immunohistochemistry
Rats were euthanized and perfused with saline followed by 4% paraformaldehyde in phosphate-buffered saline, pH 7.4.Brains were dissected and postfixed in the same fixative for two hours at room temperature and cryoprotected with 20% sucrose (Sigma-Aldrich) in 0.1 M phosphate-buffered saline.Brains were then frozen in 2-methylbutane chilled on dry ice.Serial coronal slices (30 μm thick) were cut in a cryostat model CM 1510-3 (Leica).Free-floating slices were collected in blocking solution (3% bovine serum albumin, 0.1% triton X-100 and 0.025% sodium azide in phosphate buffered saline).Slices were incubated in blocking solution at 4˚C for 12 h and then incubated overnight with primary antibodies diluted in blocking solution, at 4˚C.The primary antibodies used were rabbit anti-rat/human Shank-2 at 1/200 dilution; rabbit antihuman synaptophysin at 1/5000 dilution; rabbit anti-bovine neurofilament 150 kDa or anti medium-sized neurofilament protein at 1/1000 dilution; and mouse anti-bovine microtubule-associated protein 2 a&b (MAP2a & b) (Chemicon International) at 1/1000 dilution.Next, slices were washed three times, 15 minutes each, with 0.1% triton X-100 in phosphate-buffered saline and then incubated with biotinylated goat antirabbit IgG or with goat anti-mouse IgG antibodies (1/500 dilution) at room temperature for 1.5 hours.After three more washes, slices were incubated for 1 hour with the avidin-peroxidase complex (VECTASTAIN ® Elite ABC; Vector Laboratories, Inc. Burlingame, CA) at room temperature.Then, after three more washes slices were developed with diamobenzidine following the manufacturer's instructions (VECTASTAIN ® Elite ABC kit).Finally, slices were washed, mounted onto glass slides coated with gelatin, air dried, and sealed with a cover slip using Cytoseal 60 (Richard-Allan Scientific, Kalamazoo, MI).Images of the immunostained barrels were captured and digitized at 40× (for MAP2, and neurofilaments) or at 100× (for synaptophysin, and Shank-2) using a brightfield microscope (Axioskope, Zeiss Mi-croImaging, Carl Zeiss Microscopy, LLC, Thornwood, NY).

Western Blot
Cortices of euthanized rats were quickly removed, flattened and primary somatosensory cortex samples were obtained as previously described [23].Samples were homogenized on ice by three 15 seconds pulses of sonication (Capitol Scientific, Inc., Austin, TX) in 200 μl lyses buffer (100 mM NaCl, 1% Triton X-100, 1 mM phenylmethylsulfonyl fluoride, 10 mM Tris-HCl, pH = 8.0) supplemented with protease inhibitors (Complete ® ; Roche Indianapolis, IN).Then, samples were centrifuged at 14000 rpm for 20 minutes at 4˚C in a refrigerated microcentrifuge model 5415 (Eppendorf, Hauppauge, NY).Supernatant was collected and protein concentration determined by the Bradford method (Bio-Rad, Richmond, CA).Proteins (20 μg protein/lane) were resolved by SDS-PAGE for MAP2 and neurofilaments as described previously or by Tris-Tricine electrophoresis for synaptophysin and PSD-95 as described previously [24,25].Resolved proteins were then transferred to nitrocellulose membranes in a semidry blot chamber (Bio-Rad) for 30 minutes at room temperature.Membranes were then blocked with TTBS (0.15 M NaCl, 5% calcium-free milk, 20 mM Tris, 0.05% Tween 20, pH = 7.5) for 12 h at 4˚C.Membranes were then incubated overnight with anti-MAP2 (1/1000), anti-neurofilaments (1/5000), anti-PSD95 (1/500), or anti-synaptophysin (1/5000) antibodies at 4˚C.After washing with TTBS three times, 15 minutes each, membranes were incubated with the corresponding biotinylated secondary antibodies (1/500 dilution) for two hours at room temperature.Membranes were incubated with avidin-peroxidase complex for one hour at room temperature, washed with TTBS three times, 15 minutes each, and developed by chemiluminescence (ECL, Amersham-Pharmacia-Biotech, Buckingham-shire, UK).Images on films were digitized (Fluor-S Bio-Rad) and densitometry values were obtained using Quantity One software (Bio-Rad).Briefly, a rectangle was drawn over the background area of the film and then positive bands were outlined manually.The program calculated the absorbance values by subtracting the background from each positive band.For each neural marker, the absorbance values were expressed in the graphs relative to the absorbance at PD5 which was considered as one.

Cytochrome Oxidase Staining
Primary somatosensory cortex barrels were visualized by cytochrome oxidase staining as previously described [26].Briefly, brain free-floating slices were incubated in cytochrome oxidase solution (0.05% diamobenzidine 0.03% cytochrome c, and 0.02% catalase in 0.1 M phosphate buffer) for four hours at 37˚C.The reaction was halted by washing three times with PBS.Slices were mounted onto glass slides coated with gelatin, air dried, and sealed with a coverslip using Cytoseal 60 (Richard-Allan Scientific, Kalamazoo, MI).Images of the barrels were captured and digitized at 40× or 100× using a brightfield microscope (Axioskope, Zeiss MicroImaging).Captured images were used to evaluate the density of the staining (Scion Corp, Beta 4.0.2).Barrels were outlined manually and absorbance values were determined.A similar procedure was used for the corpus callosum.A relative estimation of the label intensity was calculated as follows: (a − b) × 100/b.Where a is the light transmittance of each barrel, and b is the light transmittance of the corpus callosum (background).

Dendrites of Spiny Stellate Neurons Increased Arbor Complexity While Remodeling within the First Month of Life
Dendrites of spiny stellate neurons within layer IV of the primary somatosensory barrel cortex are essential for receiving information from the thalamus.The number of branches in the dendritic arbor reflects the complexity of the thalamocortical pathway.Thus, in order to follow the growth and complexity changes of dendritic arbors of spiny stellate neurons at layer IV of primary somatosensory cortex (S1), neurons were identified by Golgi-Cox staining and fifty cells, per age were drawn with the help of a camera lucida.Neurons within the barrel subfield at postnatal day (PD) 5 presented very few branches in the dendritic arbors.Branches were short and single, or poorly ramified (Figure 1A), while neurons at PD15 had many more branches and more complex dendritic arbors (Figure 1A).Neurons at PD25 had longer dendritic arbors with a similar complexity as in PD15.Notably at PD25 the branches were polarized to one side of the cell, probably because the branches were oriented towards the center of the barrel to contact terminal axonal arbors of the thalamus.It is generally accepted that hollows emerge via passive displacement of cortical cells by dense thalamocortical terminal clusters in barrel centers.Additionally, it seems that the emergence of sides and septa are also related to the progressive asymmetry of dendrites from the cells of the barrel side toward the barrel hollow [27].In order to quantify dendritic arbors at different ages, morphometric analyses of these neurons was performed.Box chart whisker plots were used to show the distribution of total dendritic arbor branch number (Figure 1B), total dendritic arbor length (Figure 1C), and dendritic arbor fractal dimension (Figure 1D).Total dendritic arbor branch number increased as the neurons matured, reaching a maximum at PD15 (Figure 1B).By PD25 the total dendritic arbor branch number seemed to decrease slightly but this difference was not significant (Figure 1B).Similarly, the total length of dendritic arbors increased as neurons matured, reaching a maximum size at PD25 (Figure 1C).Also, fractal dimension of these dendritic arbors increased as the neurons matured, and reached a maximum at PD15 (Figure 1D).By PD25 the fractal dimension remained constant (Figure 1D).Since total dendritic arbor branch number and total dendritic arbor length did not fully explain complexity; we further examined the level of complexity by determining dendritic arbor branch number and length per order.These parameters reflect the complexity of a circuit and can show the construction of the whole arbor.Dendritic arbors of neurons at PD5 primarily had dendrites of first order, which branched out reaching the fourth order (Figure 1E).Dendritic arbors of neurons at PD15 showed dendritic branches that reached the ninth order (Figure 1E).The largest number of dendritic arbor branches were those of the first to the forth order (Figure 1E).Dendritic arbors of neurons at PD25 showed a similar dendritic arbor branching behavior as PD15 neurons, except that ninth order branches were not detected (Figure 1E).These data suggest that neurons expand dendrite branching within the first 15 days of postnatal life, and indicate that dendrites have more branching close to the soma.In PD5 neurons, first order dendritic branches were the longest and dendritic arbor branch length decreased as the branch order augmented (Figure 1F).In PD15 neurons, second and third order dendrites were the longest (Figure 1F).From the forth order the branch length decreased gradually as the order increased (Figure 1F).In PD25 neurons, second and third order branch dendrites were still the longest ones and also were longer than those present on PD15 neurons (Figure 1F).Similarly, as in PD15 neurons, the dendritic arbor branch length decreased gradually as the order increased (Fig- ure 1F).These results suggest that as the neurons matured, they initially increased the number and complexity of their dendrites, followed by branch loss across dendrite orders, while they simultaneously increased the length of the remaining dendrites.

Thalamocortical Axons Increased Arbor Complexity While Remodeling within the First
Another parameter that provides information on how neurons make connections is the number of branches in the terminal axonal arbors (axonal arbors).Therefore, the growth and complexity of axonal arbor branches from thalamocortical neurons was also evaluated.Axons were stained with DiI labeling and drawn with the help of a camera lucida and confocal imaging.Axonal arbors at PD5 presented few branches that ramify poorly (Figure 2A), while axonal arbors at PD15 had more branches and were more complex (Figure 2A).Axonal arbors at PD25 seemed to have fewer branches than PD15 arbors, but these branches were much longer (Figure 2A).
In order to quantify the axons, morphometric analyses of these thalamocortical axonal arbors were also performed.Box chart whisker plots were used to show the distribution of total axonal arbor branch number (Figure 2B), total axonal arbor length (Figure 2C) and axonal arbor fractal dimension (Figure 2D).Total axonal arbor branch number increased as the neurons matured, reaching a maximum at PD15 (Figure 2B).By PD25 the total axonal arbor branch number decreased (Figure 2B).In contrast, the total axonal arbor length was almost the same until PD15, but increased significantly as they reached PD25 (Figure 2C).Also, the fractal dimension of the axonal arbors increased from PD10 to PD15, and reached a maximum at PD15.There were no further significant changes in complexity by PD25 (Figure 2D).To gain additional insight into the axonal complexity, the axonal arbor branch number and length per order for these axons was determined.These parameters can show the dynamics of branching for these axonal arbors.Similarly to dendrites, first order axonal arbor branches emerge from the main axon that arrives to the barrel cortex layer IV.Axonal arbors at PD5 showed branches that reached the ninth order (Figure 2E).Axonal arbors at PD15 showed branches that reached the 17th order.The largest number of branches was present from the fifth to the ninth order (Figure 2E).Similarly, axonal arbors at PD25 showed branches up to the 14th order (Figure 2E).This is in agreement with the idea that some branches from the previous age were lost.The largest number of axonal branches at PD25 was found in the sixth to the eight orders (Figure 2E).Axonal arbor branch length per order was similar between PD5 and PD15.At PD25, the largest axonal arbor branches were measured from third to ninth order even though some branches were lost (Figure 2F).These results suggest that as thalamocortical axons matured, they increased the length and complexity of their branches.

Microtubule-Associated Proteins 2 and Neurofilaments Increased in the Primary Somatosensory Cortex during the First Month of Life
Microtubule-associated proteins 2a and 2b (MAP2) are mainly confined to cell bodies and dendrites.To visualize neural processes within layer IV of the primary somatosensory barrel cortex neurons, a stain against MAP2 proteins was used.In PD5 neurons, MAP2 proteins are expressed in cell bodies and branches arising from them (Figure 3A).In PD15 neurons, the amount of neuronal processes and cell bodies expressing MAP2 proteins increased considerably (Figure 3B).The amount of neuronal branches expressing MAP2 was similar in PD25 neurons (Figure 3C).Apical dendrites from layer V pyramidal cells showed MAP2 immunostaining at all ages (Figures 3A-C).These results suggest that dendritic arbors become more stable with age.In agreement with the immunostaining, Western blot analyses showed that the content of MAP2 in the developing the primary somatosensory cortex increased from PD5 to PD30 (Figure 3G).The ratio of expression of MAP2 relative to PD5 augmented with time reaching a two-fold increase by PD30 (Figure 3I).Neurofilaments are a type of intermediate filament that functions as major structural element of the axonal cytoskeleton.Neurofilaments antibodies were used to visualize thalamic terminal axonal arbors within layer IV of the primary somatosensory barrel cortex.Immunoreactivity for neurofilaments increased inside barrels from PD5 to PD15 (Figures 3D and E).Following this age, neurofilament staining persisted in the axons and was also observed in the cell bodies of cortical neurons (Fig- ure 3F).Western blot analyses showed an incremental increase in the content of neurofilaments in the primary somatosensory cortex from PD5 to PD15.The neurofilament content remained similar from PD15 to PD30 (Figure 3H).The ratio of neurofilament expression relative to PD5 augmented with time reaching a two-fold increase by PD30 (Figure 3J).

Synaptic Markers and Synaptic Density Increased within Layer IV during the First Month of Life
As dendritic and axonal arbors were evaluated, the functional unit of the interaction between both structures, the synapse, was also analyzed.The presence of two postsynaptic markers (Shank-2 and postsynaptic density protein 95) and a presynaptic marker (synaptophysin) were analyzed in layer IV of the primary somatosensory barrel cortex.Immunohistochemistry for Shank-2 at PD5 revealed a punctuate label mainly inside the cell bodies with some localization within the neuropil (Figure 4A).At PD15 the immunostaining remained in cell bodies but also was found around cells and in the neuropil (Figure 4B).By PD25 punctuate labeling of Shank-2 in the cell bodies was intensified and staining could also be seen all over the barrel in the neuropil (Figure 4C).Western blot analyses showed PSD-95 was barely detectable at PD5 but it increased progressively from PD15 to PD25 with a significant incremental increase from PD25 to PD30 (Figure 4G; p < 0.05).The ratio of PSD-95 expression relative to PD5 augmented with time reaching an almost three-fold increase by PD30 (Figure 4I).Label for synaptophysin was found surrounding the cell bodies of layer IV neurons.This labeling increased in PD15 and became more prominent in PD25 ( Figures 4D-F).Likewise, Western blot analyses for synaptophysin showed a gradual increase from PD5 to PD30 (Figure 4H).The ratio of synaptophysin expression, relative to PD5 was augmented over time reaching an almost four-fold increase by PD30 (Figure 4J).This observation supports the idea that synapses are first established on cell bodies and later in the surrounding processes located in the neuropil.To estimate synaptic density within the primary somatosensory barrel cortex, electron microscopy sections were collected at PD5, PD15 and PD30 (Figures 5A-C).Synaptic density doubled from PD5 to PD15, while it increased 20% from PD15 to PD30 (Figure 5D).These results suggest that as the synaptic machinery matures, synaptogenesis increased with age.

Cytochrome Oxidase Activity Increased within Layer IV during the First Month of Life
To evaluate the degree of synaptic activity in layer IV of primary somatosensory barrel cortex, mitochondrial activity was measured using cytochrome oxidase histochemistry.The histochemistry for this enzyme has been widely used to localize barrels and to determine synaptic activity [26].In coronal sections, somatosensory cortex was readily recognized by cytochrome oxidase stained barrels at PD5, PD15 and PD25.Barrel staining intensity increased over time (Figures 6A, C and E, respectively).Using a higher optical magnification, sites of mitochondrial activity were observed throughout the barrel cortex including cell bodies and neuropil (Figures 6B, D and F).In summary, mitochondrial activity, estimated by densitometry of cytochrome oxidase staining increased within layer IV during the first month of life (Figure 6G).

DISCUSSION
The mechanism(s) underlying rat somatosensory thalamocortical pathway development is still a matter of controversy.Previous observations about thalamocortical somatosensory circuitry formation were based upon studies that addressed single aspects of the development of the circuit: axonal growth, dendritic growth, or synaptogenesis.The isolated study of neuronal components could lead to the impression that a brain region uses only one mechanism to construct its final circuitry.Here, we have quantified aspects of all the components of the circuit, axons, dendrites, and synapses to better understand which mechanism is underlying the circuit formation.
Our data suggest that both selective and constructive mechanisms operate concomitantly during circuit forma-   tion in the rat barrel cortex.We propose that somatosensory thalamocortical arbors formation can be explained by a two-phase model (Figure 7).
During the first phase (growth phase), the dendritic and axonal branches are profusely added, as we observed during the first 15 days after birth.At the same time, synapses start forming in the soma.The second phase (remodeling phase) is characterized by events of elimination and elongation of both dendritic and axonal arbors branches.Synapses continue in the somas, but appear also in the neuropil.In agreement with this notion, the amount of MAP2 in the neuronal cytoskeleton of the barrel cortex increased rapidly after PD15, a phenomena that has been associated with dendritic remodeling during brain development [28].Moreover, the neurofilament content also increased over time which indicates the formation and growth of axonal terminals.Simultaneously, dendritic and axonal branches are lost throughout the whole arbor configuration, especially at higher branch orders.While this shaping process is occurring in the arbors, the amount of synapses also increased with age.The two-phase model proposed in this work is in concordance with similar events described during the development of dendritic and axonal arbors in sensory, motor and higher-order brain regions of invertebrates and vertebrates [29][30][31].This two-phase model is also contained in the concurrent formation-elimination theorical model proposed for axon development in which the selective growth of some branches and the elimination of others occur simultaneously [32].By contrast, Purkinge cell dendrites grow in a phase manner, but in this case the cells first retract branches and, then add many branches before remodeling [33].Synapse formation can also be explained under the framework of the two-phase model.Indeed they are added continuously during development: first in the soma and then also in the neuropil.The axonal presynaptic marker synaptophysin increases in an almost linear manner, while the postsynaptic marker PSD-95 takes longer to increase.Moreover, our data indicates that at the beginning of postnatal development, the pre-and postsynaptic markers in the barrel cortex were primarily associated with neuronal cell bodies and then, as development progresses, the distribution of these markers shifts to the neuropil.Thus, it seems that the synaptogenesis in barrel cortex starts with a low number of synaptic contacts associated with the neuronal cell bodies.Thereafter, an intense period of synaptic maturation at the neuropil starts after PD15 until the end of the first month of life.This pattern of synaptogenesis is similar with the centrifugal process of synapse formation in the cerebellum, where a change in synapse localization also occurs [34,35].The fact that synaptic density increased in barrel neuropil during the first month of postnatal life also is in agreement with the increase of glutamatergic and GABAergic synapse density profiles that has been described previously in barrel cortex [6].As we observed in this work, the maximum number of synapses in the cortex is reached by PD30.Interestingly, the synapse density in barrels remains relatively constant through adulthood without any evidence of synaptic elimination [7].These observations suggest that synapse formation in barrels uses a constructive mechanism, as in the case for the olfactory bulb, rather than a mechanism in which synapses are overproduced early in development and then eliminated [36][37][38][39].Previous evidence suggests that pre-and postsynaptic sites mature in an asynchronous fashion [40][41][42].Axons acquire their ability to release neurotransmitters well before they reach their targets and induce the organization of postsynaptic sites upon contacting them, a process that triggers synapse formation [43][44][45][46][47].This asynchrony likely explains why, in our western blot analyses, the content of synaptophysin increases steadily from PD5 to PD30.This suggests to us that the maturation of the pre-synapses occurs continuously.In contrast, the content of PSD-95 remains relatively low from PD5 to PD20.After PD20, PSD-95 content increases rapidly and abruptly.So, in this case the post-synapses take longer to mature and then stabilize faster as development continues.Thus, it is possible that the profiles of synapses observed before PD15 with the aid of the electron microscopy correspond to imamture synapses.The fact that the content of PSD-95 and the immunostaining for Shank-2 increase significantly after PD15 suggests that synaptic consolidation starts after this date and continues during the phase of remodeling.This might explain the difficulty to induce LTP in the somatosensory cortex after PD15 and the disappearance of silent synapses in barrels before this postnatal day [48][49][50].In support to this conclusion, it is known that PSD-95 is required for AMPA receptor incorporation and hence for synaptic consolidation [51][52][53].Also, activation of AMPA receptors inhibits actin-based spine motility, thus promoting contact stabilization [54].It is also note-worthy that synaptogenesis promotes branch elongation [51,55].Therefore, we think that this centrifugal pattern of synapse formation partly explains why, during the remodeling phase, proximal branches of dendritic arbors elongate whereas distal branches are eliminated.In addition, synaptogenesis also promotes sprouting [56].The dominance of one of these processes during the events of remodeling may determine if in some neuronal processes the length is greater than the ramification magnitude or vice versa.Once the core structure of the circuit is developed, there is evidence that the remodeling process continues throughout life.Indeed, in the adult brain circuit changes in synapse formation and elimination depend on structural plasticity which is thought to underlie aspects of long-term memory formation [57].Neuronal activity is the basis for many mechanisms underlying pattern formation, plasticity and information processing [58].One of the most established ideas in neurosciences is that neuronal activity promotes dendritogenesis, axonogenesis and synaptogenesis.Noteworthy, the stimulus-evoked cortical activity detected in S1 has been reported at the end of the first week of postnatal life using 2 deoxiglucose uptake [59][60][61][62].Additionally, differences in activity patterns before and after the onset of whisking behavior may be caused by neuronal activity induced by whisking, during active touch [61].In this regard, it has been argued that the activity of the enzyme cytochrome oxidase reflects the levels of evoked neuronal activity [63].Thus, our data showing the increase of cytochrome oxidase staining during the first two weeks likely reflected the metabolic demands associated with neuronal growth.After PD14, cytochrome oxidase staining might also reflect use dependent neural activity caused by whisking.In addition, the increase in cytochrome oxidase staining could reflect changes in the patterns and/or levels of neuronal activity associated with passive stimulation as seen in lactating rats and in immature somatosensory cortex [64,65].Indeed, metabolic indices utilized commonly to monitor indirectly the activity levels could be sensitive to changes in activity patterns [66].In the present report, we analyzed the developmental pattern of thalamic terminal axons, spiny stellate dendritic arbors, and synapses in the barrel cortex, during the first month of life.Our results suggest that this circuit does not use a single strategy to build up its connectivity, but rather follows a two phase model, concurrently elongating proximal and eliminating distal den-dritic and axonal branches as synaptogenesis proceeds.Thus, we argue both the selective and constructive mechanisms are simultaneously used during somatosensory circuit formation in specific time frames during development.

Figure 1 .
Figure 1.Dendrites of spiny stellate neurons increase complexity while remodeling within the first month of life.(A) Representative camera lucida drawings from spiny stellate neurons within layer IV of the primary somatosensory barrel cortex.The drawings represent neurons at postnatal day (PD) 5, 15 and 25.Scale bar 40 μm.(B-D) Box-and-whisker plots of dendritic arbors, indicating the mean (squares), 25th percentile (bottom line), median (middle line), 75th percentile (top line), 5th and 95th percentile (whiskers), and minimum and maximum measurements (x).Morphometric analysis of spiny stellate dendritic arbors was performed by counting total dendritic arbor branch number (B), total dendritic arbor length (C), and dendritic arbor fractal dimension (D) at PD5, PD10, PD15 and PD25.Numbers in parenthesis represent the number of dendritic arbors measured for each condition.5 rats were used per age.Dendritic arbor branch number per order (E) and dendritic arbor branch length per order (F) were also measured from the same dendritic arbors at PD5 (open squares), PD15 (close circles), and PD25 (open triangles).

Figure 2 .
Figure 2. Thalamocortical axons increase complexity while remodeling within the first month of life.(A) Representative camera lucida and confocal drawings from thalamocortical axons within layer IV of the primary somatosensory barrel cortex.The drawings represent axon terminal arbors at postnatal day (PD) 5, 15 and 25.Scale bar 100 μm.(B-D) Boxand-whisker plots of axonal arbors, indicating the mean (squares), 25th percentile (bottom line), median (middle line), 75th percentile (top line), 5th and 95th percentile (whiskers), and minimum and maximum measurements (x).Morphometric analysis of thalamocortical axonal arbors was performed by counting total axonal arbor branch number (B), total axonal branch length (C) and axonal arbor fractal dimension (D) at PD5, PD10, PD15 and PD25.Numbers in parenthesis represent the number of axonal arbors measured for each condition.5 rats were used per age.Axonal arbor branch number per order (E) and axonal arbor branch length per order (F) were also measured from the same axonal arbors at PD5 (open squares), PD15 (close circles), and PD25 (open triangles).

Figure 3 .
Figure 3. Microtubule-associated proteins 2 and neurofilaments increased in the primary somatosensory cortex during the first month of life.Brain sections of layer IV of the primary somatosensory barrel cortex (S1) of postnatal day (PD) 5 (A, D), PD15 (B, E) and PD25 (C, F), were immunostained for MAP2 (A-C) and neurofilaments (D-F).Arrowheads show dendritic and axonal stained branches.White arrows show stained cell bodies.Black arrows illustrate apical dendrites from pyramidal layer V neurons.Scale bar 20 μm. 5 rats were used per age.Representative Western blots for MAP2 (G) and for neurofilaments (H) at PD5, PD10, PD15, PD20, PD25 and PD30 of protein lysates from S1 samples.Pools of 10 S1 areas were used per age (n = 5).Band intensity was estimated by a densitometer and values were normalized to the value of PD5 who was taken as 1.Average ± SE of five Western Blots are shown for MAP2 (I), and for neurofilaments (J).* Represent statistical differences between PD5 vs. all ages.

Figure 4 .
Figure 4. Synaptic markers increase within layer IV during the first month of life.Brain sections of layer IV of the primary somatosensory barrel cortex (S1) of postnatal day (PD) 5 (A, D), PD15 (B, E) and PD25 (C, F), were immunostained for Shank-2 (A-C) and synaptophysin (D-F).Label for these pre and postsynaptic proteins are founded in cell bodies (arrow heads) and in the neuropil (white arrows).Cell bodies are dotted lined.Scale bar 10 μm.Representative Western blots for PSD-95 (G) and for synaptophysin (H) at PD5, PD10, PD15, PD20, PD25 and PD30 of protein lysates from S1 samples.Pools of 10 S1 areas were used per age (n = 5).Band intensity was estimated by a densitometer and values were normalized to the value of PD5 who was taken as 1.Average ± SE of five Western Blots are shown for PSD-95 (I), and for synaptophysin (J).* Represent statistical differences between PD5 vs. all ages.& Represent statistical differences between PD25 vs. PD30.

Figure 5 .
Figure 5. Synaptic density increase within layer IV during the first month of life.Representative electronic micrographs of layer IV primary somatosensory barrel cortex (S1) show synapses (arrows) at postnatal day (PD) 5 (A), PD15 (B) and PD30 (C).Scale bar 1μm.D) Plot showing changes in synaptic density at different ages.C3 and D3 barrels were used to count synapses, sampling of each barrel at its upper, middle and lower levels (n = 5 rats per age).* Represent statistical differences between PD5 vs. all ages; p < 0.05.

Figure 6 .
Figure 6.Cytochrome Oxidase activity increased within layer IV during the first month of life.Representative microphotographs show histochemistry against cytochrome oxidase of coronal sections of primary somatosensory cortex at postnatal day (PD) 5 (A and B), PD15 (C and D) and PD25 (E and F).At a lower magnification (panelsA, E, and E) barrels are clearly seen (encircled by dotted lines).At a higher magnification (panels B, D, and F) cell bodies are clearly seen (encircled by dotted lines).Mitochondrial activity sites are also denoted by arrow heads in the periphery of cell bodies and by white arrows in the neuropil (n = 5 rats per age).Scale bar 100 μm for lower magnification panels and 20 μm for higher magnification panels.The graph (G) shows the quantification for the cytochrome oxidase staining.* Represent statistical differences between PD5 vs. all ages; p < 0.05.

Figure 7 .
Figure 7.A two-phase model is followed by the somatosensory thalamocortical pathway to establish neural circuits in the barrel cortex.(A) Initial state of axonal and dendritic arbors development and synaptogenesis at postnatal day 5. (B) First phase of the model that consists of the addition of dendritic and axonal arbors branches at postnatal day 15.(C) Second phase or remodeling phase of elimination-elongation of dendritic and axonal arbors branches postnatal day 30.

Table 1 .
Information of primary antibodies used.