A dynamic causal model of the coupling between pulse stimulation and neural activity

We present a dynamic causal model that can explain context-dependent changes in neural responses, in the rat barrel cortex, to an electrical whisker stimulation at different frequencies. Neural responses were measured in terms of local field potentials. These were converted into current source density (CSD) data, and the time series of the CSD sink was extracted to provide a time series response train. The model structure consists of three layers (approximating the responses from the brain stem to the thalamus and then the barrel cortex), and the latter two layers contain nonlinearly coupled modules of linear second-order dynamic systems. The interaction of these modules forms a nonlinear regulatory system that determines the temporal structure of the neural response amplitude for the thalamic and cortical layers. The model is based on the measured population dynamics of neurons rather than the dynamics of a single neuron and was evaluated against CSD data from experiments with varying stimulation frequency (1–40 Hz), random pulse trains, and awake and anesthetized animals. The model parameters obtained by optimization for different physiological conditions (anesthetized or awake) were significantly different. Following Friston, Mechelli, Turner, and Price (2000), this work is part of a formal mathematical system currently being developed (Zheng et al., 2005) that links stimulation to the blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) signal through neural activity and hemodynamic variables. The importance of the model described here is that it can be used to invert the hemodynamic measurements of changes in blood flow to estimate the underlying neural activity.

A three-compartment model of the hemodynamic response and oxygen delivery to brain

We describe a mathematical model linking changes in cerebral blood flow, blood volume and the blood oxygenation state in response to stimulation. The model has three compartments to take into account the fact that the cerebral blood flow and volume as measured concurrently using laser Doppler flowmetry and optical imaging spectroscopy have contributions from the arterial, capillary as well as the venous compartments of the vasculature. It is an extension to previous one-compartment hemodynamic models which assume that the measured blood volume changes are from the venous compartment only. An important assumption of the model is that the tissue oxygen concentration is a time varying state variable of the system and is driven by the changes in metabolic demand resulting from changes in neural activity. The model takes into account the pre-capillary oxygen diffusion by flexibly allowing the saturation of the arterial compartment to be less than unity. Simulations are used to explore the sensitivity of the model and to optimise the parameters for experimental data. We conclude that the three-compartment model was better than the one-compartment model at capturing the hemodynamics of the response to changes in neural activation following stimulation.

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.

Lattice Boltzmann simulation of viscous fluid systems with elastic boundaries

A lattice Boltzmann model able to simulate viscous fluid systems with elastic and movable boundaries is proposed. By introducing the virtual distribution function at the boundary, the Galilean invariance is recovered for the full system. As examples of application, the how in elastic vessels is simulated with the pressure-radius relationship similar to that of the pulmonary blood vessels. The numerical results for steady how are in good agreement with the analytical prediction, while the simulation results for pulsative how agree with the experimental observation of the aortic flows qualitatively. The approach has potential application in the study of the complex fluid systems such as the suspension system as well as the arterial blood flow.

The effect of flavanol-rich cocoa on cerebral perfusion in healthy older adults during conscious resting state: a placebo controlled, crossover, acute trial

Rationale: There has recently been increasing interest in the potential of flavanols, plant derived compounds found in foods such as fruit and vegetables, to ameliorate age-related cognitive decline. Research suggests that cocoa flavanols improve memory and learning, possibly as a result of their anti-inflammatory and neuroprotective effects. These effects may be mediated by increased cerebral blood flow (CBF), thus stimulating neuronal function. Objectives: The present study employed arterial spin labelling (ASL) functional magnetic resonance imaging (FMRI) to explore the effect of a single acute dose of cocoa flavanols on regional CBF. Methods: CBF was measured pre and post consumption of low (23mg) or high (494mg) 330ml equicaloric flavanol drinks matched for caffeine, theobromine, taste and appearance according to a randomised counterbalanced crossover double-blind design in eight males and ten females, aged 50-65 years. Changes in perfusion from pre to post consumption were calculated as a function of each drink. Results: Significant increases in regional perfusion across the brain were observed following consumption of the high flavanol drink relative to the low flavanol drink, particularly in the anterior cingulate cortex (ACC) and the central opercular cortex of the parietal lobe. Conclusions: Consumption of cocoa flavanol improves regional cerebral perfusion in older adults. This provides evidence for a possible acute mechanism by which cocoa flavanols are associated with benefits for cognitive performance.

Subjective evaluation of experimental dyspnoea: effects of isocapnia and repeated exposure

Resistive respiratory loading is an established stimulus for the induction of experimental dyspnoea. In comparison to unloaded breathing, resistive loaded breathing alters end-tidal CO2 (PETCO2), which has independent physiological effects (e.g. upon cerebral blood flow). We investigated the subjective effects of resistive loaded breathing with stabilized PETCO2 (isocapnia) during manual control of inspired gases on varying baseline levels of mild hypercapnia increased PETCO2). Furthermore, to investigate whether perceptual habituation to dyspnoea stimuli occurs, the study was repeated over four experimental sessions. Isocapnic hypercapnia did not affect dyspnoea unpleasantness during resistive loading. A post hoc analysis revealed a small increase of respiratory unpleasantness during unloaded breathing at +0.6 kPa, the level that reliably induced isocapnia. We didnot observe perceptual habituation over the four sessions. We conclude that isocapnic respiratory loading allows stable induction of respiratory unpleasantness, making it a good stimulus for multi-session studies of dyspnoea.

Light, electron microscopy and histopathology of Myxobolus salminus n. sp., a parasite of Salminus brasiliensis from the Brazilian Pantanal

In this report, we describe the morphology and histopathology of Myxobolus salminus n. sp., a parasite of the gill filaments of wild Salminus brasiliensis (dourado) from the Brazilian Pantanal. The small polysporic plasmodia were similar to 100 mu m in diameter and the development was asynchronous. The mature spores were oval to pear shaped and had a smooth wall. The spore measurements were (mean +/- S.D., with range in parentheses): length 10.1 +/- 0.4 mu m (9.6-10.5), width 6.1 +/- 0.4 mu m (5.8-6.6) and thickness 5.0 +/- 0.6 mu m (4.7-5.3). The polar capsules were elongated and of equal size: length 4.6 +/- 0.2 mu m (4.3-4.8) and width 1.7 +/- 0.1 mu m (1.5-1.9). The histological analysis revealed numerous plasmodia in the blood vessels of the gill filaments. The site of parasite development was the wall of the large-caliber blood vessel of the gill filament, with progressive growth towards the lumen, resulting in the obstruction of blood flow, congestion and perivascular edema. The ultrastructural study revealed that the plasmodial wall was composed of two membranes, had numerous pinocytic canals and was in direct contact with the basement membrane of the vessel. The development of the parasite was asynchronous, with mature spores, immature spores and young developmental stages randomly distributed throughout the plasmodium. The prevalence of the parasite was 4.4%. with male and female fish being infected. (C) 2009 Elsevier B.V. All rights reserved.

Conceptus-mediated endometrial vascular changes during early pregnancy in mares: an anatomic, histomorphometric, and vascular endothelial growth factor receptor system immunolocalization and gene expression study

This work examined how the conceptus modulates endometrial tissue remodeling and vascular development prior to implantation in mares. A macroscopic uterine examination was completed at day 21 of pregnancy. In situ morphology revealed that the endometrium involved in encroachment is restricted to the dorsal endometrium immediately overlying the yolk sac. The amount of stromal area occupied by blood vessels and the number of endometrial glands were increased during early pregnancy. Endometrial histomorphometry as well as the endometrial mRNA abundance and immunolocalization of VEGF, VEGFR1, VEGFR2, and Ki-67 was completed at days 14 and 21 of pregnancy, at day 10 of the estrous cycle, and during estrus. No obvious differences in VEGF and VEGFR1 protein localization were detected between pregnant and cycling mares but differential staining pattern for VEGFR2 and Ki-67 was observed. VEGFR2 localized to luminal and glandular epithelium of pregnant mares, while luminal epithelium was negative in cycling mares. Ki-67 staining was weak during the luteal phase but exhibited prominent luminal epithelium staining during estrus. In pregnant mares, all endometrial layers were Ki-67 positive. Quantitative RT-PCR revealed a greater abundance of VEGF mRNA during pregnancy. VEGFR2 transcript abundance was greatest in pregnant mares on day 21. This study supports the concept that the conceptus plays an active role in directing vasculogenesis within the uterus and thereby establishing hemotrophic nutrition that supports pregnancy after implantation. Reproduction (2011) 142 593-603

N-acetylcysteine protects against renal injury following bilateral ureteral obstruction

Background. Obstructive nephropathy decreases renal blood flow (RBF) and glomerular filtration rate (GFR), causing tubular abnormalities, such as urinary concentrating defect, as well as increasing oxidative stress. This study aimed to evaluate the effects of N-acetylcysteine (NAC) on renal function, as well as on the protein expression of aquaporin 2 (AQP2) and endothelial nitric oxide synthase (eNOS), after the relief of bilateral ureteral obstruction (BUO). Methods. Adult male Wistar rats were divided into four groups: sham (sham operated); sham operated + 440 mg/kg body weight (BW) of NAC daily in drinking water, started 2 days before and maintained until 48 h after the surgery; BUO (24-h BUO only); BUO + NAC-pre (24-h BUO plus 440 mg/kg BW of NAC daily in drinking water started 2 days before BUO); and BUO + NAC-post (24-h BUO plus 440 mg/kg BW of NAC daily in drinking water started on the day of BUO relief). Experiments were conducted 48 h after BUO relief. Results. Serum levels of thiobarbituric reactive substances, which are markers of lipid peroxidation, were significantly lower in NAC-treated rats than in the BUO group rats. The administration of NAC provided significant protection against post-BUO GFR drops and reductions in RBF. Renal cortices and BUO rats presented decreased eNOS protein expression of eNOS in the renal cortex of BUO group rats, whereas it was partially recovered in BUO + NAC-pre group rats. Urine osmolality was significantly lower in BUO rats than in sham group rats or NAC-treated rats, the last also presenting less interstitial fibrosis. Post-BUO downregulation of AQP2 protein expression was averted in the BUO + NAC-pre group rats. Conclusions. This study demonstrates that NAC administration ameliorates the renal function impairment observed 48 h after the relief of 24-h BUO. Oxidative stress is important for the suppression of GFR, RBF, tissue AQP2 and eNOS in the polyuric phase after the release of BUO.

Effect of short-term creatine supplementation on markers of skeletal muscle damage after strenuous contractile activity

The protective effect of short-term creatine supplementation (CrS) upon markers of strenuous contractile activity-induced damage in human and rat skeletal muscles was investigated. Eight Ironman triathletes were randomized into the placebo (Pl; n = 4) and creatine-supplemented (CrS; n = 4) groups. Five days prior to the Ironman competition, the CrS group received creatine monohydrate (20 g day(-1)) plus maltodextrin (50 g) divided in two equal doses. The Pl group received maltodextrin (50 g day(-1)) only. The effect of CrS (5 g day(-1)/kg body weight for 5 days) was also evaluated in a protocol of strenuous contractile activity induced by electrical stimulation in rats. Blood samples were collected before and 36 and 60 h after the competition and were used to determine plasma activities of creatine kinase (CK), lactate dehydrogenase (LDH), aldolase (ALD), glutamic oxaloacetic acid transaminase (GOT), glutamic pyruvic acid transaminase (GPT), and C-reactive protein (CRP) level. In rats, plasma activities of CK and LDH, muscle vascular permeability (MVP) using Evans blue dye, muscle force and fatigue were evaluated. Activities of CK, ALD, LDH, GOT, GTP, and levels of CRP were increased in the Pl group after the competition as compared to basal values. CrS decreased plasma activities of CK, LDH, and ALD, and prevented the rise of GOT and GPT plasma activities. In rats, CrS delayed the fatigue, preserved the force, and prevented the rise of LDH and CK plasma activities and MVP in the gastrocnemius muscle. CrS presented a protective effect on muscle injury induced by strenuous contractile activities.

Black-box decomposition approach for computational hemodynamics: One-dimensional models

Increasing efforts exist in integrating different levels of detail in models of the cardiovascular system. For instance, one-dimensional representations are employed to model the systemic circulation. In this context, effective and black-box-type decomposition strategies for one-dimensional networks are needed, so as to: (i) employ domain decomposition strategies for large systemic models (1D-1D coupling) and (ii) provide the conceptual basis for dimensionally-heterogeneous representations (1D-3D coupling, among various possibilities). The strategy proposed in this article works for both of these two scenarios, though the several applications shown to illustrate its performance focus on the 1D-1D coupling case. A one-dimensional network is decomposed in such a way that each coupling point connects two (and not more) of the sub-networks. At each of the M connection points two unknowns are defined: the flow rate and pressure. These 2M unknowns are determined by 2M equations, since each sub-network provides one (non-linear) equation per coupling point. It is shown how to build the 2M x 2M non-linear system with arbitrary and independent choice of boundary conditions for each of the sub-networks. The idea is then to solve this non-linear system until convergence, which guarantees strong coupling of the complete network. In other words, if the non-linear solver converges at each time step, the solution coincides with what would be obtained by monolithically modeling the whole network. The decomposition thus imposes no stability restriction on the choice of the time step size. Effective iterative strategies for the non-linear system that preserve the black-box character of the decomposition are then explored. Several variants of matrix-free Broyden`s and Newton-GMRES algorithms are assessed as numerical solvers by comparing their performance on sub-critical wave propagation problems which range from academic test cases to realistic cardiovascular applications. A specific variant of Broyden`s algorithm is identified and recommended on the basis of its computer cost and reliability. (C) 2010 Elsevier B.V. All rights reserved.

Iterative strong coupling of dimensionally heterogeneous models

In this article we address decomposition strategies especially tailored to perform strong coupling of dimensionally heterogeneous models, under the hypothesis that one wants to solve each submodel separately and implement the interaction between subdomains by boundary conditions alone. The novel methodology takes full advantage of the small number of interface unknowns in this kind of problems. Existing algorithms can be viewed as variants of the `natural` staggered algorithm in which each domain transfers function values to the other, and receives fluxes (or forces), and vice versa. This natural algorithm is known as Dirichlet-to-Neumann in the Domain Decomposition literature. Essentially, we propose a framework in which this algorithm is equivalent to applying Gauss-Seidel iterations to a suitably defined (linear or nonlinear) system of equations. It is then immediate to switch to other iterative solvers such as GMRES or other Krylov-based method. which we assess through numerical experiments showing the significant gain that can be achieved. indeed. the benefit is that an extremely flexible, automatic coupling strategy can be developed, which in addition leads to iterative procedures that are parameter-free and rapidly converging. Further, in linear problems they have the finite termination property. Copyright (C) 2009 John Wiley & Sons, Ltd.

Arterial spin labeling of cerebral perfusion territories using a seperate labeling coil

Purpose: To obtain cerebral perfusion territories of the left, the right. and the posterior circulation in humans with high signal-to-noise ratio (SNR) and robust delineation. Materials and Methods: Continuous arterial spin labeling (CASL) was implemented using a dedicated radio frequency (RF) coil. positioned over the neck, to label the major cerebral feeding arteries in humans. Selective labeling was achieved by flow-driven adiabatic fast passage and by tilting the longitudinal labeling gradient about the Y-axis by theta = +/- 60 degrees. Results: Mean cerebral blood flow (CBF) values in gray matter (GM) and white matter (WM) were 74 +/- 13 mL center dot 100 g(-1) center dot minute(-1) and 14 +/- 13 mL center dot 100 g(-1) center dot minute(-1), respectively (N = 14). There were no signal differences between left and right hemispheres when theta = 0 degrees (P > 0.19), indicating efficient labeling of both hemispheres. When theta = +60 degrees, the signal in GM on the left hemisphere, 0.07 +/- 0.06%, was 92% lower than on the right hemisphere. 0.85 +/- 0.30% (P < 1 x 10(-9)). while for theta = -60 degrees, the signal in the right hemisphere. 0.16 +/- 0.13%, was 82% lower than on the contralateral side. 0.89 +/- 0.22% (P < 1 x 10(-10)). Similar attenuations were obtained in WM. Conclusion: Clear delineation of the left and right cerebral perfusion territories was obtained, allowing discrimination of the anterior and posterior circulation in each hemisphere.

Evidence that OGG1 glycosylase protects neurons against oxidative DNA damage and cell death under ischemic conditions

7,8-Dihydro-8-oxoguanine DNA glycosylase (OGG1) is a major DNA glycosylase involved in base-excision repair (BER) of oxidative DNA damage to nuclear and mitochondrial DNA (mtDNA). We used OGG1-deficient (OGG1(-/-)) mice to examine the possible roles of OGG1 in the vulnerability of neurons to ischemic and oxidative stress. After exposure of cultured neurons to oxidative and metabolic stress levels of OGG1 in the nucleus were elevated and mitochondria exhibited fragmentation and increased levels of the mitochondrial fission protein dynamin-related protein 1 (Drp1) and reduced membrane potential. Cortical neurons isolated from OGG1(-/-) mice were more vulnerable to oxidative insults than were OGG1(+/+) neurons, and OGG1(-/-) mice developed larger cortical infarcts and behavioral deficits after permanent middle cerebral artery occlusion compared with OGG1(+/+) mice. Accumulations of oxidative DNA base lesions (8-oxoG, FapyAde, and FapyGua) were elevated in response to ischemia in both the ipsilateral and contralateral hemispheres, and to a greater extent in the contralateral cortex of OGG1(-/-) mice compared with OGG1(+/+) mice. Ischemia-induced elevation of 8-oxoG incision activity involved increased levels of a nuclear isoform OGG1, suggesting an adaptive response to oxidative nuclear DNA damage. Thus, OGG1 has a pivotal role in repairing oxidative damage to nuclear DNA under ischemic conditions, thereby reducing brain damage and improving functional outcome. Journal of Cerebral Blood Flow & Metabolism (2011) 31, 680-692; doi:10.1038/jcbfm.2010.147; published online 25 August 2010