|
|
||||||||
Department of Neurobiology and Center for the Neural Basis of Cognition, University of Pittsburgh School of Medicine, Pittsburgh, Pennsylvania 15261
Submitted 6 May 2004; accepted in final form 16 June 2004
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
To develop a generic view of ganglionic organization, we began by considering paravertebral sympathetic ganglia 9 and 10 of the bullfrog. This model system has many features that facilitate experimental study (Smith 1994
). These advantages include 1) low levels of synaptic convergence that can be readily assessed, 2) functional criteria that permit the physiological identification and independent analysis of B and C neurons, which selectively innervate cutaneous glands and blood vessels, 3) well-defined modulatory mechanisms of muscarinic and peptidergic co-transmission, and 4) established methods that permit the study of fully differentiated neurons from dissociated ganglia in adult animals. In B-type neurons, Karila and Horn (2000)
have identified an n + 1 convergence pattern of weak and strong nicotinic synapses, where every postganglionic neuron receives one very strong, always suprathreshold, synapse and a variable number (n) of weak synapses. Arguing that a similar pattern occurs in the paravertebral ganglia of all vertebrates, they derived a mathematical theory of ganglionic integration based on the principle of stochastic coincidence detection. This approach led to the broad prediction that, paravertebral sympathetic ganglia, in general, behave as use-dependent synaptic amplifiers of preganglionic activity. In this scheme, gain is regulated by the number of converging nicotinic synapses and by their strength and plasticity.
Although the coincidence detection theory of Karila and Horn (2000)
provided a simple means for understanding how seemingly unrelated features of ganglionic organization could interact to produce amplification, it was limited by a lack of biophysical details describing ionic conductances and determinants of synaptic plasticity. Schobesberger et al. (1999
, 2000
) partially addressed this problem by adapting an earlier conductance-based model of the sympathetic B-type neuron (Yamada et al. 1989
) and examining the interaction between slow metabotropic excitation and subthreshold nicotinic excitatory postsynaptic potentials (EPSPs). The resulting analysis showed that muscarinic and peptidergic slow EPSPs could act to strengthen nicotinic synapses, but did not directly test the consequences for synaptic gain.
This study examines ganglionic integration by combining simulations of ganglionic activity with dynamic-clamp experiments. The simulations reveal how different patterns of preganglionic activity, in conjunction with synaptic mechanisms, control synaptic gain. In addition to clarifying certain limitations of the more minimal coincidence detection theory (Karila and Horn 2000
), the present simulations establish the range of synaptic transformations of activity that may be implemented by different populations of autonomic neurons. In whole cell recordings from dissociated sympathetic B neurons, the dynamic-clamp method has been used to drive postsynaptic firing with different patterns of virtual synaptic activity. By taking this approach, one has precise control over the number, strength, and occurrence of each synaptic event and can compare the response of different cells to identical stimulation, an experimental procedure that would not be possible with living synapses. The experimental results confirm key predictions of the simulations and identify a postsynaptic role for calcium in regulating synaptic gain.
| METHODS |
|---|
|
|
|---|
The part of the autonomic system we are studyingsympathetic and parasympathetic gangliaconsists of a population of spinal, preganglionic presynaptic neurons connected via a simple network of synapses to a population of purely postsynaptic ganglionic neurons. As in previous work (Karila and Horn 2000
), synaptic gain (G) has been defined as the multiplicative factor that relates the system's spike output to its input, where each is represented by their average firing rates (fpre, fpost). Gain in this system has a simple physical interpretation. When G > 1, the postsynaptic neurons, on average, are firing faster than the presynaptic neurons
![]() | (1) |
To simplify the task of simulating ganglionic integration, three central postulates were adopted from earlier work (Karila and Horn 2000
). Each postulate represents an idealization of experimental observations.
POSTULATE 1VARIATIONS IN CONVERGENCE FOLLOW AN N + 1 RULE IN SYMPATHETIC GANGLIA.
Under this rule, every postganglionic neuron receives one very strong, always suprathreshold, nicotinic synapse denoted as primary and a variable number (n) of weak nicotinic synapses denoted as secondary. In many instances, summation between pairs of secondary EPSPs is sufficient to trigger an action potential. Due to this synaptic arrangement, the spike output of postganglionic neurons reflects the activity of primary synapses combined with that driven by suprathreshold coincidences between subthreshold secondary EPSPs. The n + 1 rule describes nicotinic convergence in bullfrog sympathetic ganglia (Dodd and Horn 1983
; Karila and Horn 2000
) and also seems to operate in mammalian sympathetic ganglia (Hirst and McLachlan 1986
; Jänig and McLachlan 1992
; Skok and Ivanov 1983
).
POSTULATE 2ENTIRE POPULATIONS OF AUTONOMIC NEURONS ARE UNIFORM IN THEIR SYNAPTIC ORGANIZATION.
Functional subsets of ganglionic neurons can be differentiated by their levels of secondary synaptic convergence (n) and by their selective expression of neuromodulatory mechanisms. For example, bullfrog secretomotor B-type sympathetic neurons have between one and three secondary synapses with a mean n of 1.2 (Karila and Horn 2000
), and they selectively express an excitatory muscarinic synaptic mechanism (Smith and Weight 1986
; Tosaka et al. 1968
; Weight and Padjen 1973
). As a consequence of uniformity within the B-cell population and the general lack of lateral interactions between postganglionic neurons, one can infer population behavior of this cell type from that of a single neuron, thus simplifying the design of simulations and experiments. Using this approach, it also becomes possible to compare integration in subpopulations of sympathetic and parasympathetic neurons that differ in their synaptic organization. Such variations appear to reflect the phenotypic specialization of ganglionic neuronal subtypes that control different target tissues and the evolutionary adaptation of homologous ganglia from different vertebrate species (Gibbins 1995
; Gibbins et al. 2000
; Purves et al. 1986
).
POSTULATE 3PREGANGLIONIC NEURONAL ACTIVITY IS GENERATED BY A pOISSON PROCESS IN WHICH THE TIME INTERVALS BETWEEN SYNAPTIC EVENTS ARE EXPONENTIALLY DISTRIBUTED.
Postganglionic firing in anesthetized animals and awake humans is characteristically very irregular (Macefield et al. 2002
; McAllen and Malpas 1997
). In cases where the pattern of presynaptic activity has been recorded intracellularly from postganglionic neurons in vivo, the noisy intervals between synaptic events can be approximated by an exponential distribution (see Fig. 5A2 in McLachlan et al. 1998
and Fig. 4A in Häbler et al. 1999
). By idealizing presynaptic activity as a Poisson process, one can mathematically model the occurrence of synaptic events.
|
|
All simulations of synaptic activity followed the multi-step process outlined in Fig. 1A. In this procedure, templates of synaptic conductance were constructed by first determining the timing of primary and secondary nicotinic events. For each event, an appropriately sized postsynaptic conductance waveform was added to the template. Finally, the completed template was used to drive the conductance-based model sympathetic neuron by numerically integrating the underlying system of ordinary differential equations.
|
Psyn, an event was added to the template. For a Poisson process governing the timing of n identical synapses, Psyn(t
twin) =1 exp(nfpretwin) (Colquhoun 1971
The time course of the fast nicotinic synaptic conductance (gsyn) was derived by fitting a synaptic current, recorded from a B neuron under two-electrode voltage clamp, to the sum of two exponentials and scaling it to a peak amplitude of 1 (see Fig. 2A in Schobesberger et al. 2000
). The time constants were 1 ms for the rise and 5 ms for the decay. Unless stated otherwise, the strengths of nicotinic synapses were adjusted by scaling the unitary template to
syn = 100 nS for primary synapses and 9.61 nS for secondary synapses. The latter represents 90% of threshold-gsyn, the nicotinic conductance required to depolarize the resting model to threshold for generating an action potential (Schobesberger et al. 2000
).
|
In one set of simulations, presynaptic facilitation of nicotinic synapses was modeled as an additive process followed by exponential decay. In the presence of facilitation, each presynaptic event transiently increased the amplitudes (A) of all successive events by a factor having the following form: A =1 + 2 exp(t/
decay). The magnitude of facilitation (2) and the range of
decay (1250 ms) were chosen to reproduce experimental data from bullfrog B neurons (Shen and Horn 1995
).
Software implementation of simulations
Computational simulations employed programs written in MATLAB R1113 (MathWorks, Natick, MA) using ode15s, a stiff variable-step differential-equation solver, and a personal computer. Our software, collectively called NEUROSIM, is freely available and can be obtained by contacting the authors. The voltage data and activation/inactivation parameter values from each simulation were saved every 0.25 ms. Simulating 204,000 s of neural activity for each measurement of synaptic gain required between 2 and 100 h of real time. Illustrations were drafted with Igor Pro 4.0 (PC Edition, WaveMetrics, Lake Oswego, OR).
Cell culture
Paravertebral sympathetic ganglia 810 were dissected from pithed adult bullfrogs (Rana catesbeiana, 46 in) using a protocol approved by the University of Pittsburgh's Institutional Animal Care and Use Committee. The ganglia were desheathed, minced, and dissociated by incubation in Ringer containing 300 µg/ml Liberase Blendzyme 3 (Roche Diagnostics, Indianapolis, IN) for 3060 min at 32°C, followed by mechanical trituration with graded fire-polished Pasteur pipettes. Ringer contained (in mM) 115 NaCl, 2 KCl, 1.8 CaCl2, and 4 NaHEPES, adjusted to pH 7.3. After centrifugation and resuspension in culture medium, the dissociated cells were plated onto glass coverslips coated with poly-D-lysine (BD Biosciences, Bedford, MA). Cell cultures were maintained for
2 wk in diluted L-15 medium (65%) supplemented with 5 mM dextrose, 2 mM L-glutamine, 1 mM CaCl2, 5% fetal bovine serum, 100 U/ml penicillin, and 100 µg/ml streptomycin.
Perforated-patch clamp recording and dynamic clamp
Whole cell perforated-patch recordings were made at room temperature from neurons visualized on an inverted microscope. Fire-polished patch pipettes were filled with internal solution containing (in mM) 110 K+-gluconate, 10 NaCl, and 5 NaHEPES, adjusted to pH 7.1 with HCl. Stock solution of amphotericin-B (Sigma; 25 mg/ml in dimethyl sulfoxide) was added to the internal solution to a final concentration of 250 µg/ml. With pipette resistances of 25 M
, this amphotericin-B concentration typically resulted in final access resistances of 510 M
, after 3060 min of recording. Access resistance was monitored throughout the experiment and compensated using the bridge circuitry of the recording amplifier. The external recording solution consisted of Ringer as described above.
The dynamic current-clamp method was used to record membrane potential and create virtual nicotinic synapses. For this purpose, we used G-clamp (Kullmann et al. 2004
), which includes an AxoClamp 2B amplifier (Axon Instruments, Union City, CA), G-clamp software written in the LabVIEW-RT programming environment (National Instruments, Austin, TX), a standard Windows-based computer, an embedded controller running under a real-time operating system (National Instruments), and a multifunction I/O board for data acquisition and signal generation (PXI-6052E, National Instruments). This system has the ability to read data files containing conductance templates of synaptic activity identical to those used in the numerical simulations (Fig. 1A). These templates were used to stimulate neurons while measuring and storing their membrane potential responses and the currents injected by the clamp. The dynamic clamp's real-time loop speed and data sampling rate were both 20 kHz.
Recordings were made from the largest cells because they are generally B neurons (Dodd and Horn 1983
; Jan and Jan 1982
). Once patch perforation by amphotericin-B had reduced access resistance to 15 M
or less, we began measuring threshold-gsyn. When a stable estimate was obtained, this value was used to scale a synaptic conductance template for measuring synaptic gain. Because execution of a synaptic conductance template could transiently alter postsynaptic excitability for several minutes, we continually measured threshold-gsyn between trials, waiting until it stabilized again before execution of the next synaptic conductance template. Steady-state I-V relations were measured repeatedly throughout each experiment to monitor leak resistance.
Analysis of synaptic gain and the associated errors
Synaptic gain (G) was calculated by dividing the number of postsynaptic action potentials in a given trial by the average number of presynaptic events per synapse. To estimate the errors associated with counts of postsynaptic action potentials and with estimates of G, we assumed that postsynaptic activity was generated by an inhomogeneous Poisson process (Berry and Meister 1998
). The inhomogeneity arose from filtering of short interspike intervals at high stimulus frequencies (Fig. 3A). For all Poisson processes, the mean number of events (µ) equals the total number of events generated, and the CV = 1/
µ (Colquhoun 1971
). We therefore adjusted the duration of simulations to produce
400 postsynaptic action potentials to obtain a CV or "counting error" of
5%. The error in G, which depends on the count of presynaptic and postsynaptic spikes, was estimated as CVG
(Colquhoun 1971
). The SD of G can be calculated as
G = CVGG. Since a single simulation was used to calculate each value of G, the SE equals
G. All errors in G are plotted or reported as ±SE.
|
| RESULTS |
|---|
|
|
|---|
To examine how the presynaptic firing frequency controls gain, templates of synaptic activity were constructed and used to stimulate a model sympathetic neuron (Fig. 1A). Each synaptic template incorporated information about the number of converging nicotinic synapses, their strength, and their temporal pattern of activity. We began by considering a case with one primary synapse and nine secondary synapses, an n + 1 arrangement similar to that observed in the rat superior cervical ganglion (Purves et al. 1986
). The strength of each secondary nicotinic synapse was set to 9.61 nS, which is 90% of threshold-gsyn. The temporal pattern of synaptic events was created by assuming presynaptic activity arises from a Poisson process. By combining the timing information with appropriately scaled synaptic conductance waveforms, a complete template was constructed that described the time course of total synaptic conductance over an entire simulation. The membrane potential response to synaptic stimulation was computed by numerically integrating the Hodgkin-Huxley system of equations (see supplemental material for a full listing) that describe the model sympathetic neuron (Fig. 1, B and C).
The efficacy of secondary synapses and their ability to generate gain varied with the presynaptic frequency. In response to stimulation at 0.3 Hz, most postsynaptic action potentials were driven in one-to-one fashion by primary EPSPs (Fig. 1B). In contrast, setting the presynaptic frequency to 3 Hz changed the outcome so that most postsynaptic firing was driven by secondary synapses (Fig. 1C). Inspection of the records revealed that virtually all of the spikes triggered by secondary EPSPs resulted from binary coincidences. At 3 Hz, there were occasional triple and quadruple coincidences, but they had little impact on the overall rate of postsynaptic firing (Fig. 1C). These effects of activity on synaptic gain can be understood more easily by considering a simpler case with only three synapses (Fig. 2A). In this schematic illustration, all three presynaptic inputs fire an average of four times, and the postsynaptic cell responds by firing seven times, thus resulting in a gain of 1.75. Four of the action potentials (dark shading) originate from the primary synapse, while the remaining three (light shading) arise from temporal summation between secondary EPSPs. If the presynaptic frequency were increased in this example, the probability of secondary synaptic coincidences would also increase, eventually accounting for the majority of postsynaptic action potentials and a further increase in gain. On the other hand, reducing the presynaptic frequency would lessen the influence of secondary synapses, thereby causing gain to drop toward one.
To measure synaptic gain systematically as a function of mean presynaptic frequency, we drove the 9 + 1 model over a range of frequencies (0.120 Hz). The durations of these simulations were each adjusted to produce
400 postsynaptic spikes to reduce counting errors arising from the stochastic nature of the synaptic activity. The resulting relation shows that gain was a biphasic function of the mean presynaptic frequency (Fig. 2B), with a peak in the range of 25 Hz. The overall shape of the synaptic-gain relation constructed from the simulations was similar to that predicted by the more sparsely framed coincidence detection theory (Karila and Horn 2000
), especially when the presynaptic frequency was <1 Hz (Fig. 2B). The main differences were observed at higher frequencies where the theory overestimated the peak gain and predicted an overly steep postpeak decline. We suspected these disparities were due to the different ways in which the two approaches introduce refractory periods that limit the rate of postsynaptic firing. For this purpose, the coincidence detection theory of Karila and Horn (2000)
incorporates the arbitrary assumption that secondary EPSPs can only drive postsynaptic firing at a maximum rate of 10 Hz. In contrast, refractory properties of the model neuron emerge from its conductances and are not limited to spikes triggered by secondary EPSPs. Lowering the refractory cut-off frequency in the Karila and Horn (2000)
theory reduced the peak of the gain relation but did not improve the overall fit to the simulation results (data not shown). To better understand the behavior of the conductance-based model, we therefore examined how primary and secondary EPSPs affect one another.
Interactions between primary and secondary synapses have consequences for gain
In addition to limiting the excitation produced by the summation of secondary EPSPs, one would expect refractory periods to introduce interactions between primary and secondary synapses. If, for example, a pair of secondary EPSPs coincided during or shortly after an action potential driven by a primary EPSP, the first action potential would prevent the second. To measure the extent of such effects, we first constructed separate gain relations to describe a single primary synapse and a group of nine secondary synapses (Fig. 2C). The gain of the solitary primary synapse was unitary at low stimulation rates, but decreased by
17% when the presynaptic frequency was >10 Hz. This high-frequency roll-off suggests that refractory interactions between primary synaptic events can limit postsynaptic firing. In contrast, the gain relation for nine secondary synapses was biphasic, with a baseline that approached zero at lower frequencies, a peak of
2.1 at 310 Hz, and a roll-off at higher rates. The high-frequency roll-off in this relation suggests that refractory mechanisms can also limit the efficacy of secondary synapses in the absence of a primary synapse.
To uncover interactions between primary and secondary synapses, we added the gain responses for isolated primary and secondary synapses and compared their sum with the gain relation describing combined n + 1 convergence (Fig. 2D). Although the two curves were similar when fpre < 1 Hz, the summed gain response was greater than the combined gain response at higher frequencies. This indicates that interactions between primary and secondary synapses are minimal at low fpre, but limit gain at higher stimulation rates.
Refractory periods alter postsynaptic spike timing
If postsynaptic refractoriness is responsible for the roll-off in synaptic gain at higher presynaptic frequencies (Fig. 2, BD), frequency-dependent effects should be evident in the patterns of postsynaptic action potentials. This prediction was tested by constructing cumulative probability distributions of postsynaptic interspike intervals (Fig. 3A). When the presynaptic frequency was 0.6 Hz and gain was low, the distribution of postsynaptic interspike intervals was exponential, just like the distribution of presynaptic events. However, at higher presynaptic frequencies, the distributions of postsynaptic intervals fell short of the expected exponential relations. The disparity was clearest when the presynaptic firing rate was 20 Hz, which produced a large shortfall in postsynaptic interspike intervals <30 ms, with none <8 ms in duration. The lack of intervals <8 ms is also evident at the other presynaptic frequencies when the time scale of the distributions is expanded (Fig. 3A, insets).
To test whether refractory effects could account for the postsynaptic filtering effect, we measured two types of interactions. For pairs of primary EPSPs, the refractory period for production of an overshooting action potential was 8 ms (Fig. 3B). When the trailing stimulus was created by a pair of simultaneous secondary EPSPs, which was much weaker than a primary EPSP, the refractory period was 26 ms (Fig. 3C). The brief refractory period between primary EPSPs explains the total absence of postsynaptic interspike intervals <8 ms (Fig. 3A), and it accounts for the roll-off in the gain relation of the isolated primary synapse (Fig. 2C). When the primary synapse was driven at 20 Hz, 13% of the intervals between EPSPs were <8 ms, which agrees well with the observed 17% reduction in gain. In contrast, the longer refractory period between primary and secondary EPSPs explains the partial filtering of postsynaptic interspike intervals <30 ms (Fig. 3A) and the roll-off of the gain relation for secondary synapses (Fig. 2C).
In the absence of refractory mechanisms, one would expect gain to increase as a monotonic function of the mean presynaptic frequency, saturating when gain is n + 1 (Karila and Horn 2000
). We tested this prediction by counting the number of EPSPs separated by intervals that were <26 ms in the secondaries-only 20-Hz data. If all of these EPSPs had successfully triggered spikes, the gain would have been 9 rather than the observed 1.2. This agrees with the theoretical maximum for gain in the secondaries-only configuration, thereby supporting the conclusion that the refractory occlusion by secondary EPSPs causes most of the roll-off in the synaptic gain relation. In comparison, the effects of primary-primary synaptic occlusion (Fig. 2C) and primary-secondary synaptic occlusion (Fig. 2D) were relatively small.
Synaptic gain regulation by convergence and synaptic strength
Lowering secondary synaptic convergence to mimic the levels seen in sympathetic ganglia of amphibians and mice (n = 1.2, 3) (Karila and Horn 2000
; Purves et al. 1986
) caused a proportional reduction in maximum gain (Fig. 4A). Reducing secondary synaptic strength from 98 to 50% threshold-gsyn, while holding convergence constant at n = 9, showed that the peak of the gain relation also varied as a function of secondary synaptic strength (Fig. 4B).
Although the effects of convergence and synaptic strength on gain were generally consistent with the predictions of coincidence detection theory (Karila and Horn 2000
), it was surprising to find G > 1 when secondary synaptic strength was set to 50% threshold-gsyn (Fig. 4B). This result contradicts the theory, which would require perfect coincidences between pairs of secondary EPSPs of this strength to generate action potentials. Since perfect coincidences are rare, the theory predicts there should be little or no gain above the baseline of 1 (Fig. 4B, dashed line). We examined this further by using even weaker synapses. For fpre = 5 Hz and n = 9, G = 1.31 ± 0.06 when secondary synaptic strength was 49% threshold-gsyn and G = 1.17 ± 0.05 when secondary synaptic strength was 40% threshold-gsyn. The ability of sub-50% secondary EPSPs to trigger action potentials suggested it was very unlikely that binary coincidences were causing these effects. By examining the simulation records for triple, quadruple, and quintuple coincidences, we found that such higher order events accounted for 9799% of the spikes driven by secondary synapses at 50, 49, and 40% threshold-gsyn.
Metabotropic modulation of synaptic gain
Metabotropic muscarinic excitation of sympathetic neurons operates through a branched signaling pathway that controls two ionic conductances. The better known is the M-type potassium conductance gM, whose suppression during synaptic stimulation of muscarinic receptors by acetylcholine produces a slow EPSP that can last for tens of seconds (Adams et al. 1982
; Smith 1994
; Tosaka et al. 1968
). In addition, muscarinic receptors activate a second membrane current that arises from a nonselective cation conductance (Tsuji and Kuba 1988
). Schobesberger et al. (2000)
suggested that the cation conductance is gated by cyclic nucleotides (gCNG) and went on to show through simulations that gCNG and gM can act alone or together to increase the effective strength of subthreshold nicotinic EPSPs. These findings imply that metabotropic excitation working through this and analogous pathways could serve to increase synaptic gain.
We tested this idea in simulations with nine secondary synapses set to a relatively weak resting level of 50% threshold-gsyn. Since the kinetics of muscarinic EPSPs are much slower than the time scale of nicotinic EPSPs and action potentials, metabotropic excitation was introduced by simply adjusting
M from 40 to 30 nS and gCNG from 0 to 4 nS. These modest conductance changes reflect the effects seen during muscarinic responses in sympathetic neurons (Tsuji and Kuba 1988
). Altering
M or gCNG alone produced clear increases in synaptic gain at all but the very lowest frequencies of presynaptic stimulation (Fig. 4C). Changing both conductances together produced even larger increases. Comparing the sum of the individual effects with the combined responses reveals that the two components add linearly at lower presynaptic frequencies but combine synergistically between 2 and 10 Hz. These results are consistent with earlier experiments showing that gM and gCNG act in a supralinear way to alter the strength of nicotinic synapses and to increase the temporal window of summation between pairs of subthreshold EPSPs (Schobesberger et al. 2000
).
Presynaptic facilitation increases synaptic gain
Mechanisms of presynaptic facilitation are widespread in autonomic ganglia. When facilitation was incorporated into synaptic templates, it increased synaptic gain at frequencies of 110 Hz (Fig. 4D). The exponential rate of decay of the presynaptic facilitation factor directly affected the magnitude of the gain response.
Experimental tests of computational predictions
The simulation results in the first part of this paper show how synaptic gain can arise as a use-dependent consequence of nicotinic synaptic convergence and be modulated by synaptic strength, metabotropic excitation, and presynaptic facilitation (Fig. 4). Evaluating the physiological significance of these findings requires that appropriate experimental tests be devised. Ideally, one would want to activate independently each of the synaptic inputs to an autonomic neuron with trains of asynchronous stimuli. Due to the impracticality of such an experiment, we adopted an alternative strategy of using the dynamic-clamp method to implement virtual nicotinic synapses on dissociated bullfrog B neurons. With this approach, one can compare how real B neurons and the model B neuron respond to identical patterns of synaptic stimulation. We began by measuring synaptic strength and the gain produced by n + 1 convergence.
Dynamic-clamp estimates of synaptic strength and gain
An automated binary search paradigm was used to measure the magnitude of virtual nicotinic synaptic conductance required to depolarize B neurons to threshold (Fig. 5). The procedure begins by testing boundary values that straddle threshold and homes in on a stable result within 812 trials (Kullmann et al. 2004
). In a group of 12 neurons that were subsequently used for measurements of gain, we observed threshold-gsyn = 14.7 ± 1.7 nS (range, 7.023.7 nS), Vrest = 54.7 ± 3.0 mV (range, 40 to 70 mV), and Rleak = 739 ± 111 M
(range, 153 M
1.48 G
). All of these cells had overshooting actions potentials. In comparison, the computational model assumed Rleak = 333 M
and predicted that threshold-gsyn = 10.68 nS and Vrest = 62.7 mV.
To measure the influence of secondary nicotinic EPSPs on synaptic gain, we constructed 5-Hz templates that were 40 s in duration and contained one primary synapse and zero, three, or nine secondary synapses. Five hertz was chosen as the presynaptic frequency because it generated maximal responses in gain simulations (Fig. 4A). After recording a value of threshold-gsyn that remained stable with repeated tests, the strength of secondary synapses in the stimulus templates was scaled to 90% threshold-gsyn, and the strength of the primary synapse was set at 10 times that of the secondary synapses. In cases where conditions drifted slowly over time, synaptic conductances were adjusted to maintain constant strength. Only cells in which at least two different levels of secondary convergence could be compared were included in the analysis. Figure 6 shows an example comparing convergence levels of three and nine secondary synapses. The top traces show dynamic-clamp recordings from one cell, and the bottom traces show simulated responses to the same synaptic templates using the computational model. From this example, three key points emerge. First, secondary synapses contributed to the spike output of the cell at both levels of convergence. As consequence of this effect, the gain produced in both cases was >1. Second, the contribution of secondary synapses was greater when n = 9 than when n = 3. This effect corresponds to greater convergence producing higher gain. Third, the pattern of spike output recorded from the B neuron was very well predicted by the computational simulations. In this example, only one spike is missing from the simulated output when n = 3 (Fig. 6A), and one spike is missing from the cellular output when n = 9 (Fig. 6B).
|
|
Activity and calcium dependence of synaptic gain
Finding that gain could vary with resting potential (Fig. 7B) suggested that it might also change with activity. To test for such effects, the responses to 5-Hz stimulation were analyzed as four 10-s segments. This revealed that synaptic gain decreased over the course of individual trials when n = 9, but not when n = 3, and not in simulations. The decline was reversible when cells were given a few minutes rest between trials. When dynamic-clamp records and simulations from the last 10 s of stimulus trials were compared, secondary synapses failed to drive many of the spikes that were predicted by the simulations. This outcome differed from the close correspondence between recorded and simulated firing patterns at the beginning of the stimulus trial (Fig. 6B).
To document the use-dependent decline in gain more carefully and to examine its possible origin, we designed an experiment using a 40-s synaptic template (fpre = 5 Hz, n = 9) consisting of 10 identical 4-s repeats. In a group of seven neurons (Vrest = 63.1 ± 3.1 mV, threshold-gsyn = 10.6 ± 2.1 nS, Rleak = 732 ± 84 M
), the repeating template yielded an average gain of 1.64 ± 0.13 over the entire 40-s period. This was similar to the gain observed with the nonrepeating template (1.42), but lower than the gain (2.40 ± 0.19) produced in a simulation with the same template. The gain generated during each repeat declined gradually from a high of 2.00 ± 0.21 in the first 4 s to a low of 1.37 ± 0.09 in the last 4 s (7 cells). The time course of the decline in gain was similar when viewed for one cell (Fig. 8A) and for the group (Fig. 8B). By contrast, simulated voltage responses to the first and last segments of the repeating template were identical, with a gain of 2.40 ± 0.26.
|
; P = 0.72, paired t-test), but caused significant increases in threshold-gsyn (13.3 ± 2.6 nS; P = 0.01, paired t-test) and in the average gain over the entire 40-s template (1.94 ± 0.17; P < 0.05, paired t-test; Fig. 8C). The overall changes in gain were also reflected in the 4-s repeats. Cadmium increased gain in the first 4-s period to 2.32 ± 0.19 and in the last 4-s period to 1.76 ± 0.16. These observations suggest that increased intracellular [Ca2+] may be partially responsible for the reduction in synaptic gain seen during prolonged synaptic stimulation. | DISCUSSION |
|---|
|
|
|---|
By simulating ganglionic integration with a conductance-based model, we have shown how the stereotypic n + 1 convergence pattern of fast nicotinic synapses found in paravertebral sympathetic ganglia can serve to create a low-pass filter with unitary baseline gain and a high-frequency boost (Figs. 1, 2, and 9A). Plastic changes in several factors including convergence (Fig. 4A), secondary synaptic strength (Fig. 4B), metabotropic excitation (Fig. 4C), and presynaptic facilitation (Fig. 4D) are all capable of enhancing the filter's high-frequency boost. These findings show how the diversity in synaptic organization reported for subclasses of sympathetic neurons (Jänig 1995
; Smith 1994
) may provide a basis for specialized regulation of ganglionic gain. For example, vasomotor, secretomotor, and pilomotor sympathetic neurons may each have unique synaptic-gain functions that are optimized to drive their specific targets.
|
The architecture of the conductance-based simulations followed from a much simpler coincidence-detection theory (Karila and Horn 2000
) that rests on two basic ideaspreganglionic activity is noisy and pairs of secondary EPSPs can sum to reach threshold. Perhaps it is not surprising then that the theory successfully predicted the general form of gain relations produced by simulating n + 1 convergence and by simulating modulation through different types of synaptic plasticity (Figs. 2, B and C, and 4). Nonetheless, the theory consistently erred by predicting an overly sharp peak in the gain relation and an overly steep roll-off at high presynaptic frequencies. These effects arise from the theory's assumption that secondary synapses cannot drive firing at rates >10 Hz. In contrast to this hard cut-off, more complex refractory interactions emerge in the conductance-based model to control the peak value of gain and its high-frequency roll-off (Fig. 3). The net effect of these refractory interactions is to produce a more gradual filtering of secondary synaptic events as a function of the presynaptic firing rate. In addition, the simulations with weak secondary synapses have shown that higher-order coincidences can influence firing (Fig. 4B). Using refractory filters having arbitrary sigmoidal shapes rather than step functions, we were able to fit closely the gain curves generated by simulations of the conductance-based model (data not shown). However, we were unable to modify the theory to incorporate these mechanisms in a form that could be derived on mechanistic grounds and that would readily generalize to different conditions. One might expect such limitations to emerge, given that the underlying ordinary differential equations, based on the Hodgkin-Huxley formalism, are rarely susceptible to simple analytical solutions. In other words, coincidence-detection theory concisely captures the basic form, but not the exact form, of simulation results from the conductance-based model.
Insights from dynamic-clamp experiments
Even though the conductance-based model was designed to include many biophysical mechanisms, it omitted some details to achieve computational efficiency and to minimize the inclusion of arbitrary assumptions. For this reason, it was important to compare the responses of living neurons to the computational results. By taking such an approach, one can critically evaluate the model and focus attention on mechanisms that need further scrutiny.
The dynamic-clamp data show that the model correctly predicted the magnitude threshold-gsyn for nicotinic excitation and the effect of secondary synaptic convergence on synaptic gain during 5-Hz stimulation. In contrast, the model did not predict the activity dependence of synaptic gain, as revealed at n = 9 by the apparent correlation with resting potential (Fig. 7B) and by the decline in gain during the course of single stimulus trials (Fig. 8, A and B). Treating cells with 10 µM Cd2+ raised the level of gain during the first 4 s of stimulation (fpre = 5 Hz, n = 9) from 2.00 ± 0.21 to 2.32 ± 0.19, which is close to the level of 2.40 ± 0.26 predicted by the model. Cadmium also appeared to decrease, but not eliminate, the decline of gain (Fig. 8A). Although these results indicate that calcium-dependent processes play a role in the regulation of gain, they cannot fully account for the use-dependent decline observed during 5-Hz stimulation and should therefore be interpreted with caution. Within our records, we also looked for, but were unable to detect obvious changes in spike threshold and morphology (upstroke, peak, repolarization) that could account for the decrease in gain during stimulation. Solving this problem will require further experimental analysis.
Function follows form: other autonomic ganglia and brain relays
Outside of the paravertebral sympathetic system, the different patterns of convergence seen in other ganglia can now be interpreted, based on our results. In the simplest pattern, postganglionic neurons are innervated by one primary nicotinic synapse. This arrangement occurs in the rat parasympathetic submandibular ganglion (Lichtman 1980
), which controls salivation, and in the frog parasympathetic cardiac ganglion (Dennis et al. 1971
), which modulates contractions of the heart. In our model, neurons with this innervation pattern behave as 1:1 relays that roll off at higher presynaptic frequencies (Figs. 2C and 9B). In another motif, characterized by the absence of a primary synapse and the presence of multiple secondary synapses, autonomic neurons behave as a band-pass filter (Figs. 2C and 9C). This organization enables the parasympathetic pelvic ganglion in the cat to perform as an ON-OFF switch during the micturition reflex (DeGroat and Saum 1976
), which gates filling and emptying of the bladder. A similar convergence pattern also exists in a subpopulation of neurons in the guinea-pig celiac ganglion (Gibbins et al. 2003
; McLachlan and Meckler 1989
). These prevertebral sympathetic neurons, which receive both preganglionic and intestinofugal inputs (Gibbins et al. 2003
), exert a modulatory gating influence on gastrointestinal motility, secretions, and blood flow (Furness and Costa 1987
).
The levels of convergence seen in autonomic ganglia are sparse compared with the massive convergence found on cortical pyramidal cells, cerebellar Purkinje cells, and spinal motoneurons. However, sparse convergence operates in many brain nuclei that are traditionally viewed, like autonomic ganglia, as synaptic relays. Perhaps the best example is the lateral geniculate nucleus (LGN). Simultaneous recordings between retinal ganglion cells and the LGN have revealed that LGN neurons receive excitatory input from fewer than six ganglion cells (Mastronarde 1992
), and these connections are characterized as strong or weak (Usrey et al. 1999
). Our approach may therefore have relevance to other circuits with sparse convergence.
Importance of gain in autonomic negative-feedback loops
Negative-feedback loops are critical elements of physiological control processes, including those mediated by the neuroendocrine and autonomic systems (Khoo 2000
). The functional gain that emerges as an inherent component of a negative-feedback loop leads to faster and more accurate homeostatic control (Scher 1989
; Scher et al. 1991
). An important example of this mechanism arises in the baroreceptor reflex, an important focal point for neural control of the circulation (Scher et al. 1991
). The reflex allows blood pressure to be maintained within a relatively narrow range despite changes in body posture and cardiac output. In the process, afferent feedback from sensory receptors regulates the output of motor activity in the sympathetic and parasympathetic systems, ultimately leading to changes in cardiac function and vascular tone. Despite extensive study of the baroreceptor reflex, the neural origin of its gain remains unclear. Much of the focus has centered on synapses in the nucleus tractus solitarius and the rostral ventral lateral medulla (Dampney 1994
; Sun 1995
), but relatively little attention has been paid to the potential contributions of synapses in the spinal cord, in sympathetic ganglia, and at neuroeffectors to the overall gain of the reflex. We propose that synaptic gain in the baroreceptor reflex may be distributed throughout much of the circuit, including sympathetic ganglia, rather than being confined to one or two anatomical loci.
In vivo evidence of synaptic gain in paravertebral sympathetic ganglia
Our results indicate that synaptic gain rises above one when secondary synapses help to drive postganglionic firing. Evidence for such behavior under physiological conditions comes from in vivo intracellular recordings in the mouse, guinea pig, hamster, rabbit, rat, and bullfrog (Häbler et al. 1999
; Ivanoff and Smith 1995
; Ivanov and Purves 1989
; McLachlan et al. 1997
, 1998
; Skok and Ivanov 1983
). The data from these studies of paravertebral sympathetic neurons provide evidence that primary and secondary nicotinic synapses generally form an n + 1 pattern of convergence. In several instances, the morphology of the spike afterhypolarization has been used to identify which class of synapses triggers individual action potentials. Using this approach, it is clear that secondary synapses play a role in driving postganglionic activity. For examples, see Fig. 1B in Skok and Ivanov (1983)
, Fig. 1A in McLachlan et al. (1998)
, Fig. 4 in Ivanoff and Smith (1995
), and Fig. 3 in McLachlan (2003)
. Additional evidence for synaptic gain comes from the observation that average rates of subthreshold and suprathreshold synaptic events are very similar in a range of mammals (see Table 2 in Ivanov and Purves 1989
). When convergence is taken into account, these data indicate that individual postganglionic neurons fire at higher frequencies than preganglionic neurons. When viewed in the context of our model (Fig. 2A), this means that amplification is a common integrative feature of sympathetic ganglia.
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
Present address of D. W. Wheeler: Max-Planck Institute for Brain Research, Department of Neurophysiology, Deutschordenstrasse 46, D-60528 Frankfurt/Main, Germany (E-mail: wheeler@mpih-frankfurt.mpg.de).
| FOOTNOTES |
|---|
1 The supplementary material is available online at http://jn.physiology.org/cgi/content/full/00470.2004/DC1.
Address for reprint requests and other correspondence: D. W. Wheeler, Dept. of Neurobiology, Univ. of Pittsburgh School of Medicine, E1440 Biomedical Science Tower, Pittsburgh, PA 15261 (E-mail: wheeler{at}mpih-frankfurt.mpg.de).
| REFERENCES |
|---|
|
|
|---|
Akasu T and Nishimura T. Synaptic transmission and function of parasympathetic ganglia. Prog Neurobiol 45: 459522, 1995.[CrossRef][Web of Science][Medline]
Berry MJ and Meister M. Refractoriness and neural precision. J Neurosci 18: 22002211, 1998.
Colquhoun D. Lectures on BiostatisticsAn Introduction to Statistics with Applications in Biology and Medicine. Oxford, UK: Clarendon Press, 1971.
Dampney RA. Functional organization of central pathways regulating the cardiovascular system. Physiol Rev 74: 323364, 1994.
DeGroat WC and Saum WR. Synaptic transmission in parasympathetic ganglia in the urinary bladder of the cat. J Physiol 256: 137158, 1976.
Dennis MJ, Harris AJ, and Kuffler SW. Synaptic transmission and its duplication by focally applied acetylcholine in parasympathetic neurons in the heart of the frog. Proc R Soc Lond B Biol Sci 177: 509539, 1971.[Medline]
Dodd J and Horn JP. A reclassification of B and C neurones in the ninth and tenth paravertebral sympathetic ganglia of the bullfrog. J Physiol 334: 255269, 1983.
Furness JB and Costa M. The Enteric Nervous System. Edinburgh: Churchill Livingstone, 1987.
Gibbins I. Chemical neuroanatomy of sympathetic ganglia. In: Autonomic Ganglia, edited by McLachlan EM. Australia: Harwood Academic Publishers, 1995, p. 73121.
Gibbins IL, Jobling P, Messenger JP, Teo EH, and Morris JL. Neuronal morphology and the synaptic organisation of sympathetic ganglia. J Auton Nerv Syst 81: 104109, 2000.[CrossRef][Web of Science][Medline]
Gibbins IL, Teo EH, Jobling P, and Morris JL. Synaptic density, convergence, and dendritic complexity of prevertebral sympathetic neurons. J Comp Neurol 455: 285298, 2003.[CrossRef][Web of Science][Medline]
Häbler HJ, McLachlan EM, Jamieson J, and Davies PJ. Synaptic responses evoked by lower urinary tract stimulation in superior cervical ganglion cells in the rat. J Urol 161: 16661671, 1999.[CrossRef][Web of Science][Medline]
Hirst GD and McLachlan EM. Development of dendritic calcium currents in ganglion cells of the rat lower lumbar sympathetic chain. J Physiol 377: 349368, 1986.
Ivanoff AY and Smith PA. In vivo activity of B- and C-neurones in the paravertebral sympathetic ganglia of the bullfrog. J Physiol 485: 797815, 1995.
Ivanov A and Purves D. Ongoing electrical activity of superior cervical ganglion cells in mammals of different size. J Comp Neurol 284: 398404, 1989.[CrossRef][Web of Science][Medline]
Iversen S, Iversen L, and Saper CB. The autonomic nervous system and the hypothalamus. In: Principles of Neural Science (4th ed.), edited by Kandel ER, Schwartz JH, and Jessell TM. New York: McGraw-Hill, 2000, p. 962981.
Jan LY and Jan YN. Peptidergic transmission in sympathetic ganglia of the frog. J Physiol 327: 219246, 1982.
Jänig W and McLachlan EM. Characteristics of function-specific pathways in the sympathetic nervous system. Trends Neurosci 15: 475481, 1992.[CrossRef][Web of Science][Medline]
Jänig WM. Ganglion transmission in vivo. In: Autonomic Ganglia, edited by McLachlan EM. Luxembourg City, Luxembourg: Harwood, 1995, p. 349395.
Jones SW and Marks TN. Calcium currents in bullfrog sympathetic neurons. I. Activation kinetics and pharmacology. J Gen Physiol 94: 151167, 1989.
Karila P and Horn JP. Secondary nicotinic synapses on sympathetic B neurons and their putative role in ganglionic amplification of activity. J Neurosci 20: 908918, 2000.
Khoo MCK. Physiological Control Systems. Piscataway, NJ: IEEE Press, 2000.
Kuba K and Koketsu K. Synaptic events in sympathetic ganglia. Prog Neurobiol 11: 77169, 1978.[CrossRef][Web of Science][Medline]
Kullmann PHM, Wheeler DW, Beacom J, and Horn JP. Implementation of a fast 16-bit dynamic clamp using LabVIEW-RT. J Neurophysiol 91: 542554, 2004.
Lichtman JW. On the predominantly single innervation of submandibular ganglion cells in the rat. J Physiol 302: 121130, 1980.
Macefield VG, Elam M, and Wallin BG. Firing properties of single postganglionic sympathetic neurones recorded in awake human subjects. Auton Neurosci 95: 146159, 2002.[CrossRef][Web of Science][Medline]
Mastronarde DN. Nonlagged relay cells and interneurons in the cat lateral geniculate nucleus: receptive-field properties and retinal inputs. Vis Neurosci 8: 407441, 1992.[Web of Science][Medline]
McAllen RM and Malpas SC. Sympathetic burst activity: characteristics and significance. Clin Exp Pharmacol Physiol 24: 791799, 1997.[Web of Science][Medline]
McLachlan EM. Transmission of signals through sympathetic gangliamodulation, integration or simply distribution? Acta Physiol Scand 177: 227235, 2003.[CrossRef][Web of Science][Medline]
McLachlan EM, Davies PJ, Häbler HJ, and Jamieson J. On-going and reflex synaptic events in rat superior cervical ganglion cells. J Physiol 501: 165181, 1997.
McLachlan EM, Häbler H-J, Jamieson J, and Davies PJ. Analysis of the periodicity of synaptic events in neurones in the superior cervical ganglion of anaesthetized rats. J Physiol 511: 461478, 1998.
McLachlan EM and Meckler RL. Characteristics of synaptic input to three classes of sympathetic neurone in the coeliac ganglion of the guinea-pig. J Physiol 415: 109129, 1989.
Powley TL. Central control of autonomic functions: organization of the autonomic nervous system. In: Fundamental Neuroscience (2nd ed.), edited by Squire LR, Bloom FE, McConnell SK, Roberts JL, Spitzer NC, and Zigmond MJ. San Diego: Academic Press, 2003, p. 911933.
Purves D, Rubin E, Snider WD, and Lichtman J. Relation of animal size to convergence, divergence, and neuronal number in peripheral sympathetic pathways. J Neurosci 6: 158163, 1986.[Abstract]
Scher AM. Cardiovascular control. In: Textbook of Physiology: Circulation, Respiration, Body Fluids, Metabolism, and Endocrinology, edited by Patton HD, Fuchs AF, Hille B, Scher AM, and Steiner R. Philadelphia, PA: Saunders Company, 1989, p. 972990.
Scher AM, O'Leary DS, and Sheriff DD. Arterial baroreceptor regulation of peripheral resistance and of cardiac performance. In: Baroreceptor Reflexes: Integrative Functions and Clinical Aspects, edited by Persson PB and Kirchheim HR. Berlin: Springer-Verlag, 1991, p. 75125.
Schobesberger H, Gutkin BS, and Horn JP. A minimal model for metabotropic modulation of fast synaptic transmission and firing properties in bullfrog sympathetic B neurons. Neurocomputing 26-27: 255262, 1999.[CrossRef]
Schobesberger H, Wheeler DW, and Horn JP. A model for pleiotropic muscarinic potentiation of fast synaptic transmission. J Neurophysiol 83: 19121923, 2000.
Shen WX and Horn JP. A presynaptic mechanism accounts for the differential block of nicotinic synapses on sympathetic B and C neurons by d-Tubocurarine. J Neurosci 15: 50255035, 1995.[Abstract]
Simmons MA. The complexity and diversity of synaptic transmission in the prevertebral sympathetic ganglia. Prog Neurobiol 24: 4393, 1985.[CrossRef][Web of Science][Medline]
Skok VI and Ivanov AY. What is the ongoing activity of sympathetic neurons? J Auton Nerv Syst 7: 263270, 1983.[CrossRef][Web of Science][Medline]
Smith PA. Amphibian sympathetic ganglia: an owner's and operator's manual. Prog Neurobiol 43: 439510, 1994.[CrossRef][Web of Science][Medline]
Smith PA and Weight FF. The pathway for the slow inhibitory postsynaptic potential in bullfrog sympathetic ganglia. J Neurophysiol 56: 823834, 1986.
Sun M-K. Central neural organization and control of sympathetic nervous system in mammals. Prog Neurobiol 47: 157233, 1995.[CrossRef][Web of Science][Medline]
Tosaka T, Chichibu S, and Libet B. Intracellular analysis of slow inhibitors and excitatory postsynaptic potentials in sympathetic ganglia of the frog. J Neurophysiol 31: 396409, 1968.
Tsuji S and Kuba K. Muscarinic regulation of two ionic currents in the bullfrog sympathetic neurone. Pfluegers 411: 361370, 1988.
Usrey WM, Reppas JB, and Reid RC. Specificity and strength of retinogeniculate connections. J Neurophysiol 82: 35273540, 1999.
Weight FF and Padjen A. Acetylcholine and slow synaptic inhibition in frog sympathetic ganglion cells. Brain Res 55: 225228, 1973.[CrossRef][Web of Science][Medline]
Yamada WM, Koch C, and Adams PR. Multiple channels and calcium dynamicsmodeling bullfrog sympathetic ganglion cells. In: Methods in Neuronal Modeling from Synapses to Networks, edited by Koch C and Segev I. Cambridge, MA: MIT Press, 1989, p. 97133.
This article has been cited by other articles:
![]() |
C. Li and J. P. Horn Differential Inhibition of Ca2+ Channels by {alpha}2-Adrenoceptors in Three Functional Subclasses of Rat Sympathetic Neurons J Neurophysiol, December 1, 2008; 100(6): 3055 - 3063. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. H. M. Kullmann and J. P. Horn Excitatory Muscarinic Modulation Strengthens Virtual Nicotinic Synapses on Sympathetic Neurons and Thereby Enhances Synaptic Gain J Neurophysiol, December 1, 2006; 96(6): 3104 - 3113. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. Zaika, L. S. Lara, N. Gamper, D. W. Hilgemann, D. B. Jaffe, and M. S. Shapiro Angiotensin II regulates neuronal excitability via phosphatidylinositol 4,5-bisphosphate-dependent modulation of Kv7 (M-type) K+ channels J. Physiol., August 15, 2006; 575(1): 49 - 67. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Li and J. P. Horn Physiological Classification of Sympathetic Neurons in the Rat Superior Cervical Ganglion J Neurophysiol, January 1, 2006; 95(1): 187 - 195. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. L Morris, I. L Gibbins, and P. Jobling Post-stimulus potentiation of transmission in pelvic ganglia enhances sympathetic dilatation of guinea-pig uterine artery in vitro J. Physiol., July 1, 2005; 566(1): 189 - 203. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |