Quantifying gene network connectivity in silico: scalability and accuracy of a modular approach.

Large, complex data sets that are generated from microarray experiments, create a need for systematic analysis techniques to unravel the underlying connectivity of gene regulatory networks. A modular approach, previously proposed by Kholodenko and co-workers, helps to scale down the network complexity into more computationally manageable entities called modules. A functional module includes a gene's mRNA, promoter and resulting products, thus encompassing a large set of interacting states. The essential elements of this approach are described in detail for a three-gene model network and later extended to a ten-gene model network, demonstrating scalability. The network architecture is identified by analysing in silico steady-state changes in the activities of only the module outputs, communicating intermediates, that result from specific perturbations applied to the network modules one at a time. These steady-state changes form the system response matrix, which is used to compute the network connectivity or network interaction map. By employing a known biochemical network, the accuracy of the modular approach and its sensitivity to key assumptions are evaluated.

Scattered-light imaging in vivo tracks fast and slow processes of neurophysiological activation.

We imaged fast optical changes associated with evoked neural activation in the dorsal brainstem of anesthetized rats, using a novel imaging device. The imager consisted of a gradient-index (GRIN) lens, a microscope objective, and a miniature charged-coupled device (CCD) video camera. We placed the probe in contact with tissue above cardiorespiratory areas of the nucleus of the solitary tract and illuminated the tissue with 780-nm light through flexible fibers around the probe perimeter. The focus depth was adjusted by moving the camera and microscope objective relative to the fixed GRIN lens. Back-scattered light images were relayed through the GRIN lens to the CCD camera. Video frames were digitized at 100 frames per second, along with tracheal pressure, arterial blood pressure, and electrocardiogram signals recorded at 1 kHz per channel. A macroelectrode placed under the GRIN lens recorded field potentials from the imaged area. Aortic, vagal, and superior laryngeal nerves were dissected free of surrounding tissue within the neck. Separate shocks to each dissected nerve elicited evoked electrical responses and caused localized optical activity patterns. The optical response was modeled by four distinct temporal components corresponding to putative physical mechanisms underlying scattered light changes. Region-of-interest analysis revealed image areas which were dominated by one or more of the different time-course components, some of which were also optimally recorded at different tissue depths. Two slow optical components appear to correspond to hemodynamic responses to metabolic demand associated with activation. Two fast optical components paralleled electrical evoked responses.

Analysis and neuronal modeling of the nonlinear characteristics of a local cardiac reflex in the rat.

Previous experimental results have suggested the existence of a local cardiac reflex in the rat. In this study, the putative role of such a local reflex in cardiovascular regulation is quantitatively analyzed. A model for the local reflex is developed from anatomical experimental results and physiological data in the literature. Using this model, a systems-level analysis is conducted. Simulation results indicate that the neuromodulatory mechanism of the local reflex attenuates the nonlinearity of the relationship between cardiac vagal drive and arterial pressure. This behavior is characterized through coherence analysis. Furthermore, the modulation of phase-related characteristics of the cardiovascular system is suggested as a plausible mechanism for the nonlinear attenuation. Based on these results, it is plausible that the functional role of the local reflex is highly robust nonlinear compensation at the heart, which results in less complex dynamics in other parts of the reflex.

Response properties of baroreceptive NTS neurons.

Neurons in the nucleus of the solitary tract (NTS) responding to activation of arterial baroreceptors were recorded intracellularly using patch pipettes in an in situ arterially perfused working heart-brain stem preparation of rat. Seven of 15 (i.e., 46%) of NTS neurons showed adaptive (nonlinear) excitatory synaptic response patterns during baroreceptor stimulation followed by an "evoked hyperpolarization." This evoked hyperpolarization was stimulus intensity dependent and capable of shunting out a subsequent baroreceptor input. We suggest that this adaptive response behavior may be mediated, in part, by calcium-dependent potassium currents (IKCa) since neurons showed spike frequency adaptation during step depolarizations and an after-hyperpolarization after repetitive firing. Furthermore, in in vivo anesthetized rats, NTS microinjections of either charybdotoxin (225 fmol) or apamin (4.5 pmol) to block IKCa increased the baroreceptor reflex gain. Our data purport that the responsiveness of baroreceptive NTS neurons can be regulated by intrinsic membrane conductances such as IKCa. Modulation of such conductances during either physiological (exercise) or pathophysiological (essential hypertension) conditions may lead to changes in both the operating point and gain of the baroreceptor reflex.

Information theoretic analysis of pulmonary stretch receptor spike trains.

Primary afferent neurons transduce physical, continuous stimuli into discrete spike trains. Investigators have long been interested in interpreting the meaning of the number or pattern of action potentials in attempts to decode the spike train back into stimulus parameters. Pulmonary stretch receptors (PSRs) are visceral mechanoreceptors that respond to deformation of the lungs and pulmonary tree. They provide the brain stem with feedback that is used by cardiorespiratory control circuits. In anesthetized, paralyzed, artificially ventilated rabbits, we recorded the action potential trains of individual PSRs while continuously manipulating ventilator rate and volume. We describe an information theoretic-based analytical method for evaluating continuous stimulus and spike train data that is of general applicability to any continuous, dynamic system. After adjusting spike times for conduction velocity, we used a sliding window to discretize the stimulus (average tracheal pressure) and response (number of spikes), and constructed co-occurrence matrices. We systematically varied the number of categories into which the stimulus and response were evenly divided at 26 different sliding window widths (5, 10, 20, 30,..., 230, 240, 250 ms). Using the probability distributions defined by the co-occurrence matrices, we estimated associated stimulus, response, joint, and conditional entropies, from which we calculated information transmitted as a fraction of the maximum possible, as well as encoding and decoding efficiencies. We found that, in general, information increases rapidly as the sliding window width increases from 5 to approximately 50 ms and then saturates as observation time increases. In addition, the information measures suggest that individual PSRs transmit more "when" than "what" type of information about the stimulus, based on the finding that the maximum information at a given window width was obtained when the stimulus was divided into just a few (usually <6) categories. Our results indicate that PSRs provide quite reliable information about tracheal pressure, with each PSR conveying about 31% of the maximum possible information about the dynamic stimulus, given our analytical parameters. When the stimulus and response are divided into more categories, slightly less information is transmitted, and this quantity also saturates as a function of observation time. We consider and discuss the importance of information contained in window widths on the time scales of an excitatory postsynaptic potential and Hering-Breuer reflex central delay.

Computational modeling of the baroreflex arc: nucleus tractus solitarius.

In this study, we examine the utility of computational modeling in understanding nervous system function. We start by examining the reasons for, and major approaches to, computational modeling. We then chose a modeling approach and applied different variations to understanding nucleus tractus solitarius (NTS) neuronal responses to various baroreceptive stimuli. We examine the results in light of our objectives and with regard to the known parameters of the system under investigation. Our results demonstrate that modeling can be a useful tool in analysis of (and examination of underlying mechanisms for) NTS behavior on many levels.

Projections of the dorsal motor nucleus of the vagus to cardiac ganglia of rat atria: an anterograde tracing study.

We injected the anterograde fluorescent tracer 1,1'-dioleyl-3,3,3',3'-tetramethylindocarbocyanine methanesulfonate (DiI) into the dorsal motor nucleus of the vagus (DmnX), counterstained the cardiac ganglia with Fluorogold (FG), and used confocal microscopy to examine the distributions and different types of DmnX fibers in wholemounts of the atria. We also quantified the number of DmnX cardiac axons and the number of innervated cardiac principal neurons (PNs). Rats with unilateral DiI injections were used in three different experiments, including unilateral FG soaking of cervical vagal trunks, intracranially rhizotomizing the vagal afferent roots, or contralaterally sectioning the cervical vagus. These manipulations indicated that DiI-labeled cardiac fibers were exclusively from the DmnX. Our observations established that: (1) three major ganglionic plexuses were localized in the epicardium; (2) both sides of the DmnX supplied significant fibers to each of the plexuses; (3) these cardiac efferents formed dense basket terminals around individual PNs; (4) collaterals of individual DmnX fibers diverged, producing calyx endings on multiple PNs; (5) small intensely fluorescent (SIF) cells in the cardiac plexuses were innervated pericellularly; (6) individual axons could innervate both PNs and SIF cells; (7) the total number of DmnX fibers were in the range of [68, 96; left] and [67, 115; right]; (8) these fibers innervated 709 (left) and 494 (right), or at least 18% and 12%, of the PNs, respectively; and (9) vagal preganglionics exhibited a degree of lateralization: Significantly more PNs were contacted by fiber varicosities in the sinoatrial plexus than in the atrioventricular plexus after right DmnX injections. In summary, the present observations suggest that the DmnX plays a significant role(s) in controlling the heart.

A neuro-mimetic dynamic scheduling algorithm for control: analysis and applications.

A simple neuronal network model of the baroreceptor reflex is analyzed. From a control perspective, the analysis suggests a dynamic scheduled control mechanisms by which the baroreflex may perform regulation of the blood pressure. The main objectives of this work are to investigate the static and dynamic response characteristics of the single neurons and the network, to analyze the neuromimetic dynamic scheduled control function of the model, and to apply the algorithm to nonlinear process control problems. The dynamic scheduling activity of the network is exploited in two control architectures. Control structure I is drawn directly from the present model of the baroreceptor reflex. An application of this structure for level control in a conical tank is described. Control structure II employs an explicit set point to determine the feedback error. The performance of this control structure is illustrated on a nonlinear continuous stirred tank reactor with van de Vusse kinetics. The two case studies validate the dynamic scheduled control approach for nonlinear process control applications.

Vagal afferent innervation of the atria of the rat heart reconstructed with confocal microscopy.

We have used confocal microscopy to analyze the vagal afferent innervation of the rat heart. Afferents were labeled by injecting 1,1'-dioleyl-3,3,3',3'-tetramethylindocarbocyanine methanesulfonate (DiI) into the nodose ganglia of animals with prior supranodose de-efferentations, autonomic ganglia were stained with Fluoro-gold, and tissues were examined in whole mounts. Distinctively different fiber specializations were observed in the epi-, myo-, and endocardium: Afferents to the epicardium formed complexes associated with cardiac ganglia. These ganglia consisted of four major ganglionated plexuses, two on each atrium, at junctions of the major vessels with the atria. Ganglionic locations and sizes (left > right) were consistent across animals. In addition to principal neurons (PNs), significant numbers of small intensely fluorescent (SIF) cells were located in each of these plexuses, and vagal afferents provided dense pericellular varicose endings around the SIF cells in each ganglionic plexus, with few if any terminations on PNs. In the myocardium, vagal afferents formed close contacts with cardiac muscles, including conduction fibers. In the endocardium, vagal fibers formed "flower-spray" and "end-net" terminals in connective tissue. With three-dimensional reconstruction of confocal optical sections, a novel polymorphism was seen: Some fibers had one or more collaterals ending as endocardial flower sprays and other collaterals ending as myocardial intramuscular endings. Some unipolar or pseudounipolar neurons within each cardiac ganglionic plexus were retrogradely labeled from the nodose ganglia. In conclusion, vagal afferents form a heterogeneity of differentiated endings in the heart, including structured elements which may mediate chemoreceptor function, stretch reception, and local cardiac reflexes.

Modeling neural mechanisms for genesis of respiratory rhythm and pattern. II. Network models of the central respiratory pattern generator.

The present paper describes several models of the central respiratory pattern generator (CRPG) developed employing experimental data and current hypotheses for respiratory rhythmogenesis. Each CRPG model includes a network of respiratory neuron types (e.g., early inspiratory; ramp inspiratory; late inspiratory; decrementing expiratory; postinspiratory; stage II expiratory; stage II constant firing expiratory; preinspiratory) and simplified models of lung and pulmonary stretch receptors (PSR), which provide feedback to the respiratory network. The used models of single respiratory neurons were developed in the Hodgkin-Huxley style as described in the previous paper. The mechanism for termination of inspiration (the inspiratory off-switch) in all models operates via late-I neuron, which is considered to be the inspiratory off-switching neuron. Several two- and three-phase CRPG models have been developed using different accepted hypotheses of the mechanism for termination of expiration. The key elements in the two-phase models are the early-I and dec-E neurons. The expiratory off-switch mechanism in these models is based on the mutual inhibitory connections between early-I and dec-E and adaptive properties of the dec-E neuron. The difference between the two-phase models concerns the mechanism for ramp firing patterns of E2 neurons resulting either from the intrinsic neuronal properties of the E2 neuron or from disinhibition from the adapting dec-E neuron. The key element of the three-phase models is the pre-I neuron, which acts as the expiratory off-switching neuron. The three-phase models differ by the mechanisms used for termination of expiration and for the ramp firing patterns of E2 neurons. Additional CRPG models were developed employing a dual switching neuron that generates two bursts per respiratory cycle to terminate both inspiration and expiration. Although distinctly different each model generates a stable respiratory rhythm and shows physiologically plausible firing patterns of respiratory neurons with and without PSR feedback. Using our models, we analyze the roles of different respiratory neuron types and their interconnections for the respiratory rhythm and pattern generation. We also investigate the possible roles of intrinsic biophysical properties of different respiratory neurons in controlling the duration of respiratory phases and timing of switching between them. We show that intrinsic membrane properties of respiratory neurons are integrated with network properties of the CRPG at three hierarchical levels: at the cellular level to provide the specific firing patterns of respiratory neurons (e.g., ramp firing patterns); at the network level to provide switching between the respiratory phases; and at the systems level to control the duration of inspiration and expiration under different conditions (e.g., lack of PSR feedback).

Modeling neural mechanisms for genesis of respiratory rhythm and pattern. I. Models of respiratory neurons.

The general objectives of our research, presented in this series of papers, were to develop a computational model of the brain stem respiratory neural network and to explore possible neural mechanisms that provide the genesis of respiratory oscillations and the specific firing patterns of respiratory neurons. The present paper describes models of single respiratory neurons that have been used as the elements in our network models of the central respiratory pattern generator presented in subsequent papers. The models of respiratory neurons were developed in the Hodgkin-Huxley style employing both physiological and biophysical data obtained from brain stem neurons in mammals. Two single respiratory neuron models were developed to match the two distinct firing behaviors of respiratory neurons described in vivo: neuron type I shows an adapting firing pattern in response to synaptic excitation, and neuron type II shows a ramp firing pattern during membrane depolarization after a period of synaptic inhibition. We found that a frequency ramp firing pattern can result from intrinsic membrane properties, specifically from the combined influence of calcium-dependent K(AHP)(Ca), low-threshold Ca(T) and K(A) channels. The neuron models with these ionic channels (type II) demonstrated ramp firing patterns similar to those recorded from respiratory neurons in vivo. Our simulations show that K(AHP)(Ca) channels in combination with high-threshold Ca(L) channels produce spike frequency adaptation during synaptic excitation. However, in combination with low-threshold Ca(T) channels, they cause a frequency ramp firing response after release from inhibition. This promotes a testable hypothesis that the main difference between the respiratory neurons that adapt (for example, early inspiratory, postinspiratory, and decrementing expiratory) and those that show ramp firing patterns (for example, ramp inspiratory and augmenting expiratory) consists of a ratio between the two types of calcium channels: Ca(L) channels predominate in the former and Ca(T) channels in the latter respiratory neuron types. We have analyzed the dependence of adapting and ramp firing patterns on maximal conductances of different ionic channels and values of synaptic drive. The effect of adjusting specific membrane conductances and synaptic interactions revealed plausible neuronal mechanisms that may underlie modulatory effects on respiratory neuron firing patterns and network performances. The results of computer simulation provide useful insight into functional significance of specific intrinsic membrane properties and their interactions with phasic synaptic inputs for a better understanding of respiratory neuron firing behavior.

Modeling neural mechanisms for genesis of respiratory rhythm and pattern. III. Comparison of model performances during afferent nerve stimulation.

The goal of the present study was to evaluate the relative plausibility of the models of the central respiratory pattern generator (CRPG) proposed in our previous paper. To test the models, we compared changes in generated patterns with the experimentally observed alterations of the respiratory pattern induced by various stimuli applied to superior laryngeal (SLN), vagus and carotid sinus (CS) nerves. In all models, short-duration SLN simulation caused phase-resetting behavior consistent with experimental data. Relatively weak sustained SLN stimulation elicited a two-phase rhythm comprising inspiration and postinspiration whereas a stronger stimulation stopped oscillations in the postinspiratory phase ("postinspiratory apnea"). In all models, sustained vagus nerve stimulation produced postinspiratory apnea. A short vagal stimulus delivered during inspiration terminated this phase. The threshold for inspiratory termination decreased during the course of the inspiratory phase. The effects of short-duration vagal stimulation applied during expiration were different in different models. In model 1, stimuli delivered in the postinspiratory phase prolonged expiration whereas the late expiratory phase was insensitive to vagal stimulation. No insensitive period was found in model 2 because vagal stimuli delivered at any time during expiration prolonged this phase. Model 3 demonstrated a short period insensitive to vagal stimulation at the very end of expiration. When phasic CS nerve stimulation was applied during inspiration or the first half of expiration, the performances of all models were similar and consistent with experimental data: stimuli delivered at the beginning inspiration shortened this phase whereas stimuli applied in the middle or at the end of inspiration prolonged it and stimuli delivered in the first half of expiration prolonged the expiratory interval. Behavior of the models were different when CS stimuli were delivered during the late expiratory phase. In model 1, these stimuli were ineffective or shortened expiration initiating the next inspiration. Alternatively, in models 2 and 3, they caused a prolongation of expiration. Although all CRPG models demonstrated a number of plausible alterations in the respiratory pattern elicited by afferent nerve stimulation, the behavior of model 1 was most consistent with experimental data. Taking into account differences in the model architectures and employed neural mechanisms, we suggest that the concept of respiratory rhythmogenesis based on the essential role of postinspiratory neurons is more plausible than the concept employing specific functional properties of decrementing expiratory (dec-E) neurons and that the ramp firing pattern of the late expiratory neuron is more likely to reflect intrinsic properties than disinhibition from the dec-E neurons.

A laser confocal microscopic study of vagal afferent innervation of rat aortic arch: chemoreceptors as well as baroreceptors.

Although the aortic nerves contain vagal afferents that terminate in both the wall of the aortic arch (putative baroreceptors) and its associated glomus tissue (putative chemoreceptors) in most mammalian species, the aortic nerves of the rat have been widely assumed to contain only baro- or pressor afferents. The present study reconsidered this anomaly by characterizing vagal afferent endings and their targets in the aortic arch region of the rat, both qualitatively and quantitatively. Eight Sprague-Dawley rats received intracranial vagal motor rhizotomy unilaterally to eliminate efferents in the nerve and then, two weeks later, injections of the tracer DiI (1,1'-dioleyl-3,3,3',3'-tetramethylindocarbocyanine methanesulfonate) into the ipsilateral nodose ganglion. The aortic arch and its surrounding tissue, with the common carotid and subclavian arteries attached, were examined with both conventional epifluorescence and confocal microscopes. Consistent with earlier observations, vagal afferents formed both flower-spray and end-net terminals rather diffusely within the wall of the aortic arch. More interestingly, vagal afferents also innervated glomus or SIF (i.e., small intensely fluorescent) cell bodies at the junction areas of the common carotid and subclavian arteries. To identify the course of these fibers, six additional animals received DiI injection into the nodose unilaterally after a complete cervical vagotomy caudal to the nodose; in these animals, the aortic nerve had been separated from the vagal trunk and kept intact. There were no marked differences in innervation patterns between the nonvagotomized and the cervically vagotomized animals, indicating that the vagal axons innervating the walls of the blood vessels and the SIF cells in the aortic arch region travel through the aortic nerves. Using a stereological method, we estimated the relative number of chemo- and baroreceptor afferents innervating the aortic arch. About 16.4% (left) and 13.1% (right) of fibers in the aortic nerves innervate SIF cells. These findings challenge the general consensus that the aortic nerves of rats contain exclusively baroreceptor fibers.

Simultaneous encoding of carotid sinus pressure and dP/dt by NTS target neurons of myelinated baroreceptors.

1. We seek to understand the baroreceptor signal processing that occurs centrally, beginning with the transformation of the signal at the first stage of processing. Because quantitative descriptions of the encoding of mean arterial pressure and its derivative with respect to time by baroreceptive second-order neurons have been unavailable, we characterized the responses of nucleus tractus solitarius (NTS) neurons that receive direct myelinated baroreceptor inputs to combinations of these two stimulus variables. 2. In anesthetized, paralyzed, artificially ventilated rabbits, the carotid sinus was vascularly isolated and the carotid sinus nerve was dissected free from surrounding tissue. Single-unit extracellular recordings were made from NTS neurons that received direct (with the use of physiological criteria) synaptic inputs from carotid sinus baroreceptors with myelinated axons. The vast majority of these neurons did not receive ipsilateral aortic nerve convergent inputs. With the use of a computer-controlled linear motor, a piecewise linear pressure waveform containing 32 combinations of pressure and its rate of change with respect to time (dP/dt) was delivered to the ipsilateral carotid sinus. 3. The average NTS firing frequency during the different stimulus combinations of pressure and dP/dt was a nonlinear and interdependent function of both variables. Most notable was the "extinctive" encoding of carotid sinus pressure by these neurons. This was characterized by an increase in firing frequency going from low to medium mean pressures (analyzed at certain positive dP/dt values) followed by a decrease in activity during high-pressure stimuli. All second-order neurons analyzed had their maximal firing rates when dP/dt was positive. 4. All neurons had their maximal firing frequency locations ("receptive field centers") at just 3 of 32 possible pressure-dP/dt coordinates. The responses of a small population of neurons were used to generate a composite description of the encoding of pressure and dP/dt. When combined as a composite of individually normalized values, the encoding of carotid sinus pressure and dP/dt may be approximated with the use of two-dimensional Gaussian functions. 5. We conclude that the population of NTS neurons recorded most faithfully encodes the rate and direction of (mean) pressure change, as opposed to providing the CNS with an unambiguous encoding of absolute pressure. Instead, the activity of these neurons, individually or as a population, serves as an estimate for the first derivative of the myelinated baroreceptor signal's encoding of mean pressure. We therefore speculate that the output of these individual neurons is useful in dynamic, rather than static, arterial pressure control.

Analysis of heart rate-based control of arterial blood pressure.

Vagal control of the heart is the most rapidly responding limb of the arterial baroreflex. We created a mathematical model of the left heart and vascular system to evaluate the ability of heart rate to influence blood pressure. The results show that arterial pressure depends nonlinearly on rate and that changes in rate are of limited effectiveness, particularly when rate is increased above the basal level. A 10% change in heart rate from rest causes a change of only 2.4% in arterial pressure due to the reciprocal relation between heart rate and stroke volume; at higher rates, insufficient filling time causes stroke volume to fall. These findings agree well with published experimental data and challenge the idea that changes in heart rate alone can strongly and rapidly affect arterial pressure. Possible implications are that vagally mediated alterations in inotropic and dromotropic state, which are not included in this model, play important roles in the fast reflex control of blood pressure or that the vagal limb of the baroreflex is of rather limited effectiveness.

Dendritic morphology of cardiac related medullary neurons defined by circuit-specific infection by a recombinant pseudorabies virus expressing beta-galactosidase.

The transneuronal herpesvirus tracer, pseudorabies virus (PRV) was used to determine the dendritic architecture of cardiac-related neurons. We constructed a derivative of the Bartha strain of PRV called PRV-BaBlu, that carries the lacZ gene of E. coli. Expression of beta-galactosidase by this recombinant virus enabled us to define the dendritic morphology of motoneurons and interneurons that innervate the heart. beta-galactosidase antigen filled dendritic processes that were clearly revealed by antibodies to beta-galactosidase. In contrast, the standard enzymatic reaction for detection of beta-galactosidase activity stained the cell soma well, but was inferior for labeling dendrites. Following PRV-BaBlu cardiac injection, infected neurons were clearly defined and labeled dendrites could be traced for long distances, sometimes greater than 800 microns from the cell body. Labeled dendrites of cardiomotor neurons primarily located in the nucleus ambiguus (NA) were extensive and sometimes intertwined with dendrites from other labeled motoneurons. Dendrites of labeled neurons in the dorsal motor nucleus of the vagus (DMV) typically extended in the mediolateral direction in the transverse plane. Transynaptically labeled interneurons interposed between the cardiorespiratory region of the nucleus tractus solitarius (NTS) and the NA were primarily located in the NA region and the reticular arc, the area between the DMV and NA. These interneurons had long dendrites extending along the reticular arc in the transverse plane. The dendritic arborizations of infected cardiac-related neurons in the NTS were variable in extent. We conclude that antibody detection of beta-galactosidase expressed by PRV-BaBlu after infection of neural cardiac circuits provides a superior method to define the dendrites and dendritic fields of cardiac-related motoneurons and interneurons.

Central neuronal circuit innervating the rat heart defined by transneuronal transport of pseudorabies virus.

We defined the central circuit innervating various regions of the rat heart using a neurotropic herpesvirus as a transneuronal tracer. Location of viral antigens in the brain after cardiac injection of three strains of pseudorabies virus (PRV) provided insight into vagal preganglionic neurons and their connected interneurons. At short survival times, labeled vagal preganglionic neurons were localized in both the nucleus ambiguus (NA) and the dorsal motor nucleus of the vagus (DMV), and in an arc-like band through the reticular formation between the NA and the DMV. The amount of DMV labeling was dependent on viral strain. Similar distributions of labeled neurons were observed following either ganglionic, sinoatrial node, or ventricular injections. At intermediate survival times postcardiac injection, the virus replicated in vagal preganglionic neurons and was trans-synaptically transported to interneurons observed primarily in the NA regions and in an arc-like band through the reticular formation. Labeled neurons were also observed in ventral regions of the nucleus of the solitary tract (NTS). At longer survival times, labeled neurons were found in various regions of the NTS with the most abundant label dorsal and dorsomedial to the solitary tract. Abundant neuronal labeling was also found in the intermediolateral cell column, the raphe nuclei, the caudal and rostral ventral lateral medulla, the A5 region, the locus coeruleus, and the lateral and paraventricular hypothalamic nuclei. These data define the central circuits including the interneuronal connections that innervate various cardiac targets.

A digital brain atlas and its application to the visceral neuraxis.

We describe the details and application of a digital brain atlas for the comparison and integration of graphical neurobiological data. The atlas consists of multiple sets of high-resolution video images acquired from histological tissue sections representing a 3-dimensional (3D) volume of an exemplar rat brain. Through an interactive graphical interface running on a standard computer workstation, experimental data is brought into register with the atlas. Once in the atlas, coordinate reference frame data can be compared, analyzed, and visualized in 3 dimensions. We demonstrate the validity and usefulness of the digital brain atlas with a series of results on the visceral neuraxis in the rat.

Innervation of the heart and its central medullary origin defined by viral tracing.

The vagus nerve exerts a profound influence on the heart, regulating the heart rate and rhythm. An extensive vagal innervation of the cardiac ventricles and the central origin and extent of this innervation was demonstrated by transynaptic transport of pseudorabies virus with a virulent and two attenuated pseudorabies viral strains. The neurons that innervate the ventricles are numerous, and their distribution within the nucleus ambiguus and dorsal motor nucleus of the vagus is similar to that of neurons innervating other cardiac targets, such as the sino-atrial node. These data provide a neuroanatomical correlate to the physiological influence of the vagus nerve on ventricular function.