Neuroinformatics最新文献

Large-scale diffusion MRI tractography remains a significant challenge. Users must orchestrate a complex sequence of instructions that requires many software packages with complex dependencies and high computational costs. We developed MaPPeRTrac, an edge-centric tractography pipeline that simplifies and accelerates this process in a wide range of high-performance computing (HPC) environments. It fully automates either probabilistic or deterministic tractography, starting from a subject's magnetic resonance imaging (MRI) data, including structural and diffusion MRI images, to the edge density image (EDI) of their structural connectomes. Dependencies are containerized with Singularity (now called Apptainer) and decoupled from code to enable rapid prototyping and modification. Data derivatives are organized with the Brain Imaging Data Structure (BIDS) to ensure that they are findable, accessible, interoperable, and reusable following FAIR principles. The pipeline takes full advantage of HPC resources using the Parsl parallel programming framework, resulting in the creation of connectome datasets of unprecedented size. MaPPeRTrac is publicly available and tested on commercial and scientific hardware, so it can accelerate brain connectome research for a broader user community. MaPPeRTrac is available at: https://github.com/LLNL/mappertrac .

T1-weighted (T1w) MRI has low frequency intensity artifacts due to magnetic field inhomogeneities. Removal of these biases in T1w MRI images is a critical preprocessing step to ensure spatially consistent image interpretation. N4ITK bias field correction, the current state-of-the-art, is implemented in such a way that makes it difficult to port between different pipelines and workflows, thus making it hard to reimplement and reproduce results across local, cloud, and edge platforms. Moreover, N4ITK is opaque to optimization before and after its application, meaning that methodological development must work around the inhomogeneity correction step. Given the importance of bias fields correction in structural preprocessing and flexible implementation, we pursue a deep learning approximation / reinterpretation of the N4ITK bias fields correction to create a method which is portable, flexible, and fully differentiable. In this paper, we trained a deep learning network "DeepN4" on eight independent cohorts from 72 different scanners and age ranges with N4ITK-corrected T1w MRI and bias field for supervision in log space. We found that we can closely approximate N4ITK bias fields correction with naïve networks. We evaluate the peak signal to noise ratio (PSNR) in test dataset against the N4ITK corrected images. The median PSNR of corrected images between N4ITK and DeepN4 was 47.96 dB. In addition, we assess the DeepN4 model on eight additional external datasets and show the generalizability of the approach. This study establishes that incompatible N4ITK preprocessing steps can be closely approximated by naïve deep neural networks, facilitating more flexibility. All code and models are released at https://github.com/MASILab/DeepN4 .

Visibility graphs provide a novel approach for analysing time-series data. Graph theoretical analysis of visibility graphs can provide new features for data mining applications in fMRI. However, visibility graphs features have not been used widely in the field of neuroscience. This is likely due to a lack of understanding of their robustness in the presence of noise (e.g., motion) and their test-retest reliability. In this study, we investigated visibility graph properties of fMRI data in the human connectome project (N = 1010) and tested their sensitivity to motion and test-retest reliability. We also characterised the strength of connectivity obtained using degree synchrony of visibility graphs. We found that strong correlation (r > 0.5) between visibility graph properties, such as the number of communities and average degrees, and motion in the fMRI data. The test-retest reliability (Intraclass correlation coefficient (ICC)) of graph theoretical features was high for the average degrees (0.74, 95% CI = [0.73, 0.75]), and moderate for clustering coefficient (0.43, 95% CI = [0.41, 0.44]) and average path length (0.41, 95% CI = [0.38, 0.44]). Functional connectivity between brain regions was measured by correlating the visibility graph degrees. However, the strength of correlation was found to be moderate to low (r < 0.35). These findings suggest that even small movement in fMRI data can strongly influence robustness and reliability of visibility graph features, thus, requiring robust motion correction strategies prior to data analysis. Further studies are necessary for better understanding of the potential application of visibility graph features in fMRI.

Magnetic resonance spectroscopy (MRS) is widely used to estimate concentrations of glutamate and (gamma)-aminobutyric acid (GABA) in specific regions of the living human brain. As cytoarchitectural properties differ across the brain, interpreting these measurements can be assisted by having knowledge of such properties for the MRS region(s) studied. In particular, some knowledge of likely local neurotransmitter receptor patterns can potentially give insights into the mechanistic environment GABA- and glutamatergic neurons are functioning in. This may be of particular utility when comparing two or more regions, given that the receptor populations may differ substantially across them. At the same time, when studying MRS data from multiple participants or timepoints, the homogeneity of the sample becomes relevant, as measurements taken from areas with different cytoarchitecture may be difficult to compare. To provide insights into the likely cytoarchitectural environment of user-defined regions-of-interest, we produced an easy to use tool - InSpectro-Gadget - that interfaces with receptor mRNA expression information from the Allen Human Brain Atlas. This Python tool allows users to input masks and automatically obtain a graphical overview of the receptor population likely to be found within. This includes comparison between multiple masks or participants where relevant. The receptors and receptor subunit genes featured include GABA- and glutamatergic classes, along with a wide range of neuromodulators. The functionality of the tool is explained here and its use is demonstrated through a set of example analyses. The tool is available at https://github.com/lizmcmanus/Inspectro-Gadget.

The international neuroscience community is building the first comprehensive atlases of brain cell types to understand how the brain functions from a higher resolution, and more integrated perspective than ever before. In order to build these atlases, subsets of neurons (e.g. serotonergic neurons, prefrontal cortical neurons etc.) are traced in individual brain samples by placing points along dendrites and axons. Then, the traces are mapped to common coordinate systems by transforming the positions of their points, which neglects how the transformation bends the line segments in between. In this work, we apply the theory of jets to describe how to preserve derivatives of neuron traces up to any order. We provide a framework to compute possible error introduced by standard mapping methods, which involves the Jacobian of the mapping transformation. We show how our first order method improves mapping accuracy in both simulated and real neuron traces under random diffeomorphisms. Our method is freely available in our open-source Python package brainlit.

Decisions made during the analysis or reporting of an fMRI study influence the eligibility of that study to be entered into a meta-analysis. In a meta-analysis, results of different studies on the same topic are combined. To combine the results, it is necessary that all studies provide equivalent pieces of information. However, in task-based fMRI studies we see a large variety in reporting styles. Several specific meta-analysis methods have been developed to deal with the reporting practices occurring in task-based fMRI studies, therefore each requiring a specific type of input. In this manuscript we provide an overview of the meta-analysis methods and the specific input they require. Subsequently we discuss how decisions made during the study influence the eligibility of a study for a meta-analysis and finally we formulate some recommendations about how to report an fMRI study so that it complies with as many meta-analysis methods as possible.

Current mesoscale connectivity atlases provide limited information about the organization of thalamocortical projections in the mouse brain. Labeling the projections of spatially restricted neuron populations in thalamus can provide a functionally relevant level of connectomic analysis, but these need to be integrated within the same common reference space. Here, we present a pipeline for the segmentation, registration, integration and analysis of multiple tract-tracing experiments. The key difference with other workflows is that the data is transformed to fit the reference template. As a test-case, we investigated the axonal projections and intranuclear arrangement of seven neuronal populations of the ventral posteromedial nucleus of the thalamus (VPM), which we labeled with an anterograde tracer. Their soma positions corresponded, from dorsal to ventral, to cortical representations of the whiskers, nose and mouth. They strongly targeted layer 4, with the majority exclusively targeting one cortical area and the ones in ventrolateral VPM branching to multiple somatosensory areas. We found that our experiments were more topographically precise than similar experiments from the Allen Institute and projections to the primary somatosensory area were in agreement with single-neuron morphological reconstructions from publicly available databases. This pilot study sets the basis for a shared virtual connectivity atlas that could be enriched with additional data for studying the topographical organization of different thalamic nuclei. The pipeline is accessible with only minimal programming skills via a Jupyter Notebook, and offers multiple visualization tools such as cortical flatmaps, subcortical plots and 3D renderings and can be used with custom anatomical delineations.

To simulate whole brain dynamics with only a few equations, biophysical, mesoscopic models of local neuron populations can be connected using empirical tractography data. The development of mesoscopic mean-field models of neural populations, in particular, the Adaptive Exponential (AdEx mean-field model), has successfully summarized neuron-scale phenomena leading to the emergence of global brain dynamics associated with conscious (asynchronous and rapid dynamics) and unconscious (synchronized slow-waves, with Up-and-Down state dynamics) brain states, based on biophysical mechanisms operating at cellular scales (e.g. neuromodulatory regulation of spike-frequency adaptation during sleep-wake cycles or anesthetics). Using the Virtual Brain (TVB) environment to connect mean-field AdEx models, we have previously simulated the general properties of brain states, playing on spike-frequency adaptation, but have not yet performed detailed analyses of other parameters possibly also regulating transitions in brain-scale dynamics between different brain states. We performed a dense grid parameter exploration of the TVB-AdEx model, making use of High Performance Computing. We report a remarkable robustness of the effect of adaptation to induce synchronized slow-wave activity. Moreover, the occurrence of slow waves is often paralleled with a closer relation between functional and structural connectivity. We find that hyperpolarization can also generate unconscious-like synchronized Up and Down states, which may be a mechanism underlying the action of anesthetics. We conclude that the TVB-AdEx model reveals large-scale properties identified experimentally in sleep and anesthesia.

BOLD-based fMRI is the most widely used method for studying brain function. The BOLD signal while valuable, is beset with unique vulnerabilities. The most notable of these is the modest signal to noise ratio, and the relatively low temporal and spatial resolution. However, the high dimensional complexity of the BOLD signal also presents unique opportunities for functional discovery. Topological Data Analyses (TDA), a branch of mathematics optimized to search for specific classes of structure within high dimensional data may provide particularly valuable applications. In this investigation, we acquired fMRI data in the anterior cingulate cortex (ACC) using a basic motor control paradigm. Then, for each participant and each of three task conditions, fMRI signals in the ACC were summarized using two methods: a) TDA based methods of persistent homology and persistence landscapes and b) non-TDA based methods using a standard vectorization scheme. Finally, using machine learning (with support vector classifiers), classification accuracy of TDA and non-TDA vectorized data was tested across participants. In each participant, TDA-based classification out-performed the non-TDA based counterpart, suggesting that our TDA analytic pipeline better characterized task- and condition-induced structure in fMRI data in the ACC. Our results emphasize the value of TDA in characterizing task- and condition-induced structure in regional fMRI signals. In addition to providing our analytical tools for other users to emulate, we also discuss the unique role that TDA-based methods can play in the study of individual differences in the structure of functional brain signals in the healthy and the clinical brain.