A mesoscopic modelling approach to anaesthetic action on brain electrical activity

Relating the measurable, large scale, effects of anaesthetic agents to their molecular and cellular targets of action is necessary to better understand the principles by which they affect behavior, as well as enabling the design and evaluation of more effective agents and the better clinical monitoring of existing and future drugs. Volatile and intravenous general anaesthetic agents (GAs) are now known to exert their effects on a variety of protein targets, the most important of which seem to be the neuronal ion channels. It is hence unlikely that anaesthetic effect is the result of a unitary mechanism at the single cell level. However, by altering the behavior of ion channels GAs are believed to change the overall dynamics of distributed networks of neurons. This disruption of regular network activity can be hypothesized to cause the hypnotic and analgesic effects of GAs and may well present more stereotypical characteristics than its underlying microscopic causes. Nevertheless, there have been surprisingly few theories that have attempted to integrate, in a quantitative manner, the empirically well documented alterations in neuronal ion channel behavior with the corresponding macroscopic effects. Here we outline one such approach, and show that a range of well documented effects of anaesthetics on the electroencephalogram (EEG) may be putatively accounted for. In particular we parameterize, on the basis of detailed empirical data, the effects of halogenated volatile ethers (a clinically widely used class of general anaesthetic agent). The resulting model is able to provisionally account for a range of anaesthetically induced EEG phenomena that include EEG slowing, biphasic changes in EEG power, and the dose dependent appearance of anomalous ictal activity, as well as providing a basis for novel approaches to monitoring brain function in both health and disease.

Bifurcations and state changes in the human alpha rhythm: Theory and experiment

Despite many decades investigating scalp recordable 8–13-Hz (alpha) electroencephalographic activity, no consensus has yet emerged regarding its physiological origins nor its functional role in cognition. Here we outline a detailed, physiologically meaningful, theory for the genesis of this rhythm that may provide important clues to its functional role. In particular we find that electroencephalographically plausible model dynamics, obtained with physiological admissible parameterisations, reveals a cortex perched on the brink of stability, which when perturbed gives rise to a range of unanticipated complex dynamics that include 40-Hz (gamma) activity. Preliminary experimental evidence, involving the detection of weak nonlinearity in resting EEG using an extension of the well-known surrogate data method, suggests that nonlinear (deterministic) dynamics are more likely to be associated with weakly damped alpha activity. Thus rather than the “alpha rhythm” being an idling rhythm it may be more profitable to conceive it as a readiness rhythm.

Population based models of cortical drug response: insights from anaesthesia

A great explanatory gap lies between the molecular pharmacology of psychoactive agents and the neurophysiological changes they induce, as recorded by neuroimaging modalities. Causally relating the cellular actions of psychoactive compounds to their influence on population activity is experimentally challenging. Recent developments in the dynamical modelling of neural tissue have attempted to span this explanatory gap between microscopic targets and their macroscopic neurophysiological effects via a range of biologically plausible dynamical models of cortical tissue. Such theoretical models allow exploration of neural dynamics, in particular their modification by drug action. The ability to theoretically bridge scales is due to a biologically plausible averaging of cortical tissue properties. In the resulting macroscopic neural field, individual neurons need not be explicitly represented (as in neural networks). The following paper aims to provide a non-technical introduction to the mean field population modelling of drug action and its recent successes in modelling anaesthesia.

Understanding the effects of anesthetic agents on the EEG through neural field theory

Anesthetic and analgesic agents act through a diverse range of pharmacological mechanisms. Existing empirical data clearly shows that such "microscopic" pharmacological diversity is reflected in their "macroscopic" effects on the human electroencephalogram (EEG). Based on a detailed mesoscopic neural field model we theoretically posit that anesthetic induced EEG activity is due to selective parametric changes in synaptic efficacy and dynamics. Specifically, on the basis of physiologically constrained modeling, it is speculated that the selective modification of inhibitory or excitatory synaptic activity may differentially effect the EEG spectrum. Such results emphasize the importance of neural field theories of brain electrical activity for elucidating the principles whereby pharmacological agents effect the EEG. Such insights will contribute to improved methods for monitoring depth of anesthesia using the EEG.

Reducing inconsistencies in point observations of maximum flood inundation level

Flood simulation models and hazard maps are only as good as the underlying data against which they are calibrated and tested. However, extreme flood events are by definition rare, so the observational data of flood inundation extent are limited in both quality and quantity. The relative importance of these observational uncertainties has increased now that computing power and accurate lidar scans make it possible to run high-resolution 2D models to simulate floods in urban areas. However, the value of these simulations is limited by the uncertainty in the true extent of the flood. This paper addresses that challenge by analyzing a point dataset of maximum water extent from a flood event on the River Eden at Carlisle, United Kingdom, in January 2005. The observation dataset is based on a collection of wrack and water marks from two postevent surveys. A smoothing algorithm for identifying, quantifying, and reducing localized inconsistencies in the dataset is proposed and evaluated showing positive results. The proposed smoothing algorithm can be applied in order to improve flood inundation modeling assessment and the determination of risk zones on the floodplain.

Evaluation of coupled ocean-atmosphere simulations of northern hemisphere extratropical climates in the mid-Holocene

We have used the BIOME4 biogeography–biochemistry model and comparison with palaeovegetation data to evaluate the response of six ocean–atmosphere general circulation models to mid-Holocene changes in orbital forcing in the mid- to high-latitudes of the northern hemisphere. All the models produce: (a) a northward shift of the northern limit of boreal forest, in response to simulated summer warming in high-latitudes. The northward shift is markedly asymmetric, with larger shifts in Eurasia than in North America; (b) an expansion of xerophytic vegetation in mid-continental North America and Eurasia, in response to increased temperatures during the growing season; (c) a northward expansion of temperate forests in eastern North America, in response to simulated winter warming. The northward shift of the northern limit of boreal forest and the northward expansion of temperate forests in North America are supported by palaeovegetation data. The expansion of xerophytic vegetation in mid-continental North America is consistent with palaeodata, although the extent may be over-estimated. The simulated expansion of xerophytic vegetation in Eurasia is not supported by the data. Analysis of an asynchronous coupling of one model to an equilibrium-vegetation model suggests vegetation feedback exacerbates this mid-continental drying and produces conditions more unlike the observations. Not all features of the simulations are robust: some models produce winter warming over Europe while others produce winter cooling. As a result, some models show a northward shift of temperate forests (consistent with, though less marked than, the expansion shown by data) and others produce a reduction in temperate forests. Elucidation of the cause of such differences is a focus of the current phase of the Palaeoclimate Modelling Intercomparison Project.

Using coordinated observations in polarized white light and Faraday rotation to probe the spatial position and magnetic field of an interplanetary sheath

Coronal mass ejections (CMEs) can be continuously tracked through a large portion of the inner heliosphere by direct imaging in visible and radio wavebands. White light (WL) signatures of solar wind transients, such as CMEs, result from Thomson scattering of sunlight by free electrons and therefore depend on both viewing geometry and electron density. The Faraday rotation (FR) of radio waves from extragalactic pulsars and quasars, which arises due to the presence of such solar wind features, depends on the line-of-sight magnetic field component B ∥ and the electron density. To understand coordinated WL and FR observations of CMEs, we perform forward magnetohydrodynamic modeling of an Earth-directed shock and synthesize the signatures that would be remotely sensed at a number of widely distributed vantage points in the inner heliosphere. Removal of the background solar wind contribution reveals the shock-associated enhancements in WL and FR. While the efficiency of Thomson scattering depends on scattering angle, WL radiance I decreases with heliocentric distance r roughly according to the expression Ir –3. The sheath region downstream of the Earth-directed shock is well viewed from the L4 and L5 Lagrangian points, demonstrating the benefits of these points in terms of space weather forecasting. The spatial position of the main scattering site r sheath and the mass of plasma at that position M sheath can be inferred from the polarization of the shock-associated enhancement in WL radiance. From the FR measurements, the local B ∥sheath at r sheath can then be estimated. Simultaneous observations in polarized WL and FR can not only be used to detect CMEs, but also to diagnose their plasma and magnetic field properties.

Association of vitamin D status with arterial blood pressure and hypertension risk: a mendelian randomisation study

BACKGROUND: Low plasma 25-hydroxyvitamin D (25[OH]D) concentration is associated with high arterial blood pressure and hypertension risk, but whether this association is causal is unknown. We used a mendelian randomisation approach to test whether 25(OH)D concentration is causally associated with blood pressure and hypertension risk. METHODS: In this mendelian randomisation study, we generated an allele score (25[OH]D synthesis score) based on variants of genes that affect 25(OH)D synthesis or substrate availability (CYP2R1 and DHCR7), which we used as a proxy for 25(OH)D concentration. We meta-analysed data for up to 108 173 individuals from 35 studies in the D-CarDia collaboration to investigate associations between the allele score and blood pressure measurements. We complemented these analyses with previously published summary statistics from the International Consortium on Blood Pressure (ICBP), the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, and the Global Blood Pressure Genetics (Global BPGen) consortium. FINDINGS: In phenotypic analyses (up to n=49 363), increased 25(OH)D concentration was associated with decreased systolic blood pressure (β per 10% increase, -0·12 mm Hg, 95% CI -0·20 to -0·04; p=0·003) and reduced odds of hypertension (odds ratio [OR] 0·98, 95% CI 0·97-0·99; p=0·0003), but not with decreased diastolic blood pressure (β per 10% increase, -0·02 mm Hg, -0·08 to 0·03; p=0·37). In meta-analyses in which we combined data from D-CarDia and the ICBP (n=146 581, after exclusion of overlapping studies), each 25(OH)D-increasing allele of the synthesis score was associated with a change of -0·10 mm Hg in systolic blood pressure (-0·21 to -0·0001; p=0·0498) and a change of -0·08 mm Hg in diastolic blood pressure (-0·15 to -0·02; p=0·01). When D-CarDia and consortia data for hypertension were meta-analysed together (n=142 255), the synthesis score was associated with a reduced odds of hypertension (OR per allele, 0·98, 0·96-0·99; p=0·001). In instrumental variable analysis, each 10% increase in genetically instrumented 25(OH)D concentration was associated with a change of -0·29 mm Hg in diastolic blood pressure (-0·52 to -0·07; p=0·01), a change of -0·37 mm Hg in systolic blood pressure (-0·73 to 0·003; p=0·052), and an 8·1% decreased odds of hypertension (OR 0·92, 0·87-0·97; p=0·002). INTERPRETATION: Increased plasma concentrations of 25(OH)D might reduce the risk of hypertension. This finding warrants further investigation in an independent, similarly powered study.

The Holocene temperature conundrum

A recent temperature reconstruction of global annual temperature shows Early Holocene warmth followed by a cooling trend through the Middle to Late Holocene [Marcott SA, et al., 2013, Science 339(6124):1198–1201]. This global cooling is puzzling because it is opposite from the expected and simulated global warming trend due to the retreating ice sheets and rising atmospheric greenhouse gases. Our critical reexamination of this contradiction between the reconstructed cooling and the simulated warming points to potentially significant biases in both the seasonality of the proxy reconstruction and the climate sensitivity of current climate models.

Clonal expansion of early to mid-life mitochondrial DNA point mutations drives mitochondrial dysfunction during human ageing

Age-related decline in the integrity of mitochondria is an important contributor to the human ageing process. In a number of ageing stem cell populations, this decline in mitochondrial function is due to clonal expansion of individual mitochondrial DNA (mtDNA) point mutations within single cells. However the dynamics of this process and when these mtDNA mutations occur initially are poorly understood. Using human colorectal epithelium as an exemplar tissue with a well-defined stem cell population, we analysed samples from 207 healthy participants aged 17-78 years using a combination of techniques (Random Mutation Capture, Next Generation Sequencing and mitochondrial enzyme histochemistry), and show that: 1) non-pathogenic mtDNA mutations are present from early embryogenesis or may be transmitted through the germline, whereas pathogenic mtDNA mutations are detected in the somatic cells, providing evidence for purifying selection in humans, 2) pathogenic mtDNA mutations are present from early adulthood (<20 years of age), at both low levels and as clonal expansions, 3) low level mtDNA mutation frequency does not change significantly with age, suggesting that mtDNA mutation rate does not increase significantly with age, and 4) clonally expanded mtDNA mutations increase dramatically with age. These data confirm that clonal expansion of mtDNA mutations, some of which are generated very early in life, is the major driving force behind the mitochondrial dysfunction associated with ageing of the human colorectal epithelium.

Efficacy of condensed tannins against larval Hymenolepis diminuta (Cestoda) in vitro and in the intermediate host Tenebrio molitor (Coleoptera) in vivo

Natural anti-parasitic compounds in plants such as condensed tannins (CT) have anthelmintic properties against a range of gastrointestinal nematodes, but for other helminths such effects are unexplored. The aim of this study was to assess the effects of CT from three different plant extracts in a model system employing the rat tapeworm, Hymenolepis diminuta, in its intermediate host, Tenebrio molitor. An in vitro study examined infectivity of H. diminuta cysticercoids (excystation success) isolated from infected beetles exposed to different concentrations of CT extracts from pine bark (PB) (Pinus sps), hazelnut pericarp (HN) (Corylus avellana) or white clover flowers (WC) (Trifolium repens), in comparison with the anthelmintic drug praziquantel (positive control). In the in vitro study, praziquantel and CT from all three plant extracts had dose-dependent inhibitory effects on cysticercoid excystation. The HN extract was most effective at inhibiting excystation, followed by PB and WC. An in vivo study was carried out on infected beetles (measured as cysticercoid establishment) fed different doses of PB, HN and praziquantel. There was a highly significant inhibitory effect of HN on cysticercoid development (p = 0.0002). Overall, CT showed a promising anti-cestodal effect against the metacestode stage of H. diminuta.