Analyzing multiple spike trains with nonparametric granger causality

Simultaneous recordings of spike trains from multiple single neurons are becoming commonplace. Understanding the interaction patterns among these spike trains remains a key research area. A question of interest is the evaluation of information flow between neurons through the analysis of whether one spike train exerts causal influence on another. For continuous-valued time series data, Granger causality has proven an effective method for this purpose. However, the basis for Granger causality estimation is autoregressive data modeling, which is not directly applicable to spike trains. Various filtering options distort the properties of spike trains as point processes. Here we propose a new nonparametric approach to estimate Granger causality directly from the Fourier transforms of spike train data. We validate the method on synthetic spike trains generated by model networks of neurons with known connectivity patterns and then apply it to neurons limultaneously recorded from the thalamus and the primary somatosensory cortex of a squirrel monkey undergoing tactile stimulation.

Denoising neural data with state-space smoothing: Method and application

Neural data are inevitably contaminated by noise. When such noisy data are subjected to statistical analysis, misleading conclusions can be reached. Here we attempt to address this problem by applying a state-space smoothing method, based on the combined use of the Kalman filter theory and the Expectation–Maximization algorithm, to denoise two datasets of local field potentials recorded from monkeys performing a visuomotor task. For the first dataset, it was found that the analysis of the high gamma band (60–90 Hz) neural activity in the prefrontal cortex is highly susceptible to the effect of noise, and denoising leads to markedly improved results that were physiologically interpretable. For the second dataset, Granger causality between primary motor and primary somatosensory cortices was not consistent across two monkeys and the effect of noise was suspected. After denoising, the discrepancy between the two subjects was significantly reduced.

Fast Robust Pattern Classification Algorithms for Real Time Neuro-Motor Prosthetic Applications

In this paper, we give a brief review of pattern classification algorithms based on discriminant analysis. We then apply these algorithms to classify movement direction based on multivariate local field potentials recorded from a microelectrode array in the primary motor cortex of a monkey performing a reaching task. We obtain prediction accuracies between 55% and 90% using different methods which are significantly above the chance level of 12.5%.

Effect of amplitude correlations on coherence in the local field potential

Neural activity across the brain shows both spatial and temporal correlations at multiple scales, and understanding these correlations is a key step toward understanding cortical processing. Correlation in the local field potential (LFP) recorded from two brain areas is often characterized by computing the coherence, which is generally taken to reflect the degree of phase consistency across trials between two sites. Coherence, however, depends on two factors-phase consistency as well as amplitude covariation across trials-but the spatial structure of amplitude correlations across sites and its contribution to coherence are not well characterized. We recorded LFP from an array of microelectrodes chronically implanted in the primary visual cortex of monkeys and studied correlations in amplitude across electrodes as a function of interelectrode distance. We found that amplitude correlations showed a similar trend as coherence as a function of frequency and interelectrode distance. Importantly, even when phases were completely randomized between two electrodes, amplitude correlations introduced significant coherence. To quantify the contributions of phase consistency and amplitude correlations to coherence, we simulated pairs of sinusoids with varying phase consistency and amplitude correlations. These simulations confirmed that amplitude correlations can significantly bias coherence measurements, resulting in either over-or underestimation of true phase coherence. Our results highlight the importance of accounting for the correlations in amplitude while using coherence to study phase relationships across sites and frequencies.

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.

Mutations in WDR62, encoding a centrosomal and nuclear protein, in Indian primary microcephaly families with cortical malformations

Primary microcephaly is an autosomal recessive disorder characterized by smaller than normal brain size and mental retardation. It is genetically heterogeneous with seven loci: MCPH1-MCPH7. We have previously reported genetic analysis of 35 families, including the identification of the MCPH7 gene STIL. Of the 35 families, three families showed linkage to the MCPH2 locus. Recent whole-exome sequencing studies have shown that the WDR62 gene, located in the MCPH2 candidate region, is mutated in patients with severe brain malformations. We therefore sequenced the WDR62 gene in our MCPH2 families and identified two novel homozygous protein truncating mutations in two families. Affected individuals in the two families had pachygyria, microlissencephaly, band heterotopias, gyral thickening, and dysplastic cortex. Using immunofluorescence study, we showed that, as with other MCPH proteins, WDR62 localizes to centrosomes in A549, HepG2, and HaCaT cells. In addition, WDR62 was also localized to nucleoli. Bioinformatics analysis predicted two overlapping nuclear localization signals and multiple WD-40 repeats in WDR62. Two other groups have also recently identified WDR62 mutations in MCPH2 families. Our results therefore add further evidence that WDR62 is the MCPH2 gene. The present findings will be helpful in genetic diagnosis of patients linked to the MCPH2 locus.

Mutations in STIL, Encoding a Pericentriolar and Centrosomal Protein, Cause Primary Microcephaly

Primary microcephaly (MCPH) is an autosomal-recessive congenital disorder characterized by smaller-than-normal brain size and mental retardation. MCPH is genetically heterogeneous with six known loci: MCPH1-MCPH6. We report mapping of a novel locus, MCPH7, to chromosome 1p32.3-p33 between markers D1S2797 and D1S417, corresponding to a physical distance of 8.39 Mb. Heterogeneity analysis of 24 families previously excluded from linkage to the six known MCPH loci suggested linkage of five families (20.83%) to the MCPH7 locus. In addition, four families were excluded from linkage to the MCPH7 locus as well as all of the six previously known loci, whereas the remaining 15 families could not be conclusively excluded or included. The combined maximum two-point LOD score for the linked families was 5.96 at marker D1S386 at theta = 0.0. The combined multipoint LOD score was 6.97 between markers D1S2797 and D1S417. Previously, mutations in four genes, MCPH1, CDK5RAP2, ASPM, and CENPJ, that code for centrosomal proteins have been shown to cause this disorder. Three different homozygous mutations in STIL, which codes for a pericentriolar and centrosomal protein, were identified in patients from three of the five families linked to the MCPH7 locus; all are predicted to truncate the STIL protein. Further, another recently ascertained family was homozygous for the same mutation as one of the original families. There was no evidence for a common haplotype. These results suggest that the centrosome and its associated structures are important in the control of neurogenesis in the developing human brain.

The nature of non-equilibrium solidification of undercooled Agsingle bondCu alloys in contact with primary copper dendrites

Liquids of silver-copper alloys with near eutectic compositions embedded in a copper matrix were undercooled. The structural and microstructural investigations of these alloys solidified from undercooled temperature indicated the absence of both the eutectic reaction and diffusionless transformation below the equal free energy curve (T0). Instead the liquid maintained local equilibrium with the copper dendrites continuously until it intersected the extended solidus of the silver rich solid solution.

Molecular vibrations and conformation of primary dithiocarbamate ester. Infrared and NMR studies on S-methyl dithiocarbamate

The i.r. spectra of a primary dithiocarbamate ester namely, S-methyl dithiocarbamate (SMDTC) and its N-dideuterated compound have been measured between 4000 and 30 cm−1. Spectra in solution and at liquid nitrogen temperature have also been obtained. Assignment of all the fundamentals has been proposed and supported from a full normal coordinate analysis. The band assignments for SMDTC have been compared with those of related molecules and the characteristic bands of primary thioamides are derived. Conformational flexibility of SMDTC has been examined by i.r. and proton NMR spectroscopy. The hindered rotation around the C---N bond has been studied by a complete line shape analysis. The magnitude of ---NH2 and ---CH3 torsional barriers is also estimated from vibrational frequencies and force constants.

Primary structure of a Thomsen-Friedenreich-antigen-specific lectin, jacalin [Artocarpus integrifolia (jack fruit) agglutinin]. Evidence for the presence of an internal repeat

Jacalin [Artocarpus integrifolia (jack fruit) agglutinin] is made up of two types of chains, heavy and light, with M(r) values of 16,200 +/- 1200 and 2090 +/- 300 respectively (on the basis of gel-permeation chromatography under denaturing conditions). Its complete amino acid sequence was determined by manual degradation using a 4-dimethylaminoazobenzene 4'-isothiocyanate double-coupling method. Peptide fragments for sequence analysis were obtained by chemical cleavages of the heavy chain with CNBr, hydroxylamine hydrochloride and iodosobenzoic acid and enzymic cleavage with Staphylococcus aureus proteinase. The peptides were purified by a combination gel-permeation and reverse-phase chromatography. The light chains, being only 20 residues long, could be sequenced without fragmentation. Amino acid analyses and carboxypeptidase-Y-digestion C-terminal analyses of the subunits provided supportive evidence for their sequence. Computer-assisted alignment of the jacalin heavy-chain sequence failed to show sequence similarity to that of any lectin for which the complete sequence is known. Analyses of the sequence showed the presence of an internal repeat spanning residues 7-64 and 76-130. The internal repeat was found to be statistically significant.

Regulation of extracellular matrix genes by arecoline in primary gingival fibroblasts requires epithelial factors

Background and Objective: Oral submucous fibrosis, a disease of collagen disorder, has been attributed to arecoline present in the saliva of betel quid chewers. However, the molecular basis of the action of arecoline in the pathogenesis of oral submucous fibrosis is poorly understood. The basic aim of our study was to elucidate the mechanism underlying the action of arecoline on the expression of genes in oral fibroblasts. Material and Methods: Human keratinocytes (HaCaT cells) and primary human gingival fibroblasts were treated with arecoline in combination with various pathway inhibitors, and the expression of transforming growth factor-beta isoform genes and of collagen isoforms was assessed using reverse transcription polymerase chain reaction analysis. Results: We observed the induction of transforming growth factor-beta2 by arecoline in HaCaT cells and this induction was found to be caused by activation of the M-3 muscarinic acid receptor via the induction of calcium and the protein kinase C pathway. Most importantly, we showed that transforming growth factor-beta2 was significantly overexpressed in oral submucous fibrosis tissues (p = 0.008), with a median of 2.13 (n = 21) compared with 0.75 (n = 18) in normal buccal mucosal tissues. Furthermore, arecoline down-regulated the expression of collagens 1A1 and 3A1 in human primary gingival fibroblasts; however these collagens were induced by arecoline in the presence of spent medium of cultured human keratinocytes. Treatment with a transforming growth factor-beta blocker, transforming growth factor-beta1 latency-associated peptide, reversed this up-regulation of collagen, suggesting a role for profibrotic cytokines, such as transforming growth factor-beta, in the induction of collagens. Conclusion: Taken together, our data highlight the importance of arecoline-induced epithelial changes in the pathogenesis of oral submucous fibrosis.

Prediction of Candidate Primary Immunodeficiency Disease Genes Using a Support Vector Machine Learning Approach

Screening and early identification of primary immunodeficiency disease (PID) genes is a major challenge for physicians. Many resources have catalogued molecular alterations in known PID genes along with their associated clinical and immunological phenotypes. However, these resources do not assist in identifying candidate PID genes. We have recently developed a platform designated Resource of Asian PDIs, which hosts information pertaining to molecular alterations, protein-protein interaction networks, mouse studies and microarray gene expression profiling of all known PID genes. Using this resource as a discovery tool, we describe the development of an algorithm for prediction of candidate PID genes. Using a support vector machine learning approach, we have predicted 1442 candidate PID genes using 69 binary features of 148 known PID genes and 3162 non-PID genes as a training data set. The power of this approach is illustrated by the fact that six of the predicted genes have recently been experimentally confirmed to be PID genes. The remaining genes in this predicted data set represent attractive candidates for testing in patients where the etiology cannot be ascribed to any of the known PID genes.

Primary Structural Documentation of the Major Urinary Protein of the Indian Commensal Rat (Rattus rattus) Using a Proteomics Platform

Purpose: A number of proteome studies have been performed recently to identify pheromone-related protein expression and their molecular function using genetically modified rodents' urine. However, no such studies have used Indian commensal rodents; interestingly, in a previous investigation, we confirmed the presence of volatile molecules in commensal rodents urine and these molecules seem to be actively involved in pheromonal communication. Therefore, the present study aims to identify the major urinary protein [MUP] present in commensal rat urine, which will help us to understand the protein's expression pattern and intrinsic properties among the rodents globally. Experimental Design: Initially, the total urinary proteins were separated by 1-D and 2-D electrophoresis and then subsequently analyzed by Matrix Assisted Laser Desorption Ionization-Time of Flight and Mass Spectrometer (MALDI-TOF/MS). Furthermore, they were then fragmented with the aid of a Tandem Mass Spectrometer (TOF/TOF) and the identified sequences aligned and confirmed using similarity with the deduced primary structures of members of the lipocalin superfamily.Results: The SDS-PAGE protein profiles showed distinct proteins with molecular masses of 15, 22.4, 25, 28, 42, 50, 55, 68, and 91 kDa. Of these proteins, the 22.4 kDa protein was considered to be target candidate. When 2D electrophoresis was carried out, about similar to 50 spots were detected with different masses and various pI ranges. The 22.4 kDa protein was found to have a pI of about 4.9. This 22.4 kDa protein spot was digested and subjected to mass spectrometry; it was identified as rat MUP. The fragmented peptides from the rat MUP at 935, 1026, 1192, and 1303 m/z were further fragmented with the aid of MS/MS and generated de novo sequence and this confirmed this protein to be the MUP present in the urine of commensal rats.Conclusion: The present investigation confirms the presence of MUP with a molecular mass of 22.4 kDa in the urine of commensal rats. This protein may be involved in the binding of volatile molecules and opens up a discussion about how volatile and non-volatile molecules in the commensal rats' urine may contribute chemo-communication.

Plasticity in the primary binding site of galactose/N-acetylgalactosamine-specific lectins - Implication of the C-H center dot center dot center dot O hydrogen bond at the specificity-determining C-4 locus of the saccharide in 4-methoxygalactose recognition by jacalin and winged bean (basic) agglutinin I

It is currently believed that an unsubstituted axial hydroxyl at the specificity-determining C-4 locus of galactose is indispensable for recognition by galactose/N-acetylgalactosamine-specific lectins. Titration calorimetry demonstrates that 4-methoxygalactose retains binding allegiance to the Moraceae lectin jacalin and the Leguminosae lectin, winged bean (basic) agglutinin (WBA I). The binding reactions were driven by dominant favorable enthalpic contributions and exhibited significant enthalpy-entropy compensation. Proton NMR titration of C-methoxygalactose with jacalin and WBA I resulted in broadening of the sugar resonances without any change in chemical shift. The alpha-and beta-anomers of 4-methoxygalactose were found to be in slow exchange with free and lectin-bound states. Both the anomers experience magnetically equivalent environments at the respective binding sites. The binding constants derived from the dependence of NMR line widths on 4-methoxygalactose concentration agreed well with those obtained from titration calorimetry. The results unequivocally demonstrate that the loci corresponding to the axially oriented C-4 hydroxyl group of galactose within the primary binding site of these lectins exhibit plasticity. These analyses suggest, for the first time, the existence of C-H ... O-type hydrogen-bond(s) in protein-carbohydrate interactions in general and between the C-4 locus of galactose derivative and the lectins jacalin and WBA I in particular.