Complex spatiotemporal haemodynamic response following sensory stimulation in the awake rat

Detailed understanding of the haemodynamic changes that underlie non-invasive neuroimaging techniques such as blood oxygen level dependent functional magnetic resonance imaging is essential if we are to continue to extend the use of these methods for understanding brain function and dysfunction. The use of animal and in particular rodent research models has been central to these endeavours as they allow in-vivo experimental techniques that provide measurements of the haemodynamic response function at high temporal and spatial resolution. A limitation of most of this research is the use of anaesthetic agents which may disrupt or mask important features of neurovascular coupling or the haemodynamic response function. In this study we therefore measured spatiotemporal cortical haemodynamic responses to somatosensory stimulation in awake rats using optical imaging spectroscopy. Trained, restrained animals received non-noxious stimulation of the whisker pad via chronically implanted stimulating microwires whilst optical recordings were made from the contralateral somatosensory cortex through a thin cranial window. The responses we measure from un-anaesthetised animals are substantially different from those reported in previous studies which have used anaesthetised animals. These differences include biphasic response regions (initial increases in blood volume and oxygenation followed by subsequent decreases) as well as oscillations in the response time series of awake animals. These haemodynamic response features do not reflect concomitant changes in the underlying neuronal activity and therefore reflect neurovascular or cerebrovascular processes. These hitherto unreported hyperemic response dynamics may have important implications for the use of anaesthetised animal models for research into the haemodynamic response function.

Balanced excitation and inhibition: Model based analysis of local field potentials

We have developed a model of the local field potential (LFP) based on the conservation of charge, the independence principle of ionic flows and the classical Hodgkin–Huxley (HH) type intracellular model of synaptic activity. Insights were gained through the simulation of the HH intracellular model on the nonlinear relationship between the balance of synaptic conductances and that of post-synaptic currents. The latter is dependent not only on the former, but also on the temporal lag between the excitatory and inhibitory conductances, as well as the strength of the afferent signal. The proposed LFP model provides a method for decomposing the LFP recordings near the soma of layer IV pyramidal neurons in the barrel cortex of anaesthetised rats into two highly correlated components with opposite polarity. The temporal dynamics and the proportional balance of the two components are comparable to the excitatory and inhibitory post-synaptic currents computed from the HH model. This suggests that the two components of the LFP reflect the underlying excitatory and inhibitory post-synaptic currents of the local neural population. We further used the model to decompose a sequence of evoked LFP responses under repetitive electrical stimulation (5 Hz) of the whisker pad. We found that as neural responses adapted, the excitatory and inhibitory components also adapted proportionately, while the temporal lag between the onsets of the two components increased during frequency adaptation. Our results demonstrated that the balance between neural excitation and inhibition can be investigated using extracellular recordings. Extension of the model to incorporate multiple compartments should allow more quantitative interpretations of surface Electroencephalography (EEG) recordings into components reflecting the excitatory, inhibitory and passive ionic current flows generated by local neural populations.

Does neural input or processing play a greater role in the magnitude of neuroimaging signals?

An important constraint on how hemodynamic neuroimaging signals such as fMRI can be interpreted in terms of the underlying evoked activity is an understanding of neurovascular coupling mechanisms that actually generate hemodynamic responses. The predominant view at present is that the hemodynamic response is most correlated with synaptic input and subsequent neural processing rather than spiking output. It is still not clear whether input or processing is more important in the generation of hemodynamics responses. In order to investigate this we measured the hemodynamic and neural responses to electrical whisker pad stimuli in rat whisker barrel somatosensory cortex both before and after the local cortical injections of the GABAA agonist muscimol. Muscimol would not be expected to affect the thalamocortical input into the cortex but would inhibit subsequent intra-cortical processing. Pre-muscimol infusion whisker stimuli elicited the expected neural and accompanying hemodynamic responses to that reported previously. Following infusion of muscimol, although the temporal profile of neural responses to each pulse of the stimulus train was similar, the average response was reduced in magnitude by ∼79% compared to that elicited pre-infusion. The whisker-evoked hemodynamic responses were reduced by a commensurate magnitude suggesting that, although the neurovascular coupling relationships were similar for synaptic input as well as for cortical processing, the magnitude of the overall response is dominated by processing rather than from that produced from the thalamocortical input alone.

Long duration stimuli and nonlinearities in the neural–haemodynamic coupling

Recent studies have shown that the haemodynamic responses to brief (<2 secs) stimuli can be well characterised as a linear convolution of neural activity with a suitable haemodynamic impulse response. In this paper, we show that the linear convolution model cannot predict measurements of blood flow responses to stimuli of longer duration (>2 secs), regardless of the impulse response function chosen. Modifying the linear convolution scheme to a nonlinear convolution scheme was found to provide a good prediction of the observed data. Whereas several studies have found a nonlinear coupling between stimulus input and blood flow responses, the current modelling scheme uses neural activity as an input, and thus implies nonlinearity in the coupling between neural activity and blood flow responses. Neural activity was assessed by current source density analysis of depth-resolved evoked field potentials, while blood flow responses were measured using laser Doppler flowmetry. All measurements were made in rat whisker barrel cortex after electrical stimulation of the whisker pad for 1 to 16 secs at 5 Hz and 1.2 mA (individual pulse width 0.3 ms).

The hemodynamic impulse response to a single neural event

This article investigates the relation between stimulus-evoked neural activity and cerebral hemodynamics. Specifically, the hypothesis is tested that hemodynamic responses can be modeled as a linear convolution of experimentally obtained measures of neural activity with a suitable hemodynamic impulse response function. To obtain a range of neural and hemodynamic responses, rat whisker pad was stimulated using brief (less than or equal to2 seconds) electrical stimuli consisting of single pulses (0.3 millisecond, 1.2 mA) combined both at different frequencies and in a paired-pulse design. Hemodynamic responses were measured using concurrent optical imaging spectroscopy and laser Doppler flowmetry, whereas neural responses were assessed through current source density analysis of multielectrode recordings from a single barrel. General linear modeling was used to deconvolve the hemodynamic impulse response to a single "neural event" from the hemodynamic and neural responses to stimulation. The model provided an excellent fit to the empirical data. The implications of these results for modeling schemes and for physiologic systems coupling neural and hemodynamic activity are discussed.

Long-latency reductions in gamma power predict hemodynamic changes that underlie the negative BOLD signal

Studiesthat use prolonged periods of sensory stimulation report associations between regional reductions in neural activity and negative blood oxygenation level-dependent (BOLD) signaling. However, the neural generators of the negative BOLD response remain to be characterized. Here, we use single-impulse electrical stimulation of the whisker pad in the anesthetized rat to identify components of the neural response that are related to “negative” hemodynamic changes in the brain. Laminar multiunit activity and local field potential recordings of neural activity were performed concurrently withtwo-dimensional optical imaging spectroscopy measuring hemodynamic changes. Repeated measurements over multiple stimulation trials revealed significant variations in neural responses across session and animal datasets. Within this variation, we found robust long-latency decreases (300 and 2000 ms after stimulus presentation) in gammaband power (30 – 80 Hz) in the middle-superficial cortical layers in regions surrounding the activated whisker barrel cortex. This reduction in gamma frequency activity was associated with corresponding decreases in the hemodynamic responses that drive the negative BOLD signal. These findings suggest a close relationship between BOLD responses and neural events that operate over time scales that outlast the initiating sensory stimulus, and provide important insights into the neurophysiological basis of negative neuroimaging signals.

Neural network activity in the neonatal rat barrel cortex in vivo

Coordinated patterns of electrical activity are important for the early development of sensory systems. The spatiotemporal dynamics of these early activity patterns and the role of the peripheral sensory input for their generation are essentially unknown. There are two projects in this thesis. In project1, we performed extracellular multielectrode recordings in the somatosensory cortex of postnatal day 0 to 7 rats in vivo and observed three distinct patterns of synchronized oscillatory activity. (1) Spontaneous and periphery-driven spindle bursts of 1–2 s in duration and ~10 Hz in frequency occurred approximately every 10 s. (2) Spontaneous and sensory-driven gamma oscillations of 150–300 ms duration and 30–40 Hz in frequency occurred every 10–30 s. (3) Long oscillations appeared only every ~20 min and revealed the largest amplitude (250–750 µV) and longest duration (>40 s). These three distinct patterns of early oscillatory activity differently synchronized the neonatal cortical network. Whereas spindle bursts and gamma oscillations did not propagate and synchronized a local neuronal network of 200–400 µm in diameter, long oscillations propagated with 25–30 µm/s and synchronized 600-800 µm large ensembles. All three activity patterns were triggered by sensory activation. Single electrical stimulation of the whisker pad or tactile whisker activation elicited neocortical spindle bursts and gamma activity. Long oscillations could be only evoked by repetitive sensory stimulation. The neonatal oscillatory patterns in vivo depended on NMDAreceptor-mediated synaptic transmission and gap junctional coupling. Whereas spindle bursts and gamma oscillations may represent an early functional columnar-like pattern, long oscillations may serve as a propagating activation signal consolidating these immature neuronal networks. In project2, Using voltage-sensitive dye imaging and simultaneous multi-channel extracellular recordings in the barrel cortex and somatosensory thalamus of newborn rats in vivo, we found that spontaneous and whisker stimulation induced activity patterns were restricted to functional cortical columns already at the day of birth. Spontaneous and stimulus evoked cortical activity consisted of gamma oscillations followed by spindle bursts. Spontaneous events were mainly generated in the thalamus or by spontaneous whisker movements. Our findings indicate that during early developmental stages cortical networks self-organize in ontogenetic columns via spontaneous gamma oscillations triggered by the thalamus or sensory periphery.

TRPV4 is necessary for trigeminal irritant pain and functions as a cellular formalin receptor.

Detection of external irritants by head nociceptor neurons has deep evolutionary roots. Irritant-induced aversive behavior is a popular pain model in laboratory animals. It is used widely in the formalin model, where formaldehyde is injected into the rodent paw, eliciting quantifiable nocifensive behavior that has a direct, tissue-injury-evoked phase, and a subsequent tonic phase caused by neural maladaptation. The formalin model has elucidated many antipain compounds and pain-modulating signaling pathways. We have adopted this model to trigeminally innervated territories in mice. In addition, we examined the involvement of TRPV4 channels in formalin-evoked trigeminal pain behavior because TRPV4 is abundantly expressed in trigeminal ganglion (TG) sensory neurons, and because we have recently defined TRPV4's role in response to airborne irritants and in a model for temporomandibular joint pain. We found TRPV4 to be important for trigeminal nocifensive behavior evoked by formalin whisker pad injections. This conclusion is supported by studies with Trpv4(-/-) mice and TRPV4-specific antagonists. Our results imply TRPV4 in MEK-ERK activation in TG sensory neurons. Furthermore, cellular studies in primary TG neurons and in heterologous TRPV4-expressing cells suggest that TRPV4 can be activated directly by formalin to gate Ca(2+). Using TRPA1-blocker and Trpa1(-/-) mice, we found that both TRP channels co-contribute to the formalin trigeminal pain response. These results imply TRPV4 as an important signaling molecule in irritation-evoked trigeminal pain. TRPV4-antagonistic therapies can therefore be envisioned as novel analgesics, possibly for specific targeting of trigeminal pain disorders, such as migraine, headaches, temporomandibular joint, facial, and dental pain, and irritation of trigeminally innervated surface epithelia.

Bulk compressive properties of the heel fat pad during walking: a pilot investigation in plantar heel pain

Background: Altered mechanical properties of the heel pad have been implicated in the development of plantar heel pain. However, the in vivo properties of the heel pad during gait remain largely unexplored in this cohort. The aim of the current study was to characterise the bulk compressive properties of the heel pad in individuals with and without plantar heel pain while walking. ---------- Methods: The sagittal thickness and axial compressive strain of the heel pad were estimated in vivo from dynamic lateral foot radiographs acquired from nine subjects with unilateral plantar heel pain and an equivalent number of matched controls, while walking at their preferred speed. Compressive stress was derived from simultaneously acquired plantar pressure data. Principal viscoelastic parameters of the heel pad, including peak strain, secant modulus and energy dissipation (hysteresis), were estimated from subsequent stress–strain curves.---------- Findings: There was no significant difference in loaded and unloaded heel pad thickness, peak stress, peak strain, or secant and tangent modulus in subjects with and without heel pain. However, the fat pad of symptomatic feet had a significantly lower energy dissipation ratio (0.55 ± 0.17 vs. 0.69 ± 0.08) when compared to asymptomatic feet (P < .05).---------- Interpretation: Plantar heel pain is characterised by reduced energy dissipation ratio of the heel pad when measured in vivo and under physiologically relevant strain rates.

Plantar enthesopathy : thickening of the enthesis is correlated with energy dissipation of the plantar fat pad during walking

Background: The enthesis of the plantar fascia is thought to play an important role in stress dissipation. However, the potential link between entheseal thickening characteristic of enthesopathy and the stress-dissipating properties of the intervening plantar fat pad have not been investigated. Purpose: This study was conducted to identify whether plantar fat pad mechanics explain variance in the thickness of the fascial enthesis in individuals with and without plantar enthesopathy. Study Design: Case-control study; Level of evidence, 3. Methods: The study population consisted of 9 patients with unilateral plantar enthesopathy and 9 asymptomatic, individually matched controls. The thickness of the enthesis of the symptomatic, asymptomatic, and a matched control limb was acquired using high-resolution ultrasound. The compressive strain of the plantar fat pad during walking was estimated from dynamic lateral radiographs acquired with a multifunction fluoroscopy unit. Peak compressive stress was simultaneously acquired via a pressure platform. Principal viscoelastic parameters were estimated from subsequent stress-strain curves. Results: The symptomatic fascial enthesis (6.7 ± 2.0 mm) was significantly thicker than the asymptomatic enthesis (4.2 ± 0.4 mm), which in turn was thicker than the enthesis (3.3 ± 0.4 mm) of control limbs (P < .05). There was no significant difference in the mean thickness, peak stress, peak strain, or secant modulus of the plantar fat pad between limbs. However, the energy dissipated by the fat pad during loading and unloading was significantly lower in the symptomatic limb (0.55 ± 0.17) when compared with asymptomatic (0.69 ± 0.13) and control (0.70 ± 0.09) limbs (P < .05). The sonographic thickness of the enthesis was correlated with the energy dissipation ratio of the plantar fat pad (r = .72, P < .05), but only in the symptomatic limb. Conclusion: The energy-dissipating properties of the plantar fat pad are associated with the sonograpic appearance of the enthesis in symptomatic limbs, providing a previously unidentified link between the mechanical behavior of the plantar fat pad and enthesopathy.

An isoviscous, isothermal model investigating the influence of hydrostatic recesses on a spring-supported tilting pad thrust bearing

Tilting-pad hydrodynamic thrust bearings are used in hydroelectric power stations around the world, reliably supporting turbines weighing hundreds of tonnes, over decades of service. Newer designs incorporate hydrostatic recesses machined into the sector-shaped pads to enhance oil film thickness at low rotational speeds. External pressurisation practically eliminates wear and enhances service life and reliability. It follows that older generating plants, lacking such assistance, stand to benefit from being retrofitted with hydrostatic lubrication systems. The design process is not trivial however. The need to increase the groove size to permit spontaneous lifting of the turbine under hydrostatic pressure, conflicts with the need to preserve performance of the original plane pad design. A haphazardly designed recess can induce a significant rise in bearing temperature concomitant with reduced mechanical efficiency and risk of thermal damage. In this work, a numerical study of a sector-shaped pad is undertaken to demonstrate how recess size and shape can affect the performance of a typical bearing.

Dynamic effects of electrical pitting in steam-turbine tilting-pad thrust-bearings

This paper merges the analysis of a case history and the simplified theoretical model related to a rather singular phenomenon that may happen in rotating machinery. Starting from the first, a small industrial steam turbine experienced a very strange behavior during megawatt load. When the unit was approaching the maximum allowed power, the temperature of the babbitt metal of the pads of the thrust bearing showed constant increase with an unrecoverable drift. Bearing inspection showed that pad trailing edge had the typical aspect of electrical pitting. This kind of damage was not reparable and bearing pads had to replaced. This problem occurred several times in sequence and was solved only by adding further ground brushes to the shaft-line. Failure analysis indicated electrodischarge machining as the root fault. A specific model, able to take into consideration the effect of electrical pitting and loading capacity decreasing as a consequence of the damage of the babbitt metal, is proposed in the paper and shows that the phenomenon causes the irretrievable failure of the thrust bearing.

Force-deformation properties of the human heel pad during barefoot walking

Introduction: The plantar heel pad is a specialized fibroadipose tissue that attenuates and, in part, dissipates the impact energy associated with heel strike. Although near maximal deformation of the heel pad has been shown during running, in vivo measurement of the deformation and structural properties of the heel pad during walking remains largely unexplored. This study employed a fluoroscope, synchronized with a pressure platform, to obtain force–deformation data for the heel pad during walking. Methods: Dynamic lateral foot radiographs were acquired from 6 male and 10 female adults (age, 45 ± 10 yrs; height, 1.66 ± 0.10 m; and weight, 80.7 ± 10.8 kg), while walking barefoot at preferred speeds. The inferior aspect of the calcaneus was digitized and the sagittal thickness and deformation of the heel pad relative to the support surface calculated. Simultaneous measurement of the peak force beneath the heel was used to estimate the principal structural properties of the heel pad. Results: Transient loading profiles associated with walking induced rapidly changing deformation rates in the heel pad and resulted in irregular load–deformation curves. The initial stiffness (32 ± 11 N.mm-1) of the heel pad was an order of magnitude lower than its final stiffness (212 ± 125 N.mm-1) and on average, only 1.0 J of energy was dissipated by the heel pad with each step during walking. Peak deformation (10.3 mm) approached that predicted for the limit of pain tolerance (10.7 mm). Conclusion: These findings suggest the heel pad operates close to its pain threshold even at speeds encountered during barefoot walking and provides insight as to why barefoot runners may adopt ‘forefoot’ strike patterns that minimize heel loading.

Telephone and security pad at the entrance of the ZFA (Zentrum für Antisemitismusforschung; Center for Research on Antisemitism); Berlin.

Digital Image

Two unidentified women looking at Ludwig Levy's sketch pad.

One may be Hanna Levy, wife of Ludwig