Causality and synchronisation in complex systems with applications to neuroscience

This thesis presents an investigation, of synchronisation and causality, motivated by problems in computational neuroscience. The thesis addresses both theoretical and practical signal processing issues regarding the estimation of interdependence from a set of multivariate data generated by a complex underlying dynamical system. This topic is driven by a series of problems in neuroscience, which represents the principal background motive behind the material in this work. The underlying system is the human brain and the generative process of the data is based on modern electromagnetic neuroimaging methods . In this thesis, the underlying functional of the brain mechanisms are derived from the recent mathematical formalism of dynamical systems in complex networks. This is justified principally on the grounds of the complex hierarchical and multiscale nature of the brain and it offers new methods of analysis to model its emergent phenomena. A fundamental approach to study the neural activity is to investigate the connectivity pattern developed by the brain’s complex network. Three types of connectivity are important to study: 1) anatomical connectivity refering to the physical links forming the topology of the brain network; 2) effective connectivity concerning with the way the neural elements communicate with each other using the brain’s anatomical structure, through phenomena of synchronisation and information transfer; 3) functional connectivity, presenting an epistemic concept which alludes to the interdependence between data measured from the brain network. The main contribution of this thesis is to present, apply and discuss novel algorithms of functional connectivities, which are designed to extract different specific aspects of interaction between the underlying generators of the data. Firstly, a univariate statistic is developed to allow for indirect assessment of synchronisation in the local network from a single time series. This approach is useful in inferring the coupling as in a local cortical area as observed by a single measurement electrode. Secondly, different existing methods of phase synchronisation are considered from the perspective of experimental data analysis and inference of coupling from observed data. These methods are designed to address the estimation of medium to long range connectivity and their differences are particularly relevant in the context of volume conduction, that is known to produce spurious detections of connectivity. Finally, an asymmetric temporal metric is introduced in order to detect the direction of the coupling between different regions of the brain. The method developed in this thesis is based on a machine learning extensions of the well known concept of Granger causality. The thesis discussion is developed alongside examples of synthetic and experimental real data. The synthetic data are simulations of complex dynamical systems with the intention to mimic the behaviour of simple cortical neural assemblies. They are helpful to test the techniques developed in this thesis. The real datasets are provided to illustrate the problem of brain connectivity in the case of important neurological disorders such as Epilepsy and Parkinson’s disease. The methods of functional connectivity in this thesis are applied to intracranial EEG recordings in order to extract features, which characterize underlying spatiotemporal dynamics before during and after an epileptic seizure and predict seizure location and onset prior to conventional electrographic signs. The methodology is also applied to a MEG dataset containing healthy, Parkinson’s and dementia subjects with the scope of distinguishing patterns of pathological from physiological connectivity.

Complex temporal patterns in molecular dynamics:a direct measure of the phase-space exploration by the trajectory at macroscopic time scales

Computer simulated trajectories of bulk water molecules form complex spatiotemporal structures at the picosecond time scale. This intrinsic complexity, which underlies the formation of molecular structures at longer time scales, has been quantified using a measure of statistical complexity. The method estimates the information contained in the molecular trajectory by detecting and quantifying temporal patterns present in the simulated data (velocity time series). Two types of temporal patterns are found. The first, defined by the short-time correlations corresponding to the velocity autocorrelation decay times (â‰0.1â€ps), remains asymptotically stable for time intervals longer than several tens of nanoseconds. The second is caused by previously unknown longer-time correlations (found at longer than the nanoseconds time scales) leading to a value of statistical complexity that slowly increases with time. A direct measure based on the notion of statistical complexity that describes how the trajectory explores the phase space and independent from the particular molecular signal used as the observed time series is introduced. © 2008 The American Physical Society.

Signal processing in a loudspeaking telephone

One of the major problems associated with communication via a loudspeaking telephone (LST) is that, using analogue processing, duplex transmission is limited to low-loss lines and produces a low acoustic output. An architectural for an instrument has been developed and tested, which uses digital signal processing to provide duplex transmission between a LST and a telopnone handset over most of the B.T. network. Digital adaptive-filters are used in the duplex LST to cancel coupling between the loudspeaker and microphone, and across the transmit to receive paths of the 2-to-4-wire converter. Normal movement of a person in the acoustic path causes a loss of stability by increasing the level of coupling from the loudspeaker to the microphone, since there is a lag associated the adaptive filters learning about a non-stationary path, Control of the loop stability and the level of sidetone heard by the hadset user is by a microprocessoe, which continually monitors the system and regulates the gain. The result is a system which offers the best compromise available based on a set of measured parameters.A theory has been developed which gives the loop stability requirements based on the error between the parameters of the filter and those of the unknown path. The programme to develope a low-cost adaptive filter in LST produced a low-cost adaptive filter in LST produced a unique architecture which has a number of features not available in any similar system. These include automatic compensation for the rate of adaptation over a 36 dB range of output level, , 4 rates of adaptation (with a maximum of 465 dB/s), plus the ability to cascade up to 4 filters without loss o performance. A complex story has been developed to determine the adptation which can be achieved using finite-precision arithmatic. This enabled the development of an architecture which distributed the normalisation required to achieve optimum rate of adaptation over the useful input range. Comparison of theory and measurement for the adaptive filter show very close agreement. A single experimental LST was built and tested on connections to hanset telephones over the BT network. The LST demonstrated that duplex transmission was feasible using signal processing and produced a more comfortable means of communication beween people than methods emplying deep voice-switching to regulate the local-loop gain. Although, with the current level of processing power, it is not a panacea and attention must be directed toward the physical acoustic isolation between loudspeaker and microphone.

Spiral mechanisms are required to account for summation of complex motion components

Stimuli from one family of complex motions are defined by their spiral pitch, where cardinal axes represent signed expansion and rotation. Intermediate spirals are represented by intermediate pitches. It is well established that vision contains mechanisms that sum over space and direction to detect these stimuli (Morrone et al., Nature 376 (1995) 507) and one possibility is that four cardinal mechanisms encode the entire family. We extended earlier work (Meese & Harris, Vision Research 41 (2001) 1901) using subthreshold summation of random dot kinematograms and a two-interval forced choice technique to investigate this possibility. In our main experiments, the spiral pitch of one component was fixed and that of another was varied in steps of 15° relative to the first. Regardless of whether the fixed component was aligned with cardinal axes or an intermediate spiral, summation to-coherence-threshold between the two components declined as a function of their difference in spiral pitch. Similar experiments showed that none of the following were critical design features or stimulus parameters for our results: superposition of signal dots, limited life-time dots, the presence of speed gradients, stimulus size or the number of dots. A simplex algorithm was used to fit models containing mechanisms spaced at a pitch of either 90° (cardinal model) or 45° (cardinal+model) and combined using a fourth-root summation rule. For both models, direction half-bandwidth was equated for all mechanisms and was the only free parameter. Only the cardinal+model could account for the full set of results. We conclude that the detection of complex motion in human vision requires both cardinal and spiral mechanisms with a half-bandwidth of approximately 46°. © 2002 Elsevier Science Ltd. All rights reserved.

Broad direction bandwidths for complex motion mechanisms

Growing evidence from psychophysics and single-unit recordings suggests specialised mechanisms in the primate visual system for the detection of complex motion patterns such as expansion and rotation. Here we used a subthreshold summation technique to determine the direction tuning functions of the detecting mechanisms. We measured thresholds for discriminating noise and signal + noise for pairs of superimposed complex motion patterns (signal A and B) carried by random-dot stimuli in a circular 5° field. For expansion, rotation, deformation and translation we found broad tuning functions approximated by cos(d), where d is the difference in dot directions for signal A and B. These data were well described by models in which either: (a) cardinal mechanisms had direction bandwidths (half-widths) of around 60° or (b) the number of mechanisms was increased and their half-width was reduced to about 40°. When d = 180° we found summation to be greater than probability summation for expansion, rotation and translation, consistent with the idea that mechanisms for these stimuli are constructed from subunits responsive to relative motion. For deformation, however, we found sensitivity declined when d = 180°, suggesting antagonistic input from directional subunits in the deformation mechanism. This is a necessary property for a mechanism whose job is to extract the deformation component from the optic flow field. © 2001 Elsevier Science Ltd.

Discourse-analytic approaches to text and talk

This chapter explores the different ways in which discourse-analytic approaches reveal the ‘meaningfulness’ of text and talk. It reviews four diverse approaches to discourse analysis of particular value for current research in linguistics: Conversation Analysis (CA), Discourse Analysis (DA), Critical Discourse Analysis (CDA) and Feminist Post-structuralist Discourse Analysis (FPDA). Each approach is examined in terms of its background, motivation, key features, and possible strengths and limitations in relation to the field of linguistics. A key way to schematize discourse-analytic methodology is in terms of its relationship between microanalytical approaches, which examine the finer detail of linguistic interactions in transcripts, and macroanalytical approaches, which consider how broader social processes work through language (Heller, 2001). This chapter assesses whether there is a strength in a discourse-analytic approach that aligns itself exclusively with either a micro- or macrostrategy, or whether, as Heller suggests, the field needs to fi nd a way of ‘undoing’ the micro–macro dichotomy in order to produce richer, more complex insights within linguistic research.

Enhanced information transmission with signal dependent noise in an array of LIF neurons

Sensory cells usually transmit information to afferent neurons via chemical synapses, in which the level of noise is dependent on an applied stimulus. Taking into account such dependence, we model a sensory system as an array of LIF neurons with a common signal. We show that information transmission is enhanced by a nonzero level of noise. Moreover, we demonstrate a phenomenon similar to suprathreshold stochastic resonance with additive noise. We remark that many properties of information transmission found for the LIF neurons was predicted by us before with simple binary units [Phys. Rev. E 75, 021121 (2007)]. This confirmation of our predictions allows us to point out identical roots of the phenomena found in the simple threshold systems and more complex LIF neurons.

Ceramide and reactive oxygen species (ROS) as signal transduction molecules in inflammation

Reactive oxygen species (ROS) and the sphingolipid ceramide are each partly responsible for the intracellular signal transduction of a variety of physiological, pharmacological or environmental agents. Furthermore, the enhanced production of many of these agents, that utilise ROS and ceramide as signalling intermediates, is associated with the aetiologies of several vascular diseases (e.g. atherosclerosis) or disorders of inflammatory origin (e.g. rheumatoid arthritis; RA). Excessive monocyte recruitment and uncontrolled T cell activation are both strongly implicated in the chronic inflammatory responses that are associated with these pathologies. Therefore the aims of this thesis are (1) to further elucidate the cellular responses to modulations in intracellular ceramide/ROS levels in monocytes and T cells, in order to help resolve the mechanisms of progression of these diseases and (2) to examine both existing agents (methotrexate) and novel targets for possible therapeutic manipulation. Utilising synthetic, short chain ceramide to mimic the cellular responses to fluctuations in natural endogenous ceramide or, stimulation of CD95 to induce ceramide formation, it is described here that ceramide targets and manipulates two discrete sites responsible for ROS generation, preceding the cellular responses of growth arrest in U937 monocytes and apoptosis in Jurkat T-cells. In both cell types, transient elevations in mitochondrial ROS generation were observed. However, the prominent redox altering effects appear to be the ceramide-mediated reduction in cytosolic peroxide, the magnitude of which dictates in part the cellular response in U937 monocytes, Jurkat T-cells and primary human peripheral blood resting or PHA-activated T-cells in vitro. The application of synthetic ceramides to U937 monocytes for short (2 hours) or long (16 hours) treatment periods reduced the membrane expression of proteins associated with cell-cell interaction. Furthermore, ceramide treated U937 monocytes demonstrated reduced adhesion to 5 or 24 hour LPS activated human umbilical vein endothelial cells (HUVEC) but not resting HUVEC. Consequently it is hypothesised that the targeted treatment of monocytes from patients with cardiovascular diseases with short chain synthetic ceramide may reduce disease progression. Herein, the anti-inflammatory and immunosuppressant drug, methotrexate, is described to require ROS production for the induction of cytostasis or cytotoxicity in U937 monocytes and Jurkat T-cells respectively. Further, ROS are critical for methotrexate to abrogate monocyte interaction with activated HUVEC in vitro. The histological feature of RA of enhanced infiltration, survivability and hyporesponsiveness of T-cells within the diseased synovium has been suggested to arise from aberrant signalling. No difference in the concentrations of endogenous T-cell ceramide, the related lipid diacylglycerol (DAG) and cytosolic peroxide ex vivo was observed. TCR activation following PHA exposure in vitro for 72 hours did not induced maintained perturbations in DAG or ceramide in T-cells from RA patients or healthy individuals. However, T-cells from RA patients failed to upregulate cytosolic peroxide in response to PHA, unlike those from normals, despite expressing identical levels of the activation marker CD25. This inability to upregulate cytosolic peroxide may contribute to the T-cell pathology associated with RA by affecting the signalling capacity of redox sensitive biomolecules. These data highlight the importance of two distinctive cellular pools of ROS in mediating complex biological events associated with inflammatory disease and suggest that modulation of cellular ceramides represents a novel therapeutic strategy to minimise monocyte recruitment.

Contrast sensitivity for complex and random gratings

This thesis studied the effect of (i) the number of grating components and (ii) parameter randomisation on root-mean-square (r.m.s.) contrast sensitivity and spatial integration. The effectiveness of spatial integration without external spatial noise depended on the number of equally spaced orientation components in the sum of gratings. The critical area marking the saturation of spatial integration was found to decrease when the number of components increased from 1 to 5-6 but increased again at 8-16 components. The critical area behaved similarly as a function of the number of grating components when stimuli consisted of 3, 6 or 16 components with different orientations and/or phases embedded in spatial noise. Spatial integration seemed to depend on the global Fourier structure of the stimulus. Spatial integration was similar for sums of two vertical cosine or sine gratings with various Michelson contrasts in noise. The critical area for a grating sum was found to be a sum of logarithmic critical areas for the component gratings weighted by their relative Michelson contrasts. The human visual system was modelled as a simple image processor where the visual stimuli is first low-pass filtered by the optical modulation transfer function of the human eye and secondly high-pass filtered, up to the spatial cut-off frequency determined by the lowest neural sampling density, by the neural modulation transfer function of the visual pathways. The internal noise is then added before signal interpretation occurs in the brain. The detection is mediated by a local spatially windowed matched filter. The model was extended to include complex stimuli and its applicability to the data was found to be successful. The shape of spatial integration function was similar for non-randomised and randomised simple and complex gratings. However, orientation and/or phase randomised reduced r.m.s contrast sensitivity by a factor of 2. The effect of parameter randomisation on spatial integration was modelled under the assumption that human observers change the observer strategy from cross-correlation (i.e., a matched filter) to auto-correlation detection when uncertainty is introduced to the task. The model described the data accurately.

GP130/OSMR is the only LIF/IL-6 family receptor complex to promote osteoblast differentiation of calvaria progenitors

Leukemia inhibitory factor (LIF) and its receptor (LIFR) are "twins" of Oncostatin M (OSM) and OSMR, respectively, likely having arisen through gene duplications. We compared their effects in a bone nodule-forming model of in vitro osteogenesis, rat calvaria (RC) cell cultures. Using a dominant-negative LIF mutant (hLIF-05), we showed that in RC cell cultures mouse OSM (mOSM) activates exclusively glycoprotein 130 (gp130)/OSMR. In treatments starting at early nodule formation stage, LIF, mOSM, IL-11, and IL-6 + sIL-6R inhibit bone nodule formation, that is, osteoprogenitor differentiation. Treatment with mOSM, and no other cytokine of the family, in early cultures (day 1-3 or 1-4) increases bone colony numbers. hLIF-05 also dose dependently stimulates bone nodule formation, confirming the inhibitory action of gp130/LIFR on osteogenesis. In pulse treatments at successive stages of bone nodule formation and maturation, LIF blocks osteocalcin (OCN) expression by differentiated osteoblasts, but has no effect on bonesialoprotein (BSP) expression. Mouse OSM inhibits OCN and BSP expression in preconfluent cultures with no or progressively reduced effects at later stages, reflecting the disruption of early nodules, possibly due to the strong apoptotic action of mOSM in RC cell cultures. In summary, LIFR and OSMR display differential effects on differentiation and phenotypic expression of osteogenic cells, most likely through different signal transduction pathways. In particular, gp130/OSMR is the only receptor complex of the family to stimulate osteoprogenitor differentiation in the RC cell culture model. © 2005 Wiley-Liss, Inc.

Managing project risk using combined analytic hierarchy process and risk map

The main purpose of the study is to develop an integrated framework for managing project risks by analyzing risk across project, work package and activity levels, and developing responses. Design/methodology/approach: The study first reviews the literature of various contemporary risk management frameworks in order to identify gaps in project risk management knowledge. Then it develops a conceptual risk management framework using combined analytic hierarchy process (AHP) and risk map for managing project risks. The proposed framework has then been applied to a 1500 km oil pipeline construction project in India in order to demonstrate its effectiveness. The concerned project stakeholders were involved through focus group discussions for applying the proposed risk management framework in the project under study. Findings: The combined AHP and risk map approach is very effective to manage project risks across project, work package and activity levels. The risk factors in project level are caused because of external forces such as business environment (e.g. customers, competitors, technological development, politics, socioeconomic environment). The risk factors in work package and activity levels are operational in nature and created due to internal causes such as lack of material and labor productivity, implementation issues, team ineffectiveness, etc. Practical implications: The suggested model can be applied to any complex project and helps manage risk throughout the project life cycle. Originality/value: Both business and operational risks constitute project risks. In one hand, the conventional project risk management frameworks emphasize on managing business risks and often ignore operational risks. On the other hand, the studies that deal with operational risk often do not link them with business risks. However, they need to be addressed in an integrated way as there are a few risks that affect only the specific level. Hence, this study bridges the gaps. © 2010 Elsevier B.V. All rights reserved.