JN Information on EB 2010
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Neurophysiol 100: 327-335, 2008. First published May 7, 2008; doi:10.1152/jn.90392.2008
0022-3077/08 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Figures
Right arrow All Versions of this Article:
100/1/327    most recent
90392.2008v2
90392.2008v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via Web of Science (2)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Chomiak, T.
Right arrow Articles by Hu, B.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Chomiak, T.
Right arrow Articles by Hu, B.

Functional Architecture and Spike Timing Properties of Corticofugal Projections From Rat Ventral Temporal Cortex

T. Chomiak, S. Peters and B. Hu

Department of Clinical Neurosciences and, Hotchkiss Brain Institute, Faculty of Medicine, University of Calgary, Calgary, Alberta, Canada

Submitted 24 March 2008; accepted in final form 29 April 2008


 ABSTRACT
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
Sensory association and parahippocampal cortex in the ventral temporal lobe plays an important role in sensory object recognition and control of top-down attention. Although layer V neurons located in high-order cortical structures project to multiple cortical and subcortical regions, the architecture and functional organization of this large axonal network are poorly understood. Using a large in vitro slice preparation, we examined the functional organization and spike timing properties of the descending layer V axonal network. We found that most, if not all, layer V neurons in this region can form multiple axonal pathways that project to many brain structures, both proximal and remote. The conduction velocities of different axonal pathways are highly diverse and can vary up to more than threefold. Nevertheless for those axonal projections on the ipsilateral side, the speeds of axonal conduction appear to be tuned to their length. As such, spike delivery becomes nearly isochronic along these pathways regardless of projection distance. In contrast, axons projecting to the contralateral hemisphere are significantly slower and do not participate in this lateralized isochronicity. These structural and functional features of layer V network from the ventral temporal lobe may play an important role in top-down control of sensory cue processing and attention.


 INTRODUCTION
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
The temporal association and parahippocampal cortex (TeA) is a poly-sensory structure that plays a critical role in high-order sensory recognition and control of top down attention. According to current neurophysiological model of attention control (Desimone and Duncan 1995Go; Fuster 1997Go; Gilbert and Sigman 2007Go), sensory experience and awareness can be represented as network-driven persistent neuronal firing in the form of working memory or template of attention (Desimone and Duncan 1995Go; Fuster 1997Go; Gilbert and Sigman 2007Go). When this internally generated template activity is propagated to downstream cortical and subcortical structures, it leads to an enhancement of the salience of the attended environmental stimulus, thereby creating a selection bias in sensorimotor processing (Desimone and Duncan 1995Go; Fuster 1997Go; Gilbert and Sigman 2007Go; Komura et al. 2001Go; Miyashita 2004Go; Naya et al. 2001Go).

Although the cortical projection system is a critical element for the TeA network, the functional organization of specific feedback pathways remain poorly defined. In particular, there is insufficient evidence on whether and how individual projection cells can form a large axonal network that requires multiple target pathways transcending the boundaries of different sensorimotor modalities. It is also unknown how fast spike conduction and synchronized activation can be maintained in such a network. Large-scale axonal network with diffuse projections not only raises uncertainty in terms of the fidelity of signal conduction but also the issue of spike timing. Recent studies have shown that spike timing can encode neural information in the cortex on a scale of a few milliseconds (deCharms and Merzenich 1996Go; Steinmetz et al. 2000Go) or less (Shmiel et al. 2005Go). The timing of presynaptic spike volleys and their relative temporal order to a postsynaptic spike can also exert strong influence on neuronal plasticity and development (Hausser et al. 2001Go; Markram et al. 1997Go). Hence, large temporal jitters or delays in axonal conduction may disrupt these timing-based mechanisms of information processing.

Anatomical tracing studies in different species have thus far shown that deep layer neurons in motor, sensory, and temporal association cortex project to multiple cortical and subcortical regions (Arnault and Roger 1990Go; Bourassa and Deschenes 1995Go; Cruikshank et al. 2002Go; Deschenes et al. 1994Go, 1998Go; Shi and Cassell 1997Go; Vaudano et al. 1991Go; Winer et al. 2001Go). These studies, however, do not necessarily provide data on the prevalence of the uncovered projection pathways. That is, are they derived from highly branched axons from single neurons or primarily formed by many individual cells with fewer axonal branches? Moreover it is not feasible to demonstrate the relative spike timing properties of the cerebral cortical networks based on pathway tracing. To address these issues, we developed a large horizontal slice preparation where multiple electrodes can be placed in different brain regions for antidromic activation of axonal pathways. This slice preparation was previously used to demonstrate that layer V neurons in the rat TeA can maintain a myelinated and an nonmyelinated collateral pathway in the posterior thalamus (Chomiak and Hu 2007Go). In this study, the architecture, branching pattern and conduction velocity from individual layer V axons were evaluated. The results show that most, if not all, layer V neurons in this brain region can form multiple axonal pathways that project to many brain structures, both proximal and remote, and a striking isochronicity of spike conduction was found in the layer V axonal network that was lateralized to ipsilateral projections.

Because layer V neurons convey top-down cortical modulatory signals to many other sensorimotor structures (Bereshpolova et al. 2006Go; Guillery and Sherman 2002Go; Jones 2001Go; Ma and Suga 2001Go; Winkowski and Knudsen 2006Go), it was hypothesized that single layer V neurons may utilize multiple axonal collaterals and an isochronic conduction mechanism to synchronize its spike delivery to multiple targets located at different distances. Preliminary results have been reported previously in abstract form (Chomiak and Hu 2006Go).


 METHODS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
In vivo electrophysiology

All experimental protocols were approved by the University of Calgary Animal Care Committee. Male Long Evans or Sprague Dawley rats (~300 g) were initially anesthetized with sodium pentobarbital (Somnotol) and later maintained at stable anesthetic level via intravenous drip or intraperitoneal infusion of 85% ketamine and 15%/xylazine. Lidocaine hydrochloride (1%, Sigma) was injected subcutaneously in the scalp and lidocaine gel (Xylocaine) applied to the skin prior to surgical incisions. After a craniotomy, bipolar electrodes for local field potential (LFP) recording were lowered to the following locations relative to bregma according to the rat atlas (Paxinos and Watson 1986Go): TeA, anterior/posterior –7.0, medial/lateral 6.0, dorsal/ventral 4.5; caudate putamen (CPU), –1.4, 4.0, 4.0; medial geniculate dorsal division/lateral posterior thalamic nuclei (PT), –5.6, 3.4, 4.0; basolateral amygdala: –0.28, 5.0, 8.0; Par, –4.5, 4.0, 2.0. Single or trains of monophasic pulse (150 µA, 100–250 µs) stimulation were applied via a bipolar stimulating electrode to the white matter immediately beneath the TeA. LFPs were amplified, digitized at 4–5 kHz (National Instruments NIDAQ 6300E) and captured on a laptop running Igor software (Wavemetrics). Onset latency of the LFPs for each site was calculated by averaging 15 consecutive traces and taking the peak of the second derivative of the smoothed average. At the end of experiments, recording and stimulation sites were marked via electrical current lesion. Animals were given Somnotol (100 mg/kg ip) and perfused intracardially with saline or 4% paraformaldehyde. Brains were extracted and stored in PFA for histological inspection of recording and stimulation locations.

Slice preparation

Briefly, male Sprague-Dawley or Long-Evans rats (125–200 g) were anesthetized with halothane and decapitated and the brain was removed. For analysis of long-range connections, a horizontal slice preparation was used (Chomiak and Hu 2007Go), similar to that reported by Cruikshank et al. (2002)Go. The horizontal slices (~425–450 µm) were cut on the vibrotome (Leica) and incubated in artificial cerebrospinal fluid (ACSF) for 1–2 h before usage. Slices were submerged in and perfused with oxygenated (95% O2-5% CO2), warmed (32–34°C) ACSF. The ACSF had a final pH of {approx}7.4, an osmolality of 300 ± 3 mosM/kg and contained (in mM) 122 NaCl, 3 KCl, 1.3 MgCl2, 25.9 NaHCO3, 2.5 CaCl2, and 11 glucose. The patch electrode solution (pH 7.2) contained (in mM) 120 K gluconate, 10 KCl, 10 Na-HEPES, 0.2 Na-GTP, and 4 Mg-ATP. The pipette solution was adjusted to a final osmolality of 283 ± 3 mosM. The pipette solution was buffered to pH 7.2 ± 0.1. All electrophysiological experiments were performed at 30–33°C.

In vitro electrophysiology

Recordings were restricted to layer V neurons in the area 2 (Te2) and the adjacent parahippocampal area 36 (dorsal) of the ventral association temporal cortex. In rat, these two areas are collectively known as TeA, corresponding to the poly-sensory temporal cortex present in monkey and human (Burwell and Amaral 1998Go; Palomero-Gallagher and Zilles 2004Go; Shi and Cassell 1997Go; Vaudano et al. 1991Go). Neurons included in the database were either regular spiking (RS) or intrinsic bursting (IB) type, characteristic of layer V pyramidal cells (Hefti and Smith 2000Go). The liquid junction potential was subtracted on-line and the bridge balance was constantly monitored and compensated for. Patch electrodes were made from borosilicate glass capillary tubing (1.5 mm OD and 1.1 ID with filament, WPI) on a Flaming/Brown puller. The DC resistance of the patch electrodes was 6–8 M{Omega}. Signal acquisition and analysis was accomplished with the Multiclamp 700A and DIGIDATA 1322A 16-bit data-acquisition system digitized at 20 kHz and filtered at 10 kHz and with the Clampex 9 data-acquisition and analysis programs by Axon Instruments (Forster City, CA). Concentric or side-by-side (separated by 500 µm) bipolar stimulating electrodes (World Precision Instruments) were placed in different target areas where many terminals from the TeA projections are found. Single stimulus pulses (typically <110 µA and 200 µs) were delivered alternatively to each site through a Digitimer DS2 or A-M Systems 2100 isolated stimulator controlled by DS8000 Digital pulse generator (World Precision Instruments). The individual stimulation sites in our slice were separated by >1 mm. According to previous studies (Bagshaw and Evans 1976Go; McIntyre et al. 2004Go; Nowak and Bullier 1996Go), the estimated current spreading edge in our preparation is generally <500 µm. Compared with fast-conducting fibers, slow-conducting fibers typically required higher stimulation intensities (Chomiak and Hu 2007Go). Antidromic latencies or conduction times were determined by measuring the time interval between the end of the stimulation pulse and the onset (foot) of the antidromic evoked action potential (Lipski 1981Go) (see Fig. 3). Similar latency data were obtained using concentric and side-by-side bipolar electrodes, consistent with previous reports (Bagshaw and Evans 1976Go).


Figure 3
View larger version (17K):
[in this window]
[in a new window]

 
FIG. 3. Population latency distribution of TeA axonal pathways. Left: an example action potential evoked by antidromic stimulation of axonal terminal field regions. The axonal conduction time is defined here as the time interval (tt) between the stimulation artifact and the onset of the initial segment spike (IS-SD) or the "foot" of the AP. Right: population histograms of axonal spike latency distributions for each of the four TeA efferent pathways. Scale bars: horizontal 0.5 ms, vertical 20 mV.

 
Histology

For fiber tract tracing, rats were anesthetized with isoflurane (4% induction, 2–2.5% maintenance), and a Hamilton syringe containing 80% biotin dextran amine (BDA, Molecular Probes) was lowered stereotaxically to the temporal lobe associational cortex. After 5 min of settling time, injection proceeded manually at a rate of 200 nl every 3 min to a total volume of 1 µl BDA. Six minutes was allowed prior to syringe removal to ensure proper diffusion of the tracer. The animals were allowed to recover for 24–48 h before transcardially perfused and fixed at 4°C for over 24 h.

For pathway length measurement, horizontal slices (100 µm) were cut in the same way as those used for electrophysiology. The sections were treated with a 0.3% H2O2 solution in methanol for 20 min to block the remaining endogenous peroxidase activity. They were incubated in avidin–biotin complex solution (1:100; ABC Elite kit, Vector Laboratories) overnight at room temperature. After three washes for 10 min each in 1x phosphate-buffered solution (PBS), the sections were reacted with a 3,3V-diaminobenzidine solution, rinsed with distilled water, and covered with a coverslip. Same staining reactions were also performed as controls on uninjected side of the brain sections (n = 4) and in rats received only vehicle injections (n = 3) or with missed target injections (n = 3). There was no endogenous BDA background staining or false positive deposit of BDA in these specimens. Clearly identifiable BDA-positive axons from both gray matter and white matter where measured using the imaging software provided by Northern Eclipse (Empix Imaging).

Pathway tracking and length measurement

Our experimental design required that the full length of long axonal efferents from TeA can be measured along multiple trajectories and terminal sites in the same slices that were used for electrophysiological latency measurement. We found, under this circumstance, slicing the brains with pre-BDA injection gave the best results. In contrast, we were unable to uncover the axonal branches of a layer V cell to any significant length via intracellular dye injection (unpublished observations). The BDA-labeled pathways were first manually traced and drawn to first generate a silhouette of BDA-labeled fibers (see Fig. 2 for example). This was done by employing two-dimensional color enhancing of brown BDA products followed by template subtracting (Empix Imaging). The silhouette was then overlaid to the photograph of live slice preparation to help localize the pathway trajectory and stimulation electrode locations (Fig. 1 B). Our tracing data revealed that the axonal trajectories and target sites similar to that described in the literature (Cruikshank et al. 2002Go; Meltzer and Ryugo 2006Go; Vaudano et al. 1991Go). To localize the gray matter structures in the slice, we used the atlas by Paxinos and Waton and by Kruger et al. (Kruger et al. 1995Go; Paxinos and Watson 1986Go), the latter of which provides clear demarcations between gray and white matters based on myelin staining. To more accurately measure fiber length with curvature, the pathway images were first overlaid onto the brain slice image captured via a CCD camera (Fig. 1B). Using curvature measuring tool (Empix Imaging), the absolute pixel length, including its angular proportions, for each pathway branch was calculated between the stimulation and recording sites and converted into millimeters. As previously noticed, long contralateral branches (CLBs) unlike other pathways exhibited poor tracer retention (Parent and Parent 2006Go) and cannot be kept with full length in the in vitro horizontal slice. However, CLBs have been shown to almost exclusively course through the external capsule (EC) and corpus callosum (Hackett et al. 1999Go; Parent and Parent 2006Go; Vaudano et al. 1991Go). We therefore took the shortest route connecting the EC, callosum and then the contralateral EC (stimulation site) as the CLB pathway length (Fig. 1B; ···).


Figure 2
View larger version (43K):
[in this window]
[in a new window]

 
FIG. 2. Pathway tracing and geometrical measurements of TeA efferent axons. Middle (A): montage of TeA descending pathways prelabeled by BDA. Light microscope images were taken (x10) from a horizontal slice that was used for electrophysiology. Injection of BDA (brownish deposit at the TeA site) was made 48 h prior to tissue processing as described in METHODS. BDA containing fibers appear in the external (EC) and internal capsules as well as in the subcortical white matters. Prominent fibers with terminal like morphology are also seen in PT, CPU, ECIC, and tectum. Some amplified (x20) images are shown as insets for EC and ECIC sites. The red arrow indicates the location where descending fibers entering the PT and tectum are split from those to CPU. Left (B): TeA descending pathways were manually traced and laid over to the computer-generated silhouette as shown in gray (Northern Eclipse). Right (C): images of BDA-labeled pyramidal cells at the TeA injection site. Bottom graph: the presence of dense plexuses of labeled thin fibers or swellings (also see enlarged images in A) in the gray matter, which is taken as the indication of terminal fields. It is noteworthy that subcortical projections to nonthalamic regions arise primarily from layer V neurons. Abbreviations are given in text and Fig. 1.

 

Figure 1
View larger version (25K):
[in this window]
[in a new window]

 
FIG. 1. Characterization of functional temporal association and parahippocampal cortex (TeA) efferent pathways in vivo and in vitro. A: in vivo recordings of local field potentials (LFPs) in cortical and subcortical targets evoked by ipsilateral TeA stimulation (shown as individual traces). LFPs show short-latency (<4 ms) responses in the caudate putamen (CPU), parietal cortex (Par), amygdala (Amyg), and posterior thalamic nuclei (PT). The data were compiled from 6 rats for each recording site. Inset scale bars: vertical, 70 µV; horizontal, 5 ms. B: brain slice used for in vitro pathway tracing and electrophysiology. Axonal pathway trajectories (orange) are overlaid onto the horizontal preparation. Blue dots indicate stimulation sites. Inset (bottom): a typical layer V neuron with an approaching patch electrode. Inset (left): higher-magnification images of biotin dextran amine (BDA)-positive axons traveling within the external capsule located beneath the TeA (indicated by the red arrow). ECIC, external cortex of the inferior colliculus; HF, hippocampal formation; Ic, internal capsule; EC, external capsule; CC, corpus callosum; and CLB, contralateral branch.

 
Statistical tests

Data are expressed as means ± SD if they pass the D'Agostino and Pearson normality test. Otherwise they are expressed as geometric mean ±95% confidence interval (CI). P values were either determined from the test statistic derived from a two-sampled t-test (parametric), Kruskal-Wallis (nonparametric ANOVA) test followed by Dunn's multiple-comparison test (nonparametric).


 RESULTS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
In vivo tract tracing and orthrodromic stimulation

Our first goal was to obtain a slice preparation where the length of some of the axonal pathways from TeA can be anatomically traced and their time delays measured, preferably in the same cell and with minimal confounding factors. To ensure the tracer injection was targeted at the TeA area, we first conducted experiments in vivo to see if orthrodromic stimulation in this brain region can evoke synaptic responses in multiple subcortical sites. We stimulated the TeA and recorded orthrodromically evoked local field potentials (LFPs). In anesthetized rats, single electrical stimulus in TeA area elicited LFPs in the CPU, parietal cortex (Par), limbic region, PT (Fig. 1A), indicating these brain sites indeed contain axonal terminals from TeA. The postsynaptic nature of the recorded LFPs was indicated by their long duration (typically lasted from 10 to 20 ms), much longer than action potentials and their sensitivity to stimulation frequency (potentiation or depression; data not shown). Furthermore, we found that the average LFP latency of all ipsilateral sites were rapid and synchronized (mean: 3.1 ± 0.26 ms and comparable among almost all ipsilateral sites, coefficient of variance = 8%, P > 0.05; n = 6). However, in spite of corticotectal axons forming axonal plexuses of terminals in the external cortex of the inferior colliculus (ECIC) (Arnault and Roger 1990Go), the evoked LFPs in this region were small and variable. This is likely due to terminal field dispersion and/or its susceptibility to anesthesia. It should be noted that the waveforms of LFPs recorded in vivo may be subject to the influences of anesthetics. Moreover, their latencies include synaptic and membrane ion channel delays at postsynaptic sites, in addition to axonal delays (Hu 1995Go; Marsalek et al. 1997Go). Therefore although LFPs may reveal the sites of axonal endings, the absolute latency values of LFPs cannot be equalized to presynaptic axonal conduction time.

In a separate set of experiments, we injected a small amount of BDA into the same site as we applied electrical stimulation. Brains from these rats were cut into horizontal slices (115 µm) with a 30° angle (similar to that used for electrophysiology) to keep all the branches of TeA pathways as complete as possible (see METHODS). In these slices, BDA-labeled fibers could be followed from TeA to CPU, PT, tectum, and ECIC (Figs. 1 B and 2). The CLBs, although they could be activated in the slices by stimulating the contralateral external capsule, were not labeled by the BDA. This likely reflects the poor ability for tracer retention of callosal fibers (Hackett et al. 1999Go; Parent and Parent 2006Go).

Antidromic activation of layer V neurons in vitro

Whole cell recordings were made from layer V pyramidal neurons maintained in in vitro horizontal preparation (Fig. 1B). The majority of layer V neurons can be classified as bursting type although some neurons, especially in younger rat group, regular spiking, and initially did not seem to fire action potentials in response to current injections or antidromic stimulation (unpublished observations). Single electrical stimulus applied at threshold intensities to each of the four stimulation sites all evoked a single fast action potential (AP). The waveform and the latency of a typical AP are shown in Fig. 3 (left). As previously noted (Chomiak and Hu 2007Go), the stimulus threshold for evoking an antidromic AP from different stimulation sites can differ significantly, which is in part related to the local branch pattern and level of tissue myelination (e.g., between posterior thalamus and brain stem). During repetitive antidromic stimuli, antidromic responses exhibited small (<0.1 ms) temporal jitters in antidromic response. The antidromic nature of the evoked axonal responses were further confirmed by the presence of collision tests by pairing an antidromic AP with either an orthodromic AP or an antidromic AP initiated from another collateral and by the absence of synaptic potentials during antidromic stimulation at different membrane potentials (Chomiak and Hu 2007Go; Lipski 1981Go).

Population latency and conduction velocity

We constructed the spike latency histograms for each axonal pathway obtained from 61 neurons (Fig. 3, right). At the population level, the mean onset latency of antidromic spikes of all the ipsilateral pathways clustered around 0.57 ± 0.37 ms (Fig. 3B), demonstrating an identical synchronized timing pattern relationship seen in vivo (Fig. 1A). In contrast, the average latency from all responding fibers in the contralateral projection was significantly longer than all ipsilateral pathways (1.2 ± 0.12 ms; n = 34; P < 0.001; see Table 1). Although as a population, latency variation among ipsilateral pathways appears to be somewhat variable, intrabranch comparison in single cells showed an opposite trend (see following text).


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

 
TABLE 1. Summary of individual efferent pathways properties

 
Next we calculated the mean conduction velocities of individual axonal pathways (Table 1). The axonal length data were obtained by measuring the distance of each fiber projection from TeA to the stimulation sites (Fig. 1B). These pathways are prelabeled with an anterograde tracer deposit and cut in same way as the slices used for electrophysiology (see METHODS). As shown in Table 1, the conduction velocities of individual pathways are evidently heterogeneous.

Spike conduction among different branches from the same cell

If axonal geometrical properties reflect a self-adaptation process of individual fibers, how does this regulation manifested among individual axonal branches in the same cell? In other words, are there relative timing patterns between individual collaterals? To address this question, we evaluated the relative timing pattern differences between axonal branches in individual neurons. As shown in Fig. 4, two distinct timing patterns emerged. First, the three ipsilateral pathways exhibit a very small intrabranch latency shifts among themselves (0.24 ± 0.21 ms; n = 20). Second, and in contrast, the latency difference between the three ipsilateral branches and the CLB of the same cell is disproportionately large (0.78 ± 0.31 ms; n = 10), with the CLB delays being shifted by approximately equal to threefold (P < 0.001; Fig. 4). Note, however, that this value underestimates the true intrabranch latency shift because the entire length of the contralateral branch cannot be maintained in the preparation. Nevertheless, the timing pattern of the ipsi- and contralateral pathways showed a systematic asymmetry. The conduction delays of ipsilateral pathways tend to be minimized and tuned to isochrony, whereas the contralateral fibers tend to be significantly slower and lag the ipsilateral spike conduction, suggesting a co-existence of different axonal timing strategies at a level of single cells.


Figure 4
View larger version (11K):
[in this window]
[in a new window]

 
FIG. 4. Intra-pathway delays among individual branches from single TeA neurons. Average intra-pathway latency difference among all ipsilateral branches and the contralateral branch (i.e., the difference between the contralateral site with at least one (and ≤3) of the ipsilateral sites for an individual cell). The phase or latency shifts between contra- and ipsilateral branches are significantly different. See text for details. Scale bars: horizontal, 2 ms; vertical, 25 mV. **P < 0.001.

 
Velocity-length relationship

We constructed the velocity-length relationship for layer V axons. This method was used previously to evaluate the degree of synchrony in spike conduction time among parallel axonal projections in olivocerebellar pathways (Sugihara et al. 1993Go). We plotted the observed conduction velocities (Vobs) of individual branches against their pathway lengths (Dt) as shown in Table 1. It has been demonstrated that for a network to achieve isochronous spike conduction (ICT), a linear relationship between velocity (Vobs) and pathway length Dt must exist (Sugihara et al. 1993Go). This relationship is expressed as Vobs = k·Dt; where k (ms–1) is a rate constant and equivalent to the reciprocal of ICT (i.e., k = 1/ICT). In Fig. 5, the experimental velocity and distance data are plotted for all the fiber pathways. We found a quasi-linear relationship in velocity-distance plots for experimental data with a best-fit k value of 2.2 ± 0.2 ms–1. Strikingly, the contralateral projection did not fit into this linear plot (Fig. 5; "contralateral").


Figure 5
View larger version (10K):
[in this window]
[in a new window]

 
FIG. 5. Velocity-length relationship of layer V efferent axons. The velocity (Vobs)-distance plot is based on the data in Table 1. The isochronous spike conduction ICT) is represented by the inverse slope (1/k) of the line. The velocity distance demonstrates a linear relationship (R2 = 0.94) and can be described as Vobs = k(Dt), where the slope factor k equals 1/ICT (ms–1). The best fit value for k is 2.2 ± 0.2.

 

 DISCUSSION
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
In the present study, we examined timing properties of the axonal network formed by cortical layer V neurons. We first found that a single layer V neuron can issue multiple axonal collaterals of different length. This reticular-like axonal network, however, exhibits robust isochronic spike conduction. Moreover, the ICT appears to be asymmetrical and is strongly biased toward ipsilateral projections.

Methodological considerations

STIMULATION SITES.  The placement of our stimulation electrodes was targeted at brain regions that have been previously shown to host rich axonal terminal labelings when anterograde transport tracers were deposited in the TeA area (Arnault and Roger 1990Go; Shi and Cassell 1997Go; Vaudano et al. 1991Go; Winer et al. 2001Go). This labeling pattern has also been reconfirmed in our own tracing experiments. Although we cannot be certain of the exact site of spike initiation during stimulation (which is relevant to all stimulation studies), it is important to keep in mind that the errors associated with this are exceedingly small relative to the total path-length (i.e., a few hundred micrometers vs. several millimeters) and fall within the 95% confidence interval range of the experimental data (Table 1).

STIMULATION CURRENT SPREAD.  Our antidromic collision tests showed that individual axonal pathways can be selectively activated in slices. This result cannot be due to current spread. The individual stimulation sites in our slices were separated by distances >1 mm. We estimated that with a threshold current intensity <110 µA, the edge of spreading current in our preparation does not seem to exceed 500 µm. This is consistent with that reported in previous studies (Bagshaw and Evans 1976Go; McIntyre et al. 2004Go; Nowak and Bullier 1996Go). Furthermore, we have also shown that slow conducting fibers require higher stimulation intensities than fast conducting fibers, also confirming that different types of axonal collaterals can be activated independently in our preparation (Chomiak and Hu 2007Go).

ANTIDROMIC LATENCY MEASUREMENTS.  The latency and conduction time were determined based on the time interval between the end of the stimulation pulse and the onset (foot) of the antidromically evoked action potential (Lipski 1981Go) (see Fig. 3). This is because antidromic spike initiation in the axon follows the termination or peak of a stimulus pulse, which in turn reflects the brief time needed to charge the membrane of an axon (Davies 1968Go; McIntyre and Grill 1999Go). Indeed even in our preparation, we found that shortening the pulse duration of an activating stimulus (e.g., from 150 to 100 µs for a given intensity) led to an all-or-none activation failure, consistent with the fact that axonal excitation occurs near the terminating edge of the stimulus pulse (data not shown). Using the end of the pulse for latency determination can also help to improve the comparability of population date where different stimulus parameters are used for activating different axons (Ranck 1975Go). Furthermore, because the pulse durations we used here are similar to the apparent utilization time of axons, the latencies measured at the end of stimulus pulse is, in theory, more closely match to the actual conduction time (Jankowska and Roberts 1972Go).

Descending reticular network

TeA in rat comprises the areas corresponding to the poly-sensory temporal cortex present in monkey and human, including part of the dorsal aspect of parahippocampal area 36 (Burwell and Amaral 1998Go; Palomero-Gallagher and Zilles 2004Go). Previous axonal tracing studies in rats show that the TeA projects to multiple cortical and subcortical regions (Arnault and Roger 1990Go; Shi and Cassell 1997Go; Vaudano et al. 1991Go; Winer et al. 2001Go). These include the ipsilateral conjoint regions of the dorsal medial geniculate body and lateral posterior nucleus (Burwell and Amaral 1998Go; Shi and Cassell 1997Go), caudal portion of the CPU (Vaudano et al. 1991Go), the ECIC (Arnault and Roger 1990Go), the superior colliculus (Arnault and Roger 1990Go; Vaudano et al. 1991Go), brain stem nuclei (Arnault and Roger 1990Go; Vaudano et al. 1991Go), the amygdala and peri-/entorhinal cortices (Burwell and Amaral 1998Go; Shi and Cassell 1997Go), parietal and occipital cortices (Burwell and Amaral 1998Go; Shi and Cassell 1997Go), prefrontal and orbital cortices (Shi and Cassell 1997Go; Vaudano et al. 1991Go). TeA also sends projections to the contralateral TeA region (Vaudano et al. 1991Go). However, it remains uncertain that efferent projections to these high-order brain structures are mainly derived from the same neuron or formed predominately by many disparate cell groups without significant axonal branching.

In this study, using electrophysiological methods, we provide strong evidence that the long distance axonal projections from layer V neurons are mainly derived from collaterals arising from the same cell. Although previous studies have already alluded this structural organization in sensory and motor cortical neurons (Bourassa and Deschenes 1995Go; Deschenes et al. 1998Go; O'Leary and Koester 1993Go), to our knowledge, this is the first unambiguous demonstration that most, if not all, layer V neurons form divergent axonal pathways to multiple subcortical regions. The cortical neurons we recorded from also frequently showed antidromic spikes to ipsilateral intra-cortical stimulation (unpublished observations). However, due to shorter axonal length, the latency of these intra-cortical branches cannot be determined as precisely as other pathways. Therefore from a structural point of view, the axonal architecture of layer V axons may form a "reticular" network through which a cortically generated signal can be distributed effectively via presumably a small group of projection neurons, transcending the boundaries of many sensorimotor modalities. Interestingly, this type of descending network has already been proposed for the association and adjacent parahippocampal cortex as an important network structure for facilitating sensory recognition and memory recall (Damasio 1989Go). The highly branched innervation pattern may allow time-locked co-activation of the same geographically separate sites in early sensory pathways that were engaged during the initial experience (Damasio 1989Go; Fuster 1997Go; Gilbert and Sigman 2007Go).

Conduction velocity and isochronous spike conduction

Although layer V projections are considered fast-conducting fibers, the antidromic latency data obtained in vitro slice can be approximately equal to fivefold faster than that reported in vivo (Beloozerova et al. 2003Go; Swadlow and Weyand 1981Go). Several technical factors may account for this difference, notably the absence of anesthetics, interference of spontaneous activity (Swadlow 1998Go), difference in intracellular versus extracellular measurements of spike waveform latency (Lipski 1981Go), and the requirement of membrane depolarization for observing antidromic spikes (Chomiak and Hu 2007Go). Nevertheless, it is important to note that fast-conducting cortical antidromic responses, with comparable latencies, have also been observed in vivo (Sirota et al. 2005Go; Swadlow 1991Go).

In a number of central neural pathways such as olivocerebellar (Lang and Rosenbluth 2003Go; Sugihara et al. 1993Go), thalamocortical (Salami et al. 2003Go), and amygdalocortical projections (Pelletier and Pare 2002Go), the axonal delays appear rather uniform and isochronous. In our study, the k value (2.2 ms–1) appears to be slightly different from that found in olivocerebellar and thalamocortical fiber systems (Salami et al. 2003Go; Sugihara et al. 1993Go). This difference may be due to the fact that we used pure axonal conduction velocity to calculate ICT without implicating synaptic delays as was done in other studies. This may yield a shorter but more accurate conduction time for axon collaterals and, as a result, allow us to separate out composite conduction delays (i.e., synaptic and axonal). In support of this conjecture, we found that the average in vivo latencies of othrodromically evoked LFPs in ipsilateral pathways, although isochronous, were longer than antidromic spike latencies measured in vitro. Finally, the diameters of the axons in layer V network appear larger (up to ~3 µm; Supplementary Fig. S11 ) than those of the thalamocortical system (see Salami et al. 2003Go).

Possible mechanism of ICT

Despite a number of axonal networks that are known to exhibit ICT or strong tendencies toward synchronous spike conduction, it is not intuitively apparent how this is achieved. The conduction time of an axon is largely determined by its geometrical properties, e.g., the diameter, length, and degree of myelination (Lang and Rosenbluth 2003Go; Waxman and Bennett 1972Go; Waxman and Swadlow 1977Go). To account for our experimental data, we derived a simple model (see supplementary material) based on Rushton's quantitative theory of axon diameter and speed of conduction for myelinated and nonmyelinated fibers (Rushton 1951Go). With the use of experimentally determined structural data (Supplementary Figs. S1 and S2), we found that the appearance of ICT in this layer V network can be readily accounted for by regional myelination and axonal geometry adaptation among ipsilateral branches (Supplementary Fig. S3). Interestingly, in our model, the conduction time of CLBs from layer V neurons does not fit the isochronous spike timing pattern of the ipsilateral pathways despite having a conduction velocity comparable to the fastest ipsilateral branches (Fig. 5 and Table 1). We interpret this result as that the longer path-length of CLBs (>10 mm) incurs an extra time delay (Supplementary Fig. S3). These results are consistent with the suggestion that both regional myelination and axonal geometry play an important role in the occurrence of ICT (Salami et al. 2003Go).

Functional implications

Corticofugal projection network plays an important role in top-down attention modulation and facilitation of sensory cue responses in the basal ganglia, tectum, and posterior thalamus (Komura et al. 2001Go; Ma and Suga 2001Go; Mooney et al. 2004Go; Winkowski and Knudsen 2006Go). By having the long axonal fibers consolidated to a relatively few projection neurons, layer V network may significantly enhance its functionality. For example, the spike timing control can be executed faster and more precisely in this way, thereby facilitating orientating response to behaviorally significant events that often requires speedy sensorimotor coordination across different modalities. Fast spike conduction in the ipsilateral layer V network may also help dispatch cortical attention signals to multiple brain regions simultaneously. This may allow an internal frame of when, what, and where, to be rapidly reconstituted across different sensorimotor modalities that are physically separate in space but become closely connected in time.


 GRANTS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
This study was supported by grants from Canadian Health Research Institute and its Regenerative Medicine and Nanomedicine Initiative.


 FOOTNOTES
 
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1 The online version of this article contains supplemental data. Back

Address for reprint requests and other correspondence: T. Chomiak, Dept. of Neuroscience, HSC, Room 2042, Faculty of Medicine, University of Calgary, 3330 Hospital Dr. N.W., Calgary, Alberta T2N 4N1, Canada (E-mail: tgchomia{at}ucalgary.ca)


 REFERENCES
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
Arnault P, Roger M. Ventral temporal cortex in the rat: connections of secondary auditory areas Te2 and Te3. J Comp Neurol 302: 110–123, 1990.[CrossRef][Web of Science][Medline]

Bagshaw EV, Evans MH. Measurement of current spread from microelectrodes when stimulating within the nervous system. Exp Brain Res 25: 391–400, 1976.[Web of Science][Medline]

Beloozerova IN, Sirota MG, Swadlow HA. Activity of different classes of neurons of the motor cortex during locomotion. J Neurosci 23: 1087–1097, 2003.[Abstract/Free Full Text]

Bereshpolova Y, Stoelzel CR, Gusev AG, Bezdudnaya T, Swadlow HA. The impact of a corticotectal impulse on the awake superior colliculus. J Neurosci 26: 2250–2259, 2006.[Abstract/Free Full Text]

Bourassa J, Deschenes M. Corticothalamic projections from the primary visual cortex in rats: a single fiber study using biocytin as an anterograde tracer. Neuroscience 66: 253–263, 1995.[CrossRef][Web of Science][Medline]

Burwell RD, Amaral DG. Cortical afferents of the perirhinal, postrhinal, and entorhinal cortices of the rat. J Comp Neurol 398: 179–205, 1998.[CrossRef][Web of Science][Medline]

Chomiak T, Hu B. Isochronic spike conduction of layer V cortical neurons. Soc Neurosc Abstr 634.18/G2, 2006.

Chomiak T, Hu B. Axonal and somatic filtering of antidromically evoked cortical excitation by simulated deep brain stimulation in rat brain. J Physiol 579: 403–412, 2007.[Abstract/Free Full Text]

Cruikshank SJ, Rose HJ, Metherate R. Auditory thalamocortical synaptic transmission in vitro. J Neurophysiol 87: 361–384, 2002.[Abstract/Free Full Text]

Damasio AR. Time-locked multiregional retroactivation: a systems-level proposal for the neural substrates of recall and recognition. Cognition 33: 25–62, 1989.[CrossRef][Web of Science][Medline]

Davies PW. Classical electrophysiology. In: Medical Physiology (12th ed.), edited by Mountcastle VB. St. Louis, MO: Mosby, 1968, p. 1073–1076.

deCharms RC, Merzenich MM. Primary cortical representation of sounds by the coordination of action-potential timing. Nature 381: 610–613, 1996.[CrossRef][Medline]

Deschenes M, Bourassa J, Pinault D. Corticothalamic projections from layer V cells in rat are collaterals of long-range corticofugal axons. Brain Res 664: 215–219, 1994.[CrossRef][Web of Science][Medline]

Deschenes M, Veinante P, Zhang ZW. The organization of corticothalamic projections: reciprocity versus parity. Brain Res Brain Res Rev 28: 286–308, 1998.[CrossRef][Medline]

Desimone R, Duncan J. Neural mechanisms of selective visual attention. Annu Rev Neurosci 18: 193–222, 1995.[CrossRef][Web of Science][Medline]

Fuster JM. Network memory. Trends Neurosci 20: 451–459, 1997.[CrossRef][Web of Science][Medline]

Gilbert CD, Sigman M. Brain states: top-down influences in sensory processing. Neuron 54: 677–696, 2007.[CrossRef][Web of Science][Medline]

Guillery RW, Sherman SM. Thalamic relay functions and their role in corticocortical communication: generalizations from the visual system. Neuron 33: 163–175, 2002.[CrossRef][Web of Science][Medline]

Hackett TA, Stepniewska I, Kaas JH. Callosal connections of the parabelt auditory cortex in macaque monkeys. Eur J Neurosci 11: 856–866, 1999.[CrossRef][Web of Science][Medline]

Hausser M, Major G, Stuart GJ. Differential shunting of EPSPs by action potentials. Science 291: 138–141, 2001.[Abstract/Free Full Text]

Hefti BJ, Smith PH. Anatomy, physiology, and synaptic responses of rat layer V auditory cortical cells and effects of intracellular GABA(A) blockade. J Neurophysiol 83: 2626–2638, 2000.[Abstract/Free Full Text]

Hu B. Cellular basis of temporal synaptic signalling: an in vitro electrophysiological study in rat auditory thalamus. J Physiol 483: 167–182, 1995.[Abstract/Free Full Text]

Jankowska E, Roberts WJ. An electrophysiological demonstration of the axonal projections of single spinal interneurons in the cat. J Physiol 222: 597–622, 1972.[Abstract/Free Full Text]

Jones EG. The thalamic matrix and thalamocortical synchrony. Trends Neurosci 24: 595–601, 2001.[CrossRef][Web of Science][Medline]

Komura Y, Tamura R, Uwano T, Nishijo H, Kaga K, Ono T. Retrospective and prospective coding for predicted reward in the sensory thalamus. Nature 412: 546–549, 2001.[CrossRef][Medline]

Kruger L, Saporta S, Swanson LW. Photographic Atlas of the Rat Brain: The Cell and Fiber Architecture Illustrated in Three Planes With Stereotaxic Coordinates. Cambridge, UK: Cambridge, 1995.

Lang EJ, Rosenbluth J. Role of myelination in the development of a uniform olivocerebellar conduction time. J Neurophysiol 89: 2259–2270, 2003.[Abstract/Free Full Text]

Lipski J. Antidromic activation of neurones as an analytic tool in the study of the central nervous system. J Neurosci Methods 4: 1–32, 1981.[CrossRef][Web of Science][Medline]

Ma X, Suga N. Corticofugal modulation of duration-tuned neurons in the midbrain auditory nucleus in bats. Proc Natl Acad Sci USA 98: 14060–14065, 2001.[Abstract/Free Full Text]

Markram H, Lubke J, Frotscher M, Sakmann B. Regulation of synaptic efficacy by coincidence of postsynaptic APs and EPSPs. Science 275: 213–215, 1997.[Abstract/Free Full Text]

Marsalek P, Koch C, Maunsell J. On the relationship between synaptic input and spike output jitter in individual neurons. Proc Natl Acad Sci USA 94: 735–740, 1997.[Abstract/Free Full Text]

McIntyre CC, Grill WM. Excitation of central nervous system neurons by nonuniform electric fields. Biophys J 76: 878–888, 1999.[Web of Science][Medline]

McIntyre CC, Grill WM, Sherman DL, Thakor NV. Cellular effects of deep brain stimulation: model-based analysis of activation and inhibition. J Neurophysiol 91: 1457–1469, 2004.[Abstract/Free Full Text]

Meltzer NE, Ryugo DK. Projections from auditory cortex to cochlear nucleus: a comparative analysis of rat and mouse. Anat Rec A Discov Mol Cell Evol Biol 288: 397–408, 2006.[Medline]

Miyashita Y. Cognitive memory: cellular and network machineries and their top-down control. Science 306: 435–440, 2004.[Abstract/Free Full Text]

Mooney DM, Zhang L, Basile C, Senatorov VV, Ngsee J, Omar A, Hu B. Distinct forms of cholinergic modulation in parallel thalamic sensory pathways. Proc Natl Acad Sci USA 101: 320–324, 2004.[Abstract/Free Full Text]

Naya Y, Yoshida M, Miyashita Y. Backward spreading of memory-retrieval signal in the primate temporal cortex. Science 291: 661–664, 2001.[Abstract/Free Full Text]

Nowak LG, Bullier J. Spread of stimulating current in the cortical grey matter of rat visual cortex studied on a new in vitro slice preparation. J Neurosci Methods 67: 237–248, 1996.[CrossRef][Web of Science][Medline]

O'Leary DD, Koester SE. Development of projection neuron types, axon pathways, and patterned connections of the mammalian cortex. Neuron 10: 991–1006, 1993.[CrossRef][Web of Science][Medline]

Palomero-Gallagher N, Zilles K. The Rat Nervous System (3rd ed.), edited by Paxinos G. Amsterdam: Elsevier Academic, 2004.

Parent M, Parent A. Single-axon tracing study of corticostriatal projections arising from primary motor cortex in primates. J Comp Neurol 496: 202–213, 2006.[CrossRef][Web of Science][Medline]

Paxinos G, Watson C. The Rat Brain: In Stereotaxic Coordinates. Sydney: Academic, 1986.

Pelletier JG, Pare D. Uniform range of conduction times from the lateral amygdala to distributed perirhinal sites. J Neurophysiol 87: 1213–1221, 2002.[Abstract/Free Full Text]

Ranck JB Jr. Which elements are excited in electrical stimulation of mammalian central nervous system: a review. Brain Res 98: 417–440, 1975.[CrossRef][Web of Science][Medline]

Rushton WA. A theory of the effects of fiber size in medullated nerve. J Physiol 115: 101–122, 1951.[Free Full Text]

Salami M, Itami C, Tsumoto T, Kimura F. Change of conduction velocity by regional myelination yields constant latency irrespective of distance between thalamus and cortex. Proc Natl Acad Sci USA 100: 6174–6179, 2003.[Abstract/Free Full Text]

Shi CJ, Cassell MD. Cortical, thalamic, and amygdaloid projections of rat temporal cortex. J Comp Neurol 382: 153–175, 1997.[CrossRef][Web of Science][Medline]

Shmiel T, Drori R, Shmiel O, Ben-Shaul Y, Nadasdy Z, Shemesh M, Teicher M, Abeles M. Neurons of the cerebral cortex exhibit precise interspike timing in correspondence to behavior. Proc Natl Acad Sci USA 102: 18655–18657, 2005.[Abstract/Free Full Text]

Sirota MG, Swadlow HA, Beloozerova IN. Three channels of corticothalamic communication during locomotion. J Neurosci 25: 5915–5925, 2005.[Abstract/Free Full Text]

Steinmetz PN, Roy A, Fitzgerald PJ, Hsiao SS, Johnson KO, Niebur E. Attention modulates synchronized neuronal firing in primate somatosensory cortex. Nature 404: 187–190, 2000.[CrossRef][Medline]

Sugihara I, Lang EJ, Llinas R. Uniform olivocerebellar conduction time underlies Purkinje cell complex spike synchronicity in the rat cerebellum. J Physiol 470: 243–271, 1993.[Abstract/Free Full Text]

Swadlow HA. Efferent neurons and suspected interneurons in second somatosensory cortex of the awake rabbit: receptive fields and axonal properties. J Neurophysiol 66: 1392–1409, 1991.[Abstract/Free Full Text]

Swadlow HA. Neocortical efferent neurons with very slowly conducting axons: strategies for reliable antidromic identification. J Neurosci Methods 79: 131–141, 1998.[CrossRef][Web of Science][Medline]

Swadlow HA, Weyand TG. Efferent systems of the rabbit visual cortex: laminar distribution of the cells of origin, axonal conduction velocities, and identification of axonal branches. J Comp Neurol 203: 799–822, 1981.[CrossRef][Web of Science][Medline]

Vaudano E, Legg CR, Glickstein M. Afferent and efferent connections of temporal association cortex in the rat: a horseradish peroxidase study. Eur J Neurosci 3: 317–330, 1991.[CrossRef][Web of Science][Medline]

Waxman SG, Bennett MV. Relative conduction velocities of small myelinated and non-myelinated fibres in the central nervous system. Nat New Biol 238: 217–219, 1972.[CrossRef][Web of Science][Medline]

Waxman SG, Swadlow HA. The conduction properties of axons in central white matter. Prog Neurobiol 8: 297–324, 1977.[CrossRef][Web of Science][Medline]

Winer JA, Diehl JJ, Larue DT. Projections of auditory cortex to the medial geniculate body of the cat. J Comp Neurol 430: 27–55, 2001.[CrossRef][Web of Science][Medline]

Winkowski DE, Knudsen EI. Top-down gain control of the auditory space map by gaze control circuitry in the barn owl. Nature 439: 336–339, 2006.[CrossRef][Medline]





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Figures
Right arrow All Versions of this Article:
100/1/327    most recent
90392.2008v2
90392.2008v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via Web of Science (2)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Chomiak, T.
Right arrow Articles by Hu, B.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Chomiak, T.
Right arrow Articles by Hu, B.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2008 by the The American Physiological Society.