PenguinMan/ARXIV huggingface dataset:
A comparison of bivariate frequency domain measures of electrophysiological connectivity
Authors: Pascual-Marqui, R. D.; Faber, P.; Kinoshita, T.; Kochi, K.; Milz, P.; Nishida, K.; Yoshimura, M.
Category: neuroscience
Date: 2018-11-19

1The problem of interest here concerns electrophysiological signals from two cortical sites, acquired as invasive intracranial recordings, or from non-invasive estimates of cortical electric neuronal activity computed from EEG or MEG recordings (see e.g. https://doi.org/10.1101/269753). In the absence of other sources, these measured signals consist of an instantaneous linear mixture of the true, actual, unobserved local signals, due to low spatial resolution and volume conduction. A connectivity measure is unreliable as a true indicator of electrophysiological connectivity if it is not invariant to mixing, or if it reports a significant connection for a mixture of independent signals. In (Vinck et al 2011 Neuroimage 55:1548) it was shown that coherence, imaginary coherence, and phase locking value are not invariant to mixing, while the phase lag index (PLI) and the weighted version (wPLI) are invariant to mixing. Here we show that the lagged coherence (LagCoh) measure (2007, https://arxiv.org/abs/0711.1455), not studied in Vinck et al, is invariant to mixing. Additionally, we present here a new mixture-invariant connectivity statistic: the \"standardized imaginary covariance\" (sImCov). We also include in the comparisons the directed PLI (dPLI) by Stam et al (2012 Neuroimage 62:1415). Fourier coefficients for \"N\" trials are generated from a linear unidirectional causal time domain model with electrophysiological delay \"k\" and regression coefficient \"b\". 1000 random data sets of \"N\" trials are simulated, and for each one, and for each connectivity measure, non-parametric randomization tests are performed. The \"true positive detection rate\" is calculated as the fraction of 1000 cases that have significant connectivity at p<0.05, 0.1, and 0.2. The connectivity methods were compared in terms of detection rates, under non-mixed and mixed conditions, for small and large sample sizes \"N\", with and without jitter for \"k\", and for different values of signal to noise. Under mixing, the results show that LagCoh outperforms wPLI, PLI, dPLI, and sImCov. Without mixing, LagCoh and sImCov outperform wPLI, PLI, and dPLI. Finally, it is shown that dPLI is an invalid estimator of flow direction, i.e. it reverses and \"goes against the flow\" by merely changing the sign of one of the time series, a fact that violates the basic definition of Granger causality. For the sake of reproducible research, the supplementary material includes Delphi Pascal source code and all detailed result files in human readable format.
One Model to Rule Them All: An Integrative Approach to Matrix-Based Analyses in Neuroimaging Connectomics
Authors: Chen, G.; Bürkner, P.-C.; Taylor, P. A.; Li, Z.; Yin, L.; Glen, D. R.; Kinnison, J. K.; Cox, R. W.; Pessoa, L.
Category: neuroscience
Date: 2018-11-01

Network modeling in neuroimaging holds promise in probing the interrelationships among brain regions and potential clinical applications. Two types of matrix-based analysis (MBA) are usually seen in neuroimaging connectomics: one is the functional attribute matrix (FAM) of, for example, correlations, that measures the similarity of BOLD response patterns among a list of predefined regions of interest (ROIs). Another type of MBA involves the structural attribute matrix (SAM), e.g., describing the properties of white matter between any pair of gray-matter regions such as fractional anisotropy, mean diffusivity, axial and radial diffusivity. There are different methods that have been developed or adopted to summarize such matrices across subjects, including general linear models (GLMs) and various versions of graph theoretic analysis. We argue that these types of modeling strategies tend to be \"inefficient\" in statistical inferences and have many pitfalls, such as having strong dependence on arbitrary thresholding under conventional statistical frameworks.\n\nHere we offer an alternative approach that integrates the analyses of all the regions, region pairs (RPs) and subjects into one framework, called Bayesian multilevel (BML) modeling. In this approach, the intricate relationships across regions as well as across RPs are quantitatively characterized. This integrative approach avoids the multiple testing issue that typically plagues the conventional statistical analysis in neuroimaging, and it provides a principled way to quantify both the effect and its uncertainty at each region as well as for each RP. As a result, a unique feature of BML is that the effect at each region and the corresponding uncertainty can be estimated, revealing the relative strength or importance of each region; in addition, the effect at each RP is obtained along with its uncertainty as statistical evidence. Most importantly, the BML approach can be scrutinized for consistency through validation and comparisons with alternative assumptions or models. We demonstrate the BML methodology with a real dataset with 16 ROIs from 41 subjects, and compare it to the conventional GLM approach in terms of model efficiency, performance and inferences. Furthermore, we emphasize the notion of full results reporting through \"highlighting,\" instead of through the common practice of \"hiding.\" The associated program will be available as part of the AFNI suite for general use.
Studying Human Behaviour with Virtual Reality: The Unity Experiment Framework
Authors: Brookes, J.; Warburton, M.; Alghadier, M.; Mushtaq, F.; Mon-Williams, M.
Category: neuroscience
Date: 2018-11-01

Virtual Reality systems offer a powerful tool for human behaviour research. The ability to create three-dimensional visual scenes and measure responses to the visual stimuli enables the behavioural researcher to test hypotheses in a manner and scale that were previously unfeasible. For example, a researcher wanting to understand interceptive timing behaviour might wish to violate Newtonian mechanics, so objects move in novel 3D trajectories. The same researcher may wish to collect such data with hundreds of participants outside the laboratory, and the use of a VR headset makes this a realistic proposition. The difficulty facing the researcher is that sophisticated 3D graphics engines (e.g. Unity) have been created for game designers rather than behavioural scientists. In order to overcome this barrier, we have created a set of tools and programming syntaxes that allow logical encoding of the common experimental features required by the behavioural scientist. The Unity Experiment Framework (UXF) allows the researcher to readily implement several forms of data collection, and provides researchers with the ability to easily modify independent variables. UXF does not offer any stimulus presentation features, so the full power of the Unity game engine can be exploited. We use a case study experiment, measuring postural sway in response to an oscillating virtual room, to show how UXF can replicate and advance upon behavioural research paradigms. We show that UXF can simplify and speed up development of VR experiments created in commercial gaming software and facilitate the efficient acquisition of large quantities of behavioural research data.
High interictal connectivity within the resection zone is associated with favorable post-surgical outcomes in focal epilepsy patients
Authors: Shah, P.; Bernabei, J.; Kini, L.; Ashourvan, A.; Boccanfuso, J.; Archer, R.; Oechsel, K.; Lucas, T. H.; Bassett, D. S.; Davis, K.; Litt, B.
Category: neuroscience
Date: 2018-11-01

ObjectivePatients with drug-resistant focal epilepsy are often candidates for invasive surgical therapies. In these patients, it is necessary to accurately localize seizure-generators to ensure seizure freedom following intervention. While intracranial electroencephalography (iEEG) is the gold standard for mapping networks for surgery, this approach requires inducing and recording seizures, which may cause patient morbidity. The goal of this study is to evaluate the utility of mapping interictal (non-seizure) iEEG networks to identify targets for surgical treatment.\n\nMethodsWe analyze interictal iEEG recordings and neuroimaging from 27 focal epilepsy patients treated via surgical resection. We generate interictal functional networks by calculating pairwise correlation of iEEG signals across different frequency bands. We identify electrodes falling within surgically resected tissue (i.e. the resection zone), and compute node-level and edge-level synchrony in relation to the resection zone. We associate these metrics with post-surgical outcomes.\n\nResultsGreater overlap between resected electrodes and highly synchronous electrodes is associated with favorable post-surgical outcomes. Additionally, good outcome patients have significantly higher connectivity localized within the resection zone compared to those with poorer postoperative seizure control. This finding persists following normalization by a spatially-constrained null model.\n\nConclusionsThis study suggests that spatially-informed interictal network synchrony measures can distinguish between good and poor post-surgical outcomes. By capturing clinically relevant information during interictal periods, our method may ultimately reduce the need for prolonged invasive implants and provide insights into the pathophysiology of an epileptic brain. We discuss next steps for translating these findings into a prospectively useful clinical tool.
Applying non-parametric testing to discrete transfer entropy
Authors: Olin-Ammentorp, W.; Cady, N.
Category: neuroscience
Date: 2018-11-02

Transfer entropy (TE) is a powerful algorithm which attempts to detect the transfer of information from one system to another. In neuroscience, it has the potential to track the movement of information through complex neuronal systems, and provide powerful insights into their organization and operation. One such application is the ability to infer the existence of causal connectivity (such as synaptic pathways) between neurons in a culture being recorded by micro-electrode array (MEA).\n\nThere are several challenges, however, in applying TE to neurological data; one of these is the ability to robustly classify what experimental TE value qualifies as significant. We find that common methods in spike train analysis such as a Z-test cannot be applied, as their assumptions are not met. Instead, we utilize surrogate data to compute a sample under the null hypothesis (no causal connection), and resample experimental data through Markov chain Monte Carlo (MCMC) methods to create a sample of TE values under experimental conditions. A standard non-parametric test (Mann-Whitney U-test) is then applied to compare these samples, and determine if they represent a significant connection.\n\nWe have applied this methodology to MEA recordings of neuronal cultures developing over a period of roughly a month, and find that it provides a wealth of information regarding the cultures' maturity. This includes features such as the directed graph of causal connections across the MEA and identification of information exchange centers. These results are consistent and carry a well-defined significance level.
Intermittent abstinence potentiates oxycodone self-administration under extended access conditions and alters brain reward sensitivity
Authors: Nguyen, J.; Grant, Y.; Taffe, M.
Category: neuroscience
Date: 2018-11-02

Background and PurposeThe extra-medical use of, and addiction to, prescription opioid analgesics is a growing health problem. To characterize how prescription opioid abuse develops, this study investigated the affective consequences of escalating prescription opioid use using intracranial self-stimulation (ICSS) reward and oxycodone intravenous self-administration (IVSA) models.

Experimental ApproachMale Wistar rats were given access to oxycodone IVSA (0.15 mg/kg/infusion, i.v.) in Short Access (ShA; 1 h) or Long Access (LgA; 12 h) sessions for 5 sessions/week followed by intermittent 60 h discontinuations from drug access, a novel explicit test of the negative reinforcement hypothesis. Separate groups were first trained in the ICSS procedure and then in oxycodone IVSA in 11 h LgA sessions.

Key ResultsRats given LgA to oxycodone escalated their responding more than ShA rats, with further significant increases observed following each 60 h discontinuation. Pre-session brain reward thresholds increased with sequential daily LgA IVSA sessions, consistent with a growing negative affective state consequent to successive daily intoxication/abstinence cycles. A 1 h oxycodone IVSA interval was sufficient to normalize these elevated reward thresholds, as was, paradoxically, a 60 h weekend abstinence. The increase in ICSS thresholds was attenuated in a group treated with the long-acting kappa opioid antagonist norBNI prior to IVSA training.

Conclusions and ImplicationsChanges in brain reward function during escalation of oxycodone self-administration are driven by an interplay between kappa opioid receptor-mediated negative affective state associated with escalated oxycodone intake and dynamic restoration of brain reward status during longer periods of abstinence.
Intracellular calcium stores mediate a synapse-specific form of metaplasticity at hippocampal dendritic spines
Authors: Mahajan, G.; Nadkarni, S.
Category: neuroscience
Date: 2018-11-02

Long-term plasticity mediated by NMDA receptors supports input-specific, Hebbian forms of learning at excitatory CA3-CA1 connections in the hippocampus. An additional layer of stabilizing mechanisms that act globally as well as locally over multiple time scales may be in place to ensure that plasticity occurs in a constrained manner. Here, we investigate the potential role of calcium (Ca2+) stores associated with the endoplasmic reticulum (ER) in the local regulation of plasticity dynamics at individual CA1 synapses. Our study is spurred by (1) the curious observation that ER is sparsely distributed in dendritic spines, but over-represented in large spines that are likely to have undergone activity-dependent strengthening, and (2) evidence suggesting that ER motility within synapses can be rapid, and accompany activity-regulated spine remodeling. Based on a physiologically realistic computational model for ER-bearing CA1 spines, we characterize the contribution of IP3-sensitive Ca2+ stores to spine Ca2+ dynamics during activity patterns mimicking the induction of long-term potentiation (LTP) and depression (LTD). Our results suggest graded modulation of the NMDA receptor-dependent plasticity profile by ER, which selectively enhances LTD induction. We propose that spine ER can locally tune Ca2+-based plasticity on an as-needed basis, providing a braking mechanism to mitigate runaway strengthening at potentiated synapses. Our model suggests that the presence of ER in the CA1 spine may promote re-use of synapses with saturated strengths.
Optoacoustic brains stimulation at submillimeter spatial precision
Authors: Jiang, Y.; Lee, H. J.; Lan, L.; Tseng, H.-a.; Yang, C.; Man, H.-Y.; Han, X.; Cheng, J.-X.
Category: neuroscience
Date: 2018-11-02

Low-intensity ultrasound is an emerging modality for neuromodulation. Yet, piezo-based transducers offer poor spatial confinement of excitation volume, often bigger than a few millimeters in diameter. In addition, the bulky size limits their implementation in a wearable setting and prevents integration with other experimental modalities. Here, we report spatially confined optoacoustic neural stimulation through a novel miniaturized Fiber-Optoacoustic Converter (FOC). The FOC has a diameter of 600 m and generates omnidirectional ultrasound wave locally at the fiber tip through the optoacoustic effect. We show that the optoacoustic wave can directly activate individual cultured neurons and generate intracellular Ca2+ transients. The FOC activates neurons within a radius of 500 m around the fiber tip, delivering superior spatial resolution over conventional piezo-based low-frequency transducers. Combining FOC with electrophysiology, direct and spatially confined neural stimulation of mouse brain is achieved in vivo.
To CRUNCH or not to CRUNCH: Task Difficulty Affects Functional Brain Reorganisation during Visuospatial Working Memory Performance in Premanifest Huntington’s Disease
Authors: Soloveva, M. V.; Jamadar, S. D.; Velakoulis, D.; Poudel, G.; Georgiou-Karistianis, N.
Category: neuroscience
Date: 2018-11-02

Presymptomatic Huntingtons disease (pre-HD) individuals tend to increase functional brain activity to compensate for HD-related brain anomalies. We used a quantitative model of compensation, known as the CRUNCH (Compensation-Related Utilization of Neural Circuits Hypothesis) to explicitly characterise compensation in pre-HD. We acquired functionalmagnetic resonance imaging (fMRI) data (n = 15 pre-HD; n = 15 controls) during performance of an 18-minute fMRI visuospatial working memory task with low, intermediate-1, intermediate-2, and high memory loads. Consistent with the CRUNCH prediction, pre-HD individuals showed decreased fMRI activity in left intraparietal sulcus at high memory load, compared to healthy controls who showed increased fMRI activity in left intraparietal sulcus at high memory load. Contrary to the other CRUNCH prediction, the pre-HD group did not show compensatory increase in fMRI activity at lower levels of memory loads in left intraparietal sulcus. Our findings provide partial support for the validity of CRUNCH in pre-HD.\n\nHighlightsO_LIVisuospatial working memory deficits in pre-HD occur 25 years prior to predicted disease onset\nC_LIO_LITask demands differentially affect fMRI activity in left intraparietal sulcus\nC_LIO_LICRUNCH can partially apply in Huntingtons disease\nC_LI
Predicting time to dementia using a quantitative template of disease progression
Authors: Bilgel, M.; Jedynak, B. M.
Category: neuroscience
Date: 2018-11-02

IntroductionCharacterization of longitudinal trajectories of biomarkers implicated in sporadic Alzheimers disease (AD) in decades prior to clinical diagnosis is important for disease prevention and monitoring.\n\nMethodsWe used a multivariate Bayesian model to temporally align 1369 AD Neuroimaging Initiative participants based on the similarity of their longitudinal biomarker measures and estimated a quantitative template of the temporal evolution cerebrospinal fluid (CSF) A{beta}1-42, p-tau181p, and t-tau, hippocampal volume, brain glucose metabolism, and cognitive measurements. We computed biomarker trajectories as a function of time to AD dementia, and predicted AD dementia onset age in a disjoint sample.\n\nResultsQuantitative template showed early changes in verbal memory, CSF A{beta}1-42 and p-tau181p, and hippocampal volume. Mean error in predicted AD dementia onset age was < 1.5 years.\n\nDiscussionOur method provides a quantitative approach for characterizing the natural history of AD starting at preclinical stages despite the lack of individual-level longitudinal data spanning the entire disease timeline.
Biomarker-based patient selection improves stroke rehabilitation trial efficiency
Authors: Stinear, C. M.; Byblow, W. D.; Barber, P. A.; Ackerley, S.; Smith, M.-C.; Cramer, S.
Category: neuroscience
Date: 2018-11-02

Background and PurposeInter-subject variability complicates trials of novel stroke rehabilitation therapies, particularly in the sub-acute phase after stroke. We tested whether selecting patients using motor evoked potential (MEP) status, a physiological biomarker of motor system function, could improve trial efficiency.\n\nMethodsA retrospective analysis of data from 207 patients (103 women, mean (SD) 70.6 (15.1) years) was used to estimate sample sizes and recruitment rates required to detect a 7-point difference between hypothetical control and treatment groups in upper-limb Fugl-Meyer and Action Research Arm Test scores at 90 days post-stroke. Analyses were carried out for the full sample and for subsets defined by motor evoked potential (MEP) status.\n\nResultsSelecting patients according to MEP status reduced the required sample size by 75% compared to an unselected sample. The estimated time needed to recruit the required sample was also reduced by 72% for patients with MEPs, and was increased by 2-3-fold for patients without MEPs.\n\nConclusionsUsing biomarkers to select patients can improve stroke rehabilitation trial efficiency by reducing the sample size and recruitment time needed to detect a clinically meaningful effect of the tested intervention.
Integrative inference of brain cell similarities and differences from single-cell genomics
Authors: Welch, J.; Kozareva, V.; Ferreira, A.; Vanderburg, C.; Martin, C.; Macosko, E.
Category: neuroscience
Date: 2018-11-02

Defining cell types requires integrating diverse measurements from multiple experiments and biological contexts. Recent technological developments in single-cell analysis have enabled high-throughput profiling of gene expression, epigenetic regulation, and spatial relationships amongst cells in complex tissues, but computational approaches that deliver a sensitive and specific joint analysis of these datasets are lacking. We developed LIGER, an algorithm that delineates shared and dataset-specific features of cell identity, allowing flexible modeling of highly heterogeneous single-cell datasets. We demonstrated its broad utility by applying it to four diverse and challenging analyses of human and mouse brain cells. First, we defined both cell-type-specific and sexually dimorphic gene expression in the mouse bed nucleus of the stria terminalis, an anatomically complex brain region that plays important roles in sex-specific behaviors. Second, we analyzed gene expression in the substantia nigra of seven postmortem human subjects, comparing cell states in specific donors, and relating cell types to those in the mouse. Third, we jointly leveraged in situ gene expression and scRNA-seq data to spatially locate fine subtypes of cells present in the mouse frontal cortex. Finally, we integrated mouse cortical scRNA-seq profiles with single-cell DNA methylation signatures, revealing mechanisms of cell-type-specific gene regulation. Integrative analyses using the LIGER algorithm promise to accelerate single-cell investigations of cell-type definition, gene regulation, and disease states.