Classification of magnetic resonance brain images using wavelets as input to support vector machine and neural network

In this paper. we propose a novel method using wavelets as input to neural network self-organizing maps and support vector machine for classification of magnetic resonance (MR) images of the human brain. The proposed method classifies MR brain images as either normal or abnormal. We have tested the proposed approach using a dataset of 52 MR brain images. Good classification percentage of more than 94% was achieved using the neural network self-organizing maps (SOM) and 98% front support vector machine. We observed that the classification rate is high for a Support vector machine classifier compared to self-organizing map-based approach.

Solution structure of Am2766: A highly hydrophobic d-conotoxin from Conus amadis that inhibits inactivation of brain voltage gated sodium channels

The three-dimensional (3D) NMR solution structure (MeOH) of the highly hydrophobic δ-conotoxin δ-Am2766 from the molluscivorous snail Conus amadis has been determined. Fifteen converged structures were obtained on the basis of 262 distance constraints, 25 torsion-angle constraints, and ten constraints based on disulfide linkages and H-bonds. The root-mean-square deviations (rmsd) about the averaged coordinates of the backbone (N, Cα, C) and (all) heavy atoms were 0.62±0.20 and 1.12±0.23 Å, respectively. The structures determined are of good stereochemical quality, as evidenced by the high percentage (100%) of backbone dihedral angles that occupy favorable and additionally allowed regions of the Ramachandran map. The structure of δ-Am2766 consists of a triple-stranded antiparallel β-sheet, and of four turns. The three disulfides form the classical ‘inhibitory cysteine knot’ motif. So far, only one tertiary structure of a δ-conotoxin has been reported; thus, the tertiary structure of δ-Am2766 is the second such example.Another Conus peptide, Am2735 from C. amadis, has also been purified and sequenced. Am2735 shares 96% sequence identity with δ-Am2766. Unlike δ-Am2766, Am2735 does not inhibit the fast inactivation of Na+ currents in rat brain Nav1.2 Na+ channels at concentrations up to 200 nM.

Tonically Active Kainate Receptors (tKARs) : A Novel Mechanism Regulating Neuronal Function in the Brain

Fast excitatory transmission between neurons in the central nervous system is mainly mediated by L-glutamate acting on ligand gated (ionotropic) receptors. These are further categorized according to their pharmacological properties to AMPA (2-amino-3-(5-methyl-3-oxo-1,2- oxazol-4-yl)propanoic acid), NMDA (N-Methyl-D-aspartic acid) and kainate (KAR) subclasses. In the rat and the mouse hippocampus, development of glutamatergic transmission is most dynamic during the first postnatal weeks. This coincides with the declining developmental expression of the GluK1 subunit-containing KARs. However, the function of KARs during early development of the brain is poorly understood. The present study reveals novel types of tonically active KARs (hereafter referred to as tKARs) which play a central role in functional development of the hippocampal CA3-CA1 network. The study shows for the first time how concomitant pre- and postsynaptic KAR function contributes to development of CA3-CA1 circuitry by regulating transmitter release and interneuron excitability. Moreover, the tKAR-dependent regulation of transmitter release provides a novel mechanism for silencing and unsilencing early synapses and thus shaping the early synaptic connectivity. The role of GluK1-containing KARs was studied in area CA3 of the neonatal hippocampus. The data demonstrate that presynaptic KARs in excitatory synapses to both pyramidal cells and interneurons are tonically activated by ambient glutamate and that they regulate glutamate release differentially, depending on target cell type. At synapses to pyramidal cells these tKARs inhibit glutamate release in a G-protein dependent manner but in contrast, at synapses to interneurons, tKARs facilitate glutamate release. On the network level these mechanisms act together upregulating activity of GABAergic microcircuits and promoting endogenous hippocampal network oscillations. By virtue of this, tKARs are likely to have an instrumental role in the functional development of the hippocampal circuitry. The next step was to investigate the role of GluK1 -containing receptors in the regulation of interneuron excitability. The spontaneous firing of interneurons in the CA3 stratum lucidum is markedly decreased during development. The shift involves tKARs that inhibit medium-duration afterhyperpolarization (mAHP) in these neurons during the first postnatal week. This promotes burst spiking of interneurons and thereby increases GABAergic activity in the network synergistically with the tKAR-mediated facilitation of their excitatory drive. During development the amplitude of evoked medium afterhyperpolarizing current (ImAHP) is dramatically increased due to decoupling tKAR activation and ImAHP modulation. These changes take place at the same time when the endogeneous network oscillations disappear. These tKAR-driven mechanisms in the CA3 area regulate both GABAergic and glutamatergic transmission and thus gate the feedforward excitatory drive to the area CA1. Here presynaptic tKARs to CA1 pyramidal cells suppress glutamate release and enable strong facilitation in response to high-frequency input. Therefore, CA1 synapses are finely tuned to high-frequency transmission; an activity pattern that is common in neonatal CA3-CA1 circuitry both in vivo and in vitro. The tKAR-regulated release probability acts as a novel presynaptic silencing mechanism that can be unsilenced in response to Hebbian activity. The present results shed new light on the mechanisms modulating the early network activity that paves the way for oscillations lying behind cognitive tasks such as learning and memory. Kainate receptor antagonists are already being developed for therapeutic use for instance against pain and migraine. Because of these modulatory actions, tKARs also represent an attractive candidate for therapeutic treatment of developmentally related complications such as learning disabilities.

Functional outcome and health-related quality of life after traumatic brain injury in the framework of the International Classification of Functioning, Disability and Health (ICF)

Traumatic brain injury (TBI) affects people of all ages and is a cause of long-term disability. In recent years, the epidemiological patterns of TBI have been changing. TBI is a heterogeneous disorder with different forms of presentation and highly individual outcome regarding functioning and health-related quality of life (HRQoL). The meaning of disability differs from person to person based on the individual s personality, value system, past experience, and the purpose he or she sees in life. Understanding of all these viewpoints is needed in comprehensive rehabilitation. This study examines the epidemiology of TBI in Finland as well as functioning and HRQoL after TBI, and compares the subjective and objective assessments of outcome. The frame of reference is the International Classification of Functioning, Disability and Health (ICF). The subjects of Study I represent the population of Finnish TBI patients who experienced their first TBI between 1991 and 2005. The 55 Finnish subjects of Studies II and IV participated in the first wave of the international Quality of life after brain injury (QOLIBRI) validation study. The 795 subjects from six language areas of Study III formed the second wave of the QOLIBRI validation study. The average annual incidence of Finnish hospitalised TBI patients during the years 1991-2005 was 101:100 000 in patients who had TBI as the primary diagnosis and did not have a previous TBI in their medical history. Males (59.2%) were at considerably higher risk of getting a TBI than females. The most common external cause of the injury was falls in all age groups. The number of TBI patients ≥ 70 years of age increased by 59.4% while the number of inhabitants older than 70 years increased by 30.3% in the population of Finland during the same time period. The functioning of a sample of 55 persons with TBI was assessed by extracting information from the patients medical documents using the ICF checklist. The most common problems were found in the ICF components of Body Functions (b) and Activities and Participation (d). HRQoL was assessed with the QOLIBRI which showed the highest level of satisfaction on the Emotions, Physical Problems and Daily Life and Autonomy scales. The highest scores were obtained by the youngest participants and participants living independently without the help of other people, and by people who were working. The relationship between the functional outcome and HRQoL was not straightforward. The procedure of linking the QOLIBRI and the GOSE to the ICF showed that these two outcome measures cover the relevant domains of TBI patients functioning. The QOLIBRI provides the patients subjective view, while the GOSE summarises the objective elements of functioning. Our study indicates that there are certain domains of functioning that are not traditionally sufficiently documented but are important for the HRQoL of persons with TBI. This was the finding especially in the domains of interpersonal relationships, social and leisure activities, self, and the environment. Rehabilitation aims to optimize functioning and to minimize the experience of disability among people with health conditions, and it needs to be based on a comprehensive understanding of human functioning. As an integrative model, the ICF may serve as a frame of reference in achieving such an understanding.

Intraspecific variation in brain size and architecture : population divergence and phenotypic plasticity

Brain size and architecture exhibit great evolutionary and ontogenetic variation. Yet, studies on population variation (within a single species) in brain size and architecture, or in brain plasticity induced by ecologically relevant biotic factors have been largely overlooked. Here, I address the following questions: (i) do locally adapted populations differ in brain size and architecture, (ii) can the biotic environment induce brain plasticity, and (iii) do locally adapted populations differ in levels of brain plasticity? In the first two chapters I report large variation in both absolute and relative brain size, as well as in the relative sizes of brain parts, among divergent nine-spined stickleback (Pungitius pungitius) populations. Some traits show habitat-dependent divergence, implying natural selection being responsible for the observed patterns. Namely, marine sticklebacks have relatively larger bulbi olfactorii (chemosensory centre) and telencephala (involved in learning) than pond sticklebacks. Further, I demonstrate the importance of common garden studies in drawing firm evolutionary conclusions. In the following three chapters I show how the social environment and perceived predation risk shapes brain development. In common frog (Rana temporaria) tadpoles, I demonstrate that under the highest per capita predation risk, tadpoles develop smaller brains than in less risky situations, while high tadpole density results in enlarged tectum opticum (visual brain centre). Visual contact with conspecifics induces enlarged tecta optica in nine-spined sticklebacks, whereas when only olfactory cues from conspecifics are available, bulbus olfactorius become enlarged.Perceived predation risk results in smaller hypothalami (complex function) in sticklebacks. Further, group-living has a negative effect on relative brain size in the competition-adapted pond sticklebacks, but not in the predation-adapted marine sticklebacks. Perceived predation risk induces enlargement of bulbus olfactorius in pond sticklebacks, but not in marine sticklebacks who have larger bulbi olfactorii than pond fish regardless of predation. In sum, my studies demonstrate how applying a microevolutionary approach can help us to understand the enormous variation observed in the brains of wild animals a point-of-view which I high-light in the closing review chapter of my thesis.

High level stable expression of rat brain type IIA sodium channel ?-subunit in CHO cells

The neuronal sodium channels are responsible for the rising phase of action potential and are composed of three subunits, of which the alpha-subunit has been shown to be adequate for most of its functional properties. We have stably expressed the rat brain type IIA sodium channel alpha-subunit in CHO cell tine using a CMV promoter-based vector. The expression was confirmed by detecting a 6.5 kb RNA corresponding to sodium channel alpha-subunit using Northern hybridization. The cells stably expressing the alpha-subunit, yield isolated sodium currents of amplitudes greater than 4nA when studied in whole-cell configuration of the patch-clamp technique. The sodium currents are characterized by activation and inactivation properties similar to neuronal sodium channels, and are blocked by the voltage gated sodium channel blocker tetrodotoxin (TTX).

Mesh Simplification Based on Edge Collapsing Could Improve Computational Efficiency in Near Infrared Optical Tomographic Imaging

The diffusion equation-based modeling of near infrared light propagation in tissue is achieved by using finite-element mesh for imaging real-tissue types, such as breast and brain. The finite-element mesh size (number of nodes) dictates the parameter space in the optical tomographic imaging. Most commonly used finite-element meshing algorithms do not provide the flexibility of distinct nodal spacing in different regions of imaging domain to take the sensitivity of the problem into consideration. This study aims to present a computationally efficient mesh simplification method that can be used as a preprocessing step to iterative image reconstruction, where the finite-element mesh is simplified by using an edge collapsing algorithm to reduce the parameter space at regions where the sensitivity of the problem is relatively low. It is shown, using simulations and experimental phantom data for simple meshes/domains, that a significant reduction in parameter space could be achieved without compromising on the reconstructed image quality. The maximum errors observed by using the simplified meshes were less than 0.27% in the forward problem and 5% for inverse problem.

Multivariate Granger causality: an estimation framework based on factorization of the spectral density matrix

Granger causality is increasingly being applied to multi-electrode neurophysiological and functional imaging data to characterize directional interactions between neurons and brain regions. For a multivariate dataset, one might be interested in different subsets of the recorded neurons or brain regions. According to the current estimation framework, for each subset, one conducts a separate autoregressive model fitting process, introducing the potential for unwanted variability and uncertainty. In this paper, we propose a multivariate framework for estimating Granger causality. It is based on spectral density matrix factorization and offers the advantage that the estimation of such a matrix needs to be done only once for the entire multivariate dataset. For any subset of recorded data, Granger causality can be calculated through factorizing the appropriate submatrix of the overall spectral density matrix.

Connectivity analysis of multichannel EEG signals using recurrence based phase synchronization technique

Real world biological systems such as the human brain are inherently nonlinear and difficult to model. However, most of the previous studies have either employed linear models or parametric nonlinear models for investigating brain function. In this paper, a novel application of a nonlinear measure of phase synchronization based on recurrences, correlation between probabilities of recurrence (CPR), to study connectivity in the brain has been proposed. Being non-parametric, this method makes very few assumptions, making it suitable for investigating brain function in a data-driven way. CPR's utility with application to multichannel electroencephalographic (EEG) signals has been demonstrated. Brain connectivity obtained using thresholded CPR matrix of multichannel EEG signals showed clear differences in the number and pattern of connections in brain connectivity between (a) epileptic seizure and pre-seizure and (b) eyes open and eyes closed states. Corresponding brain headmaps provide meaningful insights about synchronization in the brain in those states. K-means clustering of connectivity parameters of CPR and linear correlation obtained from global epileptic seizure and pre-seizure showed significantly larger cluster centroid distances for CPR as opposed to linear correlation, thereby demonstrating the superior ability of CPR for discriminating seizure from pre-seizure. The headmap in the case of focal epilepsy clearly enables us to identify the focus of the epilepsy which provides certain diagnostic value. (C) 2013 Elsevier Ltd. All rights reserved.

Simultaneously targeting inflammatory response and parasite sequestration in brain to treat Experimental Cerebral Malaria

Malaria afflicts around 200 million people annually, with a mortality number close to 600,000. The mortality rate in Human Cerebral Malaria (HCM) is unacceptably high (15-20%), despite the availability of artemisinin-based therapy. An effective adjunct therapy is urgently needed. Experimental Cerebral Malaria (ECM) in mice manifests many of the neurological features of HCM. Migration of T cells and parasite-infected RBCs (pRBCs) into the brain are both necessary to precipitate the disease. We have been able to simultaneously target both these parameters of ECM. Curcumin alone was able to reverse all the parameters investigated in this study that govern inflammatory responses, CD8(+) T cell and pRBC sequestration into the brain and blood brain barrier (BBB) breakdown. But the animals eventually died of anemia due to parasite build-up in blood. However, arteether-curcumin (AC) combination therapy even after the onset of symptoms provided complete cure. AC treatment is a promising therapeutic option for HCM.

Processing of EEG Signals for Study of Coupling in Brain Regions for Eyes Open and Eyes Closed Conditions

This paper deals with processing the EEG signals obtained from 16 spatially arranged electrodes to measure coupling or synchrony between the frontal, parietal, occipital and temporal lobes of the cerebrum under the eyes open and eyes closed conditions. This synchrony was measured using magnitude squared coherence, Short Time Fourier Transform and wavelet based coherences. We found a pattern in the time-frequency coherence as we moved from the nasion to the inion of the subject's head. The coherence pattern obtained from the wavelet approach was found to be far more capable of picking up peaks in coherence with respect to frequency when compared to the regular Fourier based coherence. We detected high synchrony between frontal polar electrodes that is missing in coherence plots between other electrode pairs. The study has potential applications in healthcare.

A Bacterial Foraging Optimization and Learning Automata Based Feature Selection for Motor Imagery EEG Classification

Selection of relevant features is an open problem in Brain-computer interfacing (BCI) research. Sometimes, features extracted from brain signals are high dimensional which in turn affects the accuracy of the classifier. Selection of the most relevant features improves the performance of the classifier and reduces the computational cost of the system. In this study, we have used a combination of Bacterial Foraging Optimization and Learning Automata to determine the best subset of features from a given motor imagery electroencephalography (EEG) based BCI dataset. Here, we have employed Discrete Wavelet Transform to obtain a high dimensional feature set and classified it by Distance Likelihood Ratio Test. Our proposed feature selector produced an accuracy of 80.291% in 216 seconds.

Comparison of Matching Pursuit Algorithm with Other Signal Processing Techniques for Computation of the Time-Frequency Power Spectrum of Brain Signals

Signals recorded from the brain often show rhythmic patterns at different frequencies, which are tightly coupled to the external stimuli as well as the internal state of the subject. In addition, these signals have very transient structures related to spiking or sudden onset of a stimulus, which have durations not exceeding tens of milliseconds. Further, brain signals are highly nonstationary because both behavioral state and external stimuli can change on a short time scale. It is therefore essential to study brain signals using techniques that can represent both rhythmic and transient components of the signal, something not always possible using standard signal processing techniques such as short time fourier transform, multitaper method, wavelet transform, or Hilbert transform. In this review, we describe a multiscale decomposition technique based on an over-complete dictionary called matching pursuit (MP), and show that it is able to capture both a sharp stimulus-onset transient and a sustained gamma rhythm in local field potential recorded from the primary visual cortex. We compare the performance of MP with other techniques and discuss its advantages and limitations. Data and codes for generating all time-frequency power spectra are provided.

Mathematical study of complex networks : brain, internet, and power grid

The dissertation is concerned with the mathematical study of various network problems. First, three real-world networks are considered: (i) the human brain network (ii) communication networks, (iii) electric power networks. Although these networks perform very different tasks, they share similar mathematical foundations. The high-level goal is to analyze and/or synthesis each of these systems from a “control and optimization” point of view. After studying these three real-world networks, two abstract network problems are also explored, which are motivated by power systems. The first one is “flow optimization over a flow network” and the second one is “nonlinear optimization over a generalized weighted graph”. The results derived in this dissertation are summarized below.

Brain Networks: Neuroimaging data reveals the coordinated activity of spatially distinct brain regions, which may be represented mathematically as a network of nodes (brain regions) and links (interdependencies). To obtain the brain connectivity network, the graphs associated with the correlation matrix and the inverse covariance matrix—describing marginal and conditional dependencies between brain regions—have been proposed in the literature. A question arises as to whether any of these graphs provides useful information about the brain connectivity. Due to the electrical properties of the brain, this problem will be investigated in the context of electrical circuits. First, we consider an electric circuit model and show that the inverse covariance matrix of the node voltages reveals the topology of the circuit. Second, we study the problem of finding the topology of the circuit based on only measurement. In this case, by assuming that the circuit is hidden inside a black box and only the nodal signals are available for measurement, the aim is to find the topology of the circuit when a limited number of samples are available. For this purpose, we deploy the graphical lasso technique to estimate a sparse inverse covariance matrix. It is shown that the graphical lasso may find most of the circuit topology if the exact covariance matrix is well-conditioned. However, it may fail to work well when this matrix is ill-conditioned. To deal with ill-conditioned matrices, we propose a small modification to the graphical lasso algorithm and demonstrate its performance. Finally, the technique developed in this work will be applied to the resting-state fMRI data of a number of healthy subjects.

Communication Networks: Congestion control techniques aim to adjust the transmission rates of competing users in the Internet in such a way that the network resources are shared efficiently. Despite the progress in the analysis and synthesis of the Internet congestion control, almost all existing fluid models of congestion control assume that every link in the path of a flow observes the original source rate. To address this issue, a more accurate model is derived in this work for the behavior of the network under an arbitrary congestion controller, which takes into account of the effect of buffering (queueing) on data flows. Using this model, it is proved that the well-known Internet congestion control algorithms may no longer be stable for the common pricing schemes, unless a sufficient condition is satisfied. It is also shown that these algorithms are guaranteed to be stable if a new pricing mechanism is used.

Electrical Power Networks: Optimal power flow (OPF) has been one of the most studied problems for power systems since its introduction by Carpentier in 1962. This problem is concerned with finding an optimal operating point of a power network minimizing the total power generation cost subject to network and physical constraints. It is well known that OPF is computationally hard to solve due to the nonlinear interrelation among the optimization variables. The objective is to identify a large class of networks over which every OPF problem can be solved in polynomial time. To this end, a convex relaxation is proposed, which solves the OPF problem exactly for every radial network and every meshed network with a sufficient number of phase shifters, provided power over-delivery is allowed. The concept of “power over-delivery” is equivalent to relaxing the power balance equations to inequality constraints.

Flow Networks: In this part of the dissertation, the minimum-cost flow problem over an arbitrary flow network is considered. In this problem, each node is associated with some possibly unknown injection, each line has two unknown flows at its ends related to each other via a nonlinear function, and all injections and flows need to satisfy certain box constraints. This problem, named generalized network flow (GNF), is highly non-convex due to its nonlinear equality constraints. Under the assumption of monotonicity and convexity of the flow and cost functions, a convex relaxation is proposed, which always finds the optimal injections. A primary application of this work is in the OPF problem. The results of this work on GNF prove that the relaxation on power balance equations (i.e., load over-delivery) is not needed in practice under a very mild angle assumption.

Generalized Weighted Graphs: Motivated by power optimizations, this part aims to find a global optimization technique for a nonlinear optimization defined over a generalized weighted graph. Every edge of this type of graph is associated with a weight set corresponding to the known parameters of the optimization (e.g., the coefficients). The motivation behind this problem is to investigate how the (hidden) structure of a given real/complex valued optimization makes the problem easy to solve, and indeed the generalized weighted graph is introduced to capture the structure of an optimization. Various sufficient conditions are derived, which relate the polynomial-time solvability of different classes of optimization problems to weak properties of the generalized weighted graph such as its topology and the sign definiteness of its weight sets. As an application, it is proved that a broad class of real and complex optimizations over power networks are polynomial-time solvable due to the passivity of transmission lines and transformers.

Altered expression of Alzheimer's disease-related genes in the cerebellum of autistic patients: a model for disrupted brain connectome and therapy

Autism and Alzheimer's disease (AD) are, respectively, neurodevelopmental and degenerative diseases with an increasing epidemiological burden. The AD-associated amyloid-beta precursor protein-alpha has been shown to be elevated in severe autism, leading to the 'anabolic hypothesis' of its etiology. Here we performed a focused microarray analysis of genes belonging to NOTCH and WNT signaling cascades, as well as genes related to AD and apoptosis pathways in cerebellar samples from autistic individuals, to provide further evidence for pathological relevance of these cascades for autism. By using the limma package from R and false discovery rate, we demonstrated that 31% (116 out of 374) of the genes belonging to these pathways displayed significant changes in expression (corrected P-values <0.05), with mitochondria- related genes being the most downregulated. We also found upregulation of GRIN1, the channel-forming subunit of NMDA glutamate receptors, and MAP3K1, known activator of the JNK and ERK pathways with anti-apoptotic effect. Expression of PSEN2 (presinilin 2) and APBB1 (or F65) were significantly lower when compared with control samples. Based on these results, we propose a model of NMDA glutamate receptor-mediated ERK activation of alpha-secretase activity and mitochondrial adaptation to apoptosis that may explain the early brain overgrowth and disruption of synaptic plasticity and connectome in autism. Finally, systems pharmacology analyses of the model that integrates all these genes together (NOWADA) highlighted magnesium (Mg2+) and rapamycin as most efficient drugs to target this network model in silico. Their potential therapeutic application, in the context of autism, is therefore discussed.