PenguinMan/ARXIV huggingface dataset:
The vigor paradox: saccade velocity during deliberation encodes utility of effortful actions
Authors: Korbisch, C. C.; Apuan, D.; Shadmehr, R.; Ahmed, A. A.
Category: neuroscience
Date: 2022-03-11

During deliberation, as the brain considers its options, the neural activity representing the goodness of each option rises toward a threshold, and the choice is often dictated by the option for which the rise is fastest. Here we report a surprising correlate of these activities: saccade vigor. We engaged human subjects in a decision-making task in which they considered effortful options, each requiring walking various durations and inclines. As they deliberated, they made saccades between the symbolic representations of those options. These saccades had no bearing on the effort that they would later expend, yet as they deliberated, saccade velocities increased. The rate of rise in vigor was faster for saccades toward the option that they later indicated as their choice, and encoded the difference in the subjective value of the two effortful options. Once deliberation ended, following a brief delay the subjects indicated their choice by making another saccade. Remarkably, vigor for this saccade dropped to baseline and no longer encoded subjective value. These results are consistent with an urgency model of decision-making in which a global signal in the brain drives both the neural circuits that make decisions, and the neural circuits that make movements. Paradoxically, this common drive is shared between the oculomotor circuits and the decision-making circuits, even when the decision involves effortful expenditure during a future event.

SignificanceThere is a link between the decisions we make and the movements that follow. Not only do we prefer options of greater value, but we also move faster to acquire them. When deliberating between options, neural activity rises to a threshold and the option that wins this race is the one chosen. We report a potential correlate of this in the motor control circuits; during deliberation, saccade vigor to both options rise, but faster for the option ultimately chosen. Thus, our movements appear to mirror the neural activity conducting the decision-making process. Paradoxically, this is true even when the movements have no direct bearing on the decision at hand.
NLGC: Network Localized Granger Causality with Application to MEG Directional Functional Connectivity Analysis
Authors: Soleimani, B.; Das, P.; Karunathilake, I. M. D.; Kuchinsky, S. E.; Simon, J. Z.; Babadi, B.
Category: neuroscience
Date: 2022-05-11

Identifying the directed connectivity that underlie networked activity between different cortical areas is critical for understanding the neural mechanisms behind sensory processing. Granger causality (GC) is widely used for this purpose in functional magnetic resonance imaging analysis, but there the temporal resolution is low, making it difficult to capture the millisecond-scale interactions underlying sensory processing. Magnetoencephalography (MEG) has millisecond resolution, but only provides low-dimensional sensor-level linear mixtures of neural sources, which makes GC inference challenging. Conventional methods proceed in two stages: First, cortical sources are estimated from MEG using a source localization technique, followed by GC inference among the estimated sources. However, the spatiotemporal biases in estimating sources propagate into the subsequent GC analysis stage, may result in both false alarms and missing true GC links. Here, we introduce the Network Localized Granger Causality (NLGC) inference paradigm, which models the source dynamics as latent sparse multivariate autoregressive processes and estimates their parameters directly from the MEG measurements, integrated with source localization, and employs the resulting parameter estimates to produce a precise statistical characterization of the detected GC links. We offer several theoretical and algorithmic innovations within NLGC and further examine its utility via comprehensive simulations and application to MEG data from an auditory task involving tone processing from both younger and older participants. Our simulation studies reveal that NLGC is markedly robust with respect to model mismatch, network size, and low signal-to-noise ratio, whereas the conventional two-stage methods result in high false alarms and mis-detections. We also demonstrate the advantages of NLGC in revealing the cortical network-level characterization of neural activity during tone processing and resting state by delineating task- and age-related connectivity changes.
Odors drive feeding through gustatory receptor neurons in Drosophila
Authors: Wei, H.-p.; Lam, K. C.; Kazama, H.
Category: neuroscience
Date: 2022-03-12

Odors are intimately tied to the taste system to aid food selection and determine the sensory experience of food. However, how smell and taste are integrated in the nervous system to drive feeding behavior remains largely unknown. We show in Drosophila that odors alone activate gustatory receptor neurons (GRNs) and trigger proboscis extension reflex (PER), a canonical taste-evoked feeding behavior. Odor-evoked PER requires the function of sugar-sensing GRNs but not the olfactory organs. Calcium imaging shows that GRNs directly respond to odors. Odor-evoked PER is mediated by the Gr5a receptor, and is bidirectionally modulated by specific olfactory binding proteins. Finally, odors and sucrose co-applied to GRNs synergistically enhance PER. These results reveal a cell-intrinsic mechanism for odor-taste multimodal integration that takes place as early as in GRNs, indicating that unified chemosensory experience is a product of layered integration in peripheral neurons and in the brain.
BLMM: Parallelised Computing for Big Linear Mixed Models
Authors: Maullin-Sapey, T.; Nichols, T.
Category: neuroscience
Date: 2022-03-12

Within neuroimaging large-scale, shared datasets are becoming increasingly commonplace, challenging existing tools both in terms of overall scale and complexity of the study designs. As sample sizes grow, researchers are presented with new opportunities to detect and account for grouping factors and covariance structure present in large experimental designs. In particular, standard linear model methods cannot account for the covariance and grouping structures present in large datasets, and the existing linear mixed models (LMM) tools are neither scalable nor exploit the computational speed-ups afforded by vectorisation of computations over voxels. Further, nearly all existing tools for imaging (fixed or mixed effect) do not account for variability in the patterns of missing data near cortical boundaries and the edge of the brain, and instead omit any voxels with any missing data. Yet in the large-n setting, such a voxel-wise deletion missing data strategy leads to severe shrinkage of the final analysis mask. To counter these issues, we describe the "Big" Linear Mixed Models (BLMM) toolbox, an efficient Python package for large-scale fMRI LMM analyses. BLMM is designed for use on high performance computing clusters and utilizes a Fisher Scoring procedure made possible by derivations for the LMM Fisher information matrix and score vectors derived in our previous work, Maullin-Sapey and Nichols (2021).
Cellular stress in brain organoids is limited to a distinct and bioinformatically removable subpopulation.
Authors: Vertesy, A.; Eichmueller, O. L.; Naas, J.; Novatchkova, M.; Esk, C.; Balmana, M.; Ladstaetter, S.; Bock, C.; von Haeseler, A.; Knoblich, J. A.
Category: neuroscience
Date: 2022-03-12

Organoids enable disease modeling in complex and structured human tissue, in vitro. Like most 3D models, they lack sufficient oxygen supply, leading to cellular stress. These negative effects are particularly prominent in complex models, like brain organoids, where they can prevent proper lineage commitment. Here, we analyze brain organoid and fetal single cell RNA sequencing (scRNAseq) data from published and new datasets totaling over 190,000 cells. We describe a unique stress signature found in all organoid samples, but not in fetal samples. We demonstrate that cell stress is limited to a defined organoid cell population, and present Gruffi, an algorithm that uses granular functional filtering to identify and remove stressed cells from any organoid scRNAseq dataset in an unbiased manner. Our data show that adverse effects of cell stress can be corrected by bioinformatic analysis, improving developmental trajectories and resemblance to fetal data.
A Trem2*R47H mouse model without cryptic splicing drives age- and disease-dependent tissue damage and synaptic loss in response to plaques
Authors: Tran, K. M.; Kawauchi, S.; Rezaie, N.; Liang, H. Y.; Arreola, M.; Cunda, C. D.; Phan, J.; Collins, S.; Walker, A.; Milinkeviciute, G.; Neumann, J.; Gomez-Arboledas, A.; Javonillo, D. I.; Tran, K.; Gantuz, M.; Forner, S.; Swarup, V.; Tenner, A. J.; LeFerla, F.; Wood, M. A.; Mortazavi, A.; MacGregor, G.; Green, K. N.
Category: neuroscience
Date: 2022-03-12

Genome-Wide Association Studies revealed the TREM2 R47H variant as one of the strongest genetic risk factors for late-onset Alzheimers Disease (AD). Unfortunately, many current TREM2*R47H mouse models are associated with cryptic mRNA splicing of the mutant allele that produces a confounding reduction in protein product. We have developed the Trem2R47H NSS (Normal Splice Site) mouse model where the Trem2 allele is expressed at a similar level to the wild-type Trem2 allele, without evidence of cryptic splicing products, and appropriate inflammatory responses to cuprizone challenge. Utilizing the 5xFAD mouse model, we report age- and disease-dependent changes in response to pathology. At an early disease stage (4 mo), homozygous Trem2R47H NSS; hemizygous 5xFAD (Trem2R47H NSS ; 5xFAD) mice have reduced size and number of microglia plus impaired interaction with plaques, that is associated with increased dystrophic neurites and axonal damage detected through plasma neurofilament light chain (NfL) level and suppressed inflammation. However, homozygosity for Trem2R47H NSS suppressed LTP deficits and presynaptic puncta loss caused by the 5xFAD transgene array. At a more advanced disease stage (12 mo,) Trem2R47H NSS ; 5xFAD mice no longer display impaired plaque-microglia interaction or suppressed inflammatory gene expression, although NfL levels remain elevated, and a unique interferon-related gene expression signature is seen. Furthermore, Trem2R47H NSS ; 5xFAD mice also display robust LTP deficits and exacerbated presynaptic loss. Collectively, we provide a Trem2R47H variant mouse without cryptic splicing, and demonstrate it has disease stage dependent effects when combined with a plaque bearing model, with an initial loss of function that ultimately resolves, giving rise to a unique interferon signature and associated tissue damage.
Organization of thalamocortical structural covariance and a corresponding 3D atlas of the mouse thalamus
Authors: Yee, Y.; Ellegood, J.; French, L.; Lerch, J. P.
Category: neuroscience
Date: 2022-03-12

For information from sensory organs to be processed by the brain, it is usually passed to appropriate areas of the cerebral cortex. Almost all of this information passes through the thalamus, a relay structure that reciprocally connects to the vast majority of the cortex. The thalamus facilitates this information transfer through a set of thalamocortical connections that vary in cellular structure, molecular profiles, innervation patterns, and firing rates. Additionally, corticothalamic connections allow for intracortical information transfer through the thalamus. These efferent and afferent connections between the thalamus and cortex have been the focus of many studies, and the importance of cortical connectivity in defining thalamus anatomy is demonstrated by multiple studies that parcellate the thalamus based on cortical connectivity profiles.

Here, we examine correlated morphological variation between the thalamus and cortex, or thalamocortical structural covariance. For each voxel in the thalamus as a seed, we construct a cortical structural covariance map that represents correlated cortical volume variation, and examine whether high structural covariance is observed in cortical areas that are functionally relevant to the seed. Then, using these cortical structural covariance maps as features, we subdivide the thalamus into six non-overlapping regions (clusters of voxels), and assess whether cortical structural covariance is associated with cortical connectivity that specifically originates from these regions.

We show that cortical structural covariance is high in areas of the cortex that are functionally related to the seed voxel, cortical structural covariance varies along cortical depth, and sharp transitions in cortical structural covariance profiles are observed when varying seed locations in the thalamus. Subdividing the thalamus based on structural covariance, we additionally demonstrate that the six thalamic clusters of voxels stratify cortical structural covariance along the dorsal-ventral, medial-lateral, and anterior-posterior axes. These cluster-associated structural covariance patterns are prominently detected in cortical regions innervated by fibers projecting out of their related thalamic subdivisions.

Together, these results advance our understanding of how the thalamus and the cortex couple in their volumes. Our results indicate that these volume correlations reflect functional organization and structural connectivity, and further provides a novel segmentation of the mouse thalamus that can be used to examine thalamic structural variation and thalamocortical structural covariation in disease models.
Neuromodulation of striatal D1 cells shapes BOLD fluctuations in anatomically connected thalamic and cortical regions
Authors: Markicevic, M.; Sturman, O.; Bohacek, J.; Rudin, M.; Zerbi, V.; Fulcher, B. D.; Wenderoth, N.
Category: neuroscience
Date: 2022-03-12

Understanding how the brains macroscale dynamics are shaped by underlying microscale mechanisms is a key problem in neuroscience. In animal models, we can now investigate this relationship in unprecedented detail by directly manipulating cellular-level properties while measuring the whole-brain response using resting-state fMRI. Here we focused on understanding how blood-oxygen-level-dependent (BOLD) dynamics, measured within a structurally well-defined striato-thalamo-cortical circuit, are shaped by chemogenetically exciting or inhibiting D1 medium spiny neurons (MSNs) of the right dorsomedial striatum (CPdm). We characterize changes in both the BOLD dynamics of individual cortical and subcortical brain areas, and patterns of inter-regional coupling (functional connectivity) between pairs of areas. Using a classification approach based on a large and diverse set of time-series properties, we found that CPdm neuromodulation alters BOLD dynamics within thalamic subregions that project back to dorsomedial striatum. In the cortex, the strongest changes in local dynamics were observed in unimodal regions, i.e., regions that process information from a single sensory modality, while changes in the local dynamics weakened along a putative cortical hierarchical gradient towards transmodal regions. In contrast, a decrease in functional connectivity was observed only for cortico-striatal connections after D1 excitation. Our results provide a comprehensive understanding of how targeted cellular-level manipulations affect local BOLD dynamics at the macroscale, including the role of a circuits structural characteristics and hierarchical cortical level in shaping those dynamics. These findings contribute to ongoing attempts to understand the influence of structure-function relationships in shaping inter-regional communication at subcortical and cortical levels.
Development of a PET radioligand for α2δ-1 subunit of calcium channels for imaging neuropathic pain
Authors: Zhou, Y.-P.; Sun, Y.; Takahashi, K.; Belov, V.; Andrews, N.; Woolf, C. J.; Brugarolas, P.
Category: neuroscience
Date: 2022-03-12

Neuropathic pain affects 7-10% of the adult population. Being able to accurately monitor biological changes underlying neuropathic pain will improve our understanding of neuropathic pain mechanisms and facilitate the development of novel therapeutics. Positron emission tomography (PET) is a noninvasive molecular imaging technique that can provide quantitative information of biochemical changes at the whole-body level by using radiolabeled ligands. One important biological change underlying the development of neuropathic pain is the overexpression of 2{delta}-1 subunit of voltage-dependent calcium channels (the target of gabapentin). Thus, we hypothesized that a radiolabeled form of gabapentin may allow imaging changes in 2{delta}-1 for monitoring the underlying pathophysiology of neuropathic pain. Here, we report the development of two 18F-labeled derivatives of gabapentin (trans-4-[18F]fluorogabapentin and cis-4-[18F]fluorogabapentin) and their evaluation in healthy rats and a rat model of neuropathic pain (spinal nerve ligation model). Both isomers were found to selectively bind to the 2{delta}-1 receptor with trans-4-[18F]fluorogabapentin having a higher affinity. Both tracers displayed around 1.5- to 2-fold increased uptake in injured nerves over the contralateral uninjured nerves when measured by gamma counting ex vivo. Although the small size of the nerves and the signal from surrounding muscle prevented visualizing these changes using PET, this work demonstrates that fluorinated derivatives of gabapentin retain binding to 2{delta}-1 and that their radiolabeled forms can be used to detect pathological changes in vitro and ex vivo. Furthermore, this work confirms that 2{delta}-1 is a promising target for imaging specific features of neuropathic pain.
Electroencephalographic dynamics of rhythmicbreath-based meditation
Authors: Tripathi, V.; Bhasker, L.; Kharya, C.; Bhatia, M.; Kochupillai, V.
Category: neuroscience
Date: 2022-03-24

Meditation has been practised for millennia but the neuroscientific understanding of the dynamics is still lacking. Sudarshan Kriya Yoga (SKY) is an evidence based breathing based meditation technique that utilizes rhythmic breathing to induce a deep state of relaxation and calm. Multiple studies have found benefits of the SKY technique from genetic, physiological, psychological to behaviour levels. We collected Electroencephalographic (EEG) data in 43 subjects who underwent the SKY technique and analysed the brain rhythms at different stages of the technique namely preparatory breathing (Pranayama), rhythmic breathing (Kriya) and meditation (Yoga Nidra) using newly developed methods to analyse periodic and aperiodic components. Alpha waves amplitude in the parieto-occipital region decreased as the rhythmic breathing progressed and dropped sharply during the meditation period. Theta amplitudes and peak frequency increased in the centro-frontal region during the rhythmic breathing period but were marked by sustained low theta waves during the meditation period. The delta wave amplitude was not affected by breathing but both delta band power and peak frequency increased during the meditation period in the centro-frontal region. We also saw a decrease in the 1/f aperiodic signal across the brain during the meditation period suggesting a modification of excitation-inhibition balance. We see an overall slowing down of brain oscillations from alpha to theta to delta as the meditation progressed. The paper studies in depth the transitional dynamics of the SKY technique analysing the alpha, theta, delta waves and aperiodic signals and demonstrates that each phase in a breathing based meditation has a unique electrophysiological signature.
The role of UBE3A in the autism and epilepsy-related Dup15q syndrome using patient-derived, CRISPR-corrected neurons
Authors: Elamin, M.; Dumarchey, A.; Stoddard, C.; Robinson, T. M.; Cowie, C.; Gorka, D.; Chamberlain, S. J.; Levine, E. S.
Category: neuroscience
Date: 2022-03-12

Chromosome 15q11-q13 duplication syndrome (Dup15q) is a neurodevelopmental disorder caused by maternal duplications of this region. Autism and epilepsy are key features of Dup15q, but affected individuals also exhibit intellectual disability and developmental delay. UBE3A, the gene encoding the ubiquitin protein ligase E3A, is likely a major driver of Dup15q because individuals with maternal, but not paternal 15q duplications have the disorder, and UBE3A is the only imprinted gene expressed solely from the maternal allele in mature neurons. Nevertheless, the exact role of UBE3A is yet to be determined. To establish whether UBE3A overexpression is required for Dup15q neuronal deficits, UBE3A expression was manipulated in patient-derived induced pluripotent stem cell (iPSC) lines using antisense oligonucleotides. Dup15q neurons exhibited hyperexcitability features compared to genome-edited isogenic control neurons, and this phenotype was generally prevented by normalizing UBE3A levels throughout in vitro development. Overexpression of UBE3A in an iPSC line with a paternal duplication resulted in a profile similar to that of Dup15q neurons except for synaptic phenotypes. These results indicate that UBE3A overexpression is necessary for most Dup15q cellular phenotypes. However, the inability of UBE3A overexpression to recapitulate synaptic phenotypes suggests an important role for non-imprinted genes in the disorder.