Neuroinformatics

Three-dimensional (3D) whole-organ imaging and analysis at cellular resolution (termed 3D histology) provide profound insights into the organization and interactions of cells throughout organs. However, the quantitative analysis of these massive datasets remains a significant bottleneck due to the lack of integrated, user-friendly tools. Here, we present 3DBrainOne, an end-to-end ImageJ plugin that streamlines the essential 3D histological analysis of the mouse brain--from raw image preprocessing to region-wise quantification--within a single platform. 3DBrainOne features a robust whole-brain cell-counting module that uses a Difference-of-Gaussians (DoG) blob detection algorithm followed by a ResNet18-based deep learning classifier, enabling high-fidelity automatic whole-brain cell counting with a graphical user interface (GUI) for visual inspection and manual curation of analysis results. 3DBrainOne also supports multi-channel colocalization analysis. Furthermore, this platform includes modules for atlas alignment and brain-region-wise volumetric quantification, enabling brain region-resolved cell counting and structural analyses. As an ImageJ plugin, 3DBrainOne is compatible with a range of operating systems and hardware. In summary, 3DBrainOne is an integrated, versatile, and easy-to-use platform that will facilitate 3D histological analyses in experimental neuroscience.

We present a free and open-source, semi-automated, topologically robust pipeline for fitting cable models to 3D surface mesh morphology data of neuronal membranes, particularly suited to structures with complex shapes and topological holes. The motivation for this work is the discovery of morphologically complex neural spines on the auditory space-specific neurons of the barn owl (Tyto alba, Tyto furcata), dubbed "toric spines", notable for their high curvature, branching density, and holes/loops. Multicompartmental simulation software requires morphology to be represented as cable models (e.g., SWC format), yet existing software tools for fitting cable models to complex 3D surface meshes have not produced satisfactory results for toric spines, and loops are generally unsupported. We present the Mesh and Skeleton Cable Fitting (MASCAF) pipeline and software, which fits a cable model (e.g., SWC format) to a surface mesh using mean-curvature flow skeletonization. In this paper, we demonstrate how MASCAF is applied to fit cable models, how loops can be reconstructed in simulations with the Arbor and NEURON simulation software, and how the results can be validated using geometry and simulator-based methods. While non-tree morphologies such as toric spines are neuroanatomically special, our software pipeline provides a cable-model fitting approach for surface mesh data that is topologically robust, deterministic, open-source, and applicable to general morphologies, thereby closing a crucial gap between neuronal imaging and high-resolution simulation.

One of the prominent methods in neuroimaging data processing is SSM-PCA, which is based on principal component analysis and allows for the identification of diagnostically significant patterns in the form of statistical maps. We developed software, PIE Toolbox, employs SSM-PCA and classification based on the obtained diagnostic patterns revealed from functional and structural tomographic brain imaging. The program supports the entire analysis pipeline including preprocessing of brain images, diagnostic patterns extraction, building classification models, and prediction based on them. The resulting diagnostic patterns are weighted principal components obtained through SSM-PCA, or their linear combinations. PIE Toolbox allows selection of relevant structural and functional brain patterns, computation of their expression values in regions of interest, classification using support vector machines, and evaluation of model performance via cross-validation. This approach enables the use of patterns as features of intergroup differences for individual diagnosis. The software has been validated on both simulated and ADNI datasets.

Human brain tissue preserved in biorepositories is foundational for the structural, cellular, and biomolecular research necessary for a mechanistic understanding of neurological diseases. Realizing the research potential of these valuable resources requires well-characterized research-relevant tissue that can be efficiently identified by investigators and incorporated into the conceptual and computational frameworks of interdisciplinary research. Several large-scale efforts to improve research reliability and reproducibility have sought to characterize and annotate the processes by which these samples are collected, yet limited progress has been made on standardizing spatial information for these samples. Biorepositories systematically collect brain tissue according to a brain sampling protocol (BSP) that differs between institutions, yet explicit spatial information regarding the samples may not be documented in standard operating procedures (SOPs). The amount of anatomical location details available to investigators are inconsistent across biorepositories and typically lack sufficient anatomical precision to ensure correspondence with samples from other biorepositories or research relevant brain regions specified by neuroimaging, functional, or disease-susceptibility criteria. Here, we introduce a pipeline for developing a Spatial Atlas for Mapping Protocol Locations of Ex vivo Samples (SAMPLES), which uses a neuroimaging framework to create a 3D representation of a BSP through a metrically precise digital instantiation of the procedures for brain extraction, segmentation, slicing, and sampling on a modern digital brain template. SAMPLES incorporates modern neuroinformatics conventions to create explicit 3D labels of BSP-defined samples that can be interactively visualized with freely available neuroimaging software. We illustrate the pipeline by developing an atlas for the protocol from the University of Washington BioRepository and Integrated Neuropathology laboratory (UW BRaIN SAMPLES). By providing an explicit, computable reference, SAMPLES atlases can support the efficient identification, referencing, and utilization of postmortem samples for interdisciplinary research. These capabilities enable biorepository workflows, data harmonization across biorepositories, and integration with antemortem neuroimaging.

Despite growing understanding of the benefits of having Findable, Accessible, Interoperable, and Reusable (FAIR) data, many datasets still cannot be shared. Federated analysis methods can enable multisite studies that do not require the sharing of participant-level information. However, there are many practical hurdles that prevent the large-scale adoption of federated methods. We discuss challenges related to cross-site data preparation for federated learning, present solutions offered by recent neuroinformatics projects, and showcase an example of tool integration applied to neurodegenerative disease data.

Currently, there is no consensus about how investigators should format their NIBS data for sharing. This presents a barrier to the advancement of big data analyses because it requires time-consuming operations to generate consistent formats across different shared datasets. Recently, we launched Big non-invasive brain stimulation data (Big NIBS data), an open-access platform and repository for NIBS data (https://www.bignibsdata.com/), providing a structured mechanism for researchers to share NIBS data. However, the reusability and interoperability of data uploaded to Big NIBS data is restricted by the absence of a common data structure. The current paper addresses this problem by creating the NIBS data analysis structure (NIBS-DAS), a template pipeline for the layout, management, and analysis of collated NIBS outcome data. While its primary purpose is to provide a template layout for uploading collated data to the Big NIBS data repository, NIBS-DAS also offers guidelines for the management and analysis of collated NIBS data, thereby forming a data analysis pipeline that can be freely used by the NIBS field in general. We anticipate that NIBS-DAS will serve to facilitate data sharing on the Big NIBS data platform and promote greater standardisation of data management and analytical practices in the NIBS field.

BackgroundThe relevance of spontaneous activity has been unlocked thanks to recent large scale recordings that revealed, via Shared Variance Component Analysis (SVCA), the high-dimensional nature of the ongoing activity. A fundamental problem is how the dimension modifies when more neurons are included in the analysis. Contradictory results have been reported on this subject based on SVCA and Principal Component Analysis (PCA). New MethodWe investigate pro et contra of SVCA and PCA for the identification of reliable responses encoding underlying state variables. We focus on common features of the spectra of the reliable variances (RVs) and on their dimensionality. The analysis is demonstrated on previously published Ca2+ data from the visual and the dorsal cortex in head fixed mice during spontaneous behavior. ResultsRVs grow proportionally to the number N of neurons and show a power-law decay k- with the k-th SVC dimension over a range bounded by a maximal dimension kc, initially diverging as N 1/ and then saturating at sufficiently large N. The reliable dimensionality, estimated with different methodologies, also shows a clear saturation to an asymptotic value for large N. Furthermore, its value decreases when becomes larger, as demonstrated by employing experimental data as well as theoretical predictions. ConclusionWe have shown that SVCA is an extremely effective tool to extract reliable features from the neural signals, and that the exponent represents a biomarker able to reveal the level of correlation of the neurons as well as the dimensionality of the reliable space. HighlightsO_LIAdvantages and drawbacks of Shared Variance Component Analysis to extract reliable signals from neural data C_LIO_LIComparison of different methods to estimate reliable neural dimensionality associated to spontaneous activity C_LIO_LIAnalytical expressions of embedding dimensionality for power-law decaying reliable variances C_LIO_LIBounded growth of the dimensionality with the number of neurons C_LI

Electroencephalography (EEG) interpretation in clinical practice relies on the analysis of energy distribution across standard frequency bands. The Weed Plot framework encodes band-wise spectral energy, computed using Fourier-based methods, into a symbolic representation that preserves the interpretability of traditional EEG analysis. In this study, we propose a wavelet-based extension of this framework, where the energy of predefined clinical EEG bands is estimated using the Discrete Wavelet Transform instead of Power Spectral Density. Unlike Fourier-based approaches, wavelets provide a time-frequency representation that captures transient and non-stationary dynamics while remaining consistent with clinically defined bands. From these estimates, symbolic patterns are constructed based on the relative ordering of frequency bands within short temporal windows. Their empirical distribution is used to extract entropy-based features for epilepsy detection using multiple machine learning classifiers. From an Artificial Intelligence perspective, the main contribution is a structured symbolic encoding that enhances feature discriminability. From an engineering perspective, the contribution lies in an automated framework for EEG-based epilepsy detection. Experimental results show that wavelet-based representations improve classification performance compared to raw entropy and Fourier-based features. This improvement arises from the interaction between time-frequency localization and symbolic encoding, producing more discriminative feature distributions. These findings support wavelet-based symbolic representations as a robust and interpretable framework for EEG analysis, bridging clinical interpretation and data-driven methods.

Sex-related differences in the aging of the human brain were studied using large array of experimental data. The open archive CamCan was used as a source of data: the magnetic encephalograms, co-registered with magnetic resonance images of the head, were obtained for each of 434 subjects (ages 18-87 years, mean age 54.7 {+/-}18.4): 217 females (ages 18-87 years, mean age 54.5 {+/-}18.4) and 217 males (ages 18-84 years, mean age 54.8 {+/-}18.3). Recordings were split in 10-year age cohorts, each cohort consisted of equal number of men and women to calculate average intersex characteristics correctly. By massively solving the inverse problem, functional tomograms were calculated - the spatial distribution of elementary spectral components. Physiological noise was eliminated by joint analysis of MEG-based functional tomogram and magnetic resonance image for each subject. Then multichannel spectra were transformed into time series of the power of elementary current dipoles. Summary electric powers were calculated in six conventional frequency bands (1-4 Hz - delta; 4-8 Hz - theta; 8-13 Hz - alpha; 13-21 Hz - beta1; 21-30 Hz - beta2; 30-48 Hz - gamma), and sex differences in age-related changes were examined. It was found that in the youngest age cohort (18-29 years) the summary electrical power of the brain for males is 1.5 times greater than such power for females. For adults (30-69 years), male and female powers are approximately equal, while in older cohorts (70-87 years), male total brain power is greater. Age dependencies in various frequency bands are generally different for men and women, excluding higher frequencies 21-48 Hz. Basic conclusion can be made that after intersex averaging total electric power of the human brain is invariant through the lifespan from 18 to 87 years. The proposed method of joint MEG and MRI analysis can be used for further study of the sex-related details of brain sources in their connection with age changes.

Neuroimaging data management and processing are tedious and error-prone, prompting reproducibility concerns. Globally, studies with heterogeneous infrastructure and governance policies lead to eclectic data processing and sharing, necessitating standardization of data workflows to ensure reusability and comparability of multi-centric datasets. The Nipoppy neuroinformatics framework facilitates such standardization by combining specification, protocol, and software to manage study-level data workflows. With its adoption, researchers can share standardized, derived datasets enabling efficient, reproducible, and inclusive research.

BackgroundIntegrating multimodal data into medical artificial intelligence (AI) tools and evaluating whether they outperform human experts remains a critical challenge. Epilepsy surgery offers a unique paradigm for this evaluation, as it provides an expert-independent measure (Engel score) of post-surgical outcome. Currently, evaluation for epilepsy surgery relies on the visual interpretation and human synthesis of multimodal data. While clinical evaluations are individualized and account for complex anatomical variability, integrating these diverse, high-dimensional modalities to generate a probability of surgical success remains challenging. Here, we leverage this objective outcome score to investigate the feasibility of a data-driven, phenotype-based model against the current clinical gold standard. MethodsThe evaluation was performed on an epilepsy-type controlled cohort of 57 patients from six tertiary epilepsy surgery centers who underwent resective/ablative surgery in the mesiotemporal lobe. Multimodal data, namely, patient demographics, semiology, invasive electrophysiology monitoring, and neuroimaging, were utilized. We first estimated how human experts perceive surgery success. Subsequently, we developed a data-driven model integrating these modalities to predict surgery outcomes. The model performance was compared to the current clinical gold standard (three independent human experts) and published outcome calculators. Finally, modality-level phenotypes were derived based on the models predictions. ResultsPredictions by human experts correlated poorly with post-surgical outcomes, and published outcome calculators did not perform better than the experts (DeLongs p = 0.367). Our model incorporating multimodal data achieved an area under the receiver operating characteristic curve (AUROC) of 0.801. It performed statistically better than the best human expert (DeLongs p = 0.043) and achieved a higher AUROC than the best published surgical outcome calculator (0.801 vs. 0.694). ConclusionsWe demonstrated the proof-of-concept that data-driven multimodal phenotypes can inform personalized surgery planning in epilepsy. Furthermore, we provide a framework for integrating multimodal data and benchmarking medical AI performance against human experts.

The neuroscience literature contains thousands of studies localizing cognitive, sensory, and motor functions to specific brain regions, yet this knowledge remains fragmented across experimental modalities, naming conventions, and spatial reference systems. Consequently, relating reported activations, lesions, or stimulation sites to the broader functional literature often requires substantial manual synthesis. The Brain Encyclopedia Atlas Project (BEAP) was developed to address this challenge by providing a spatially grounded framework for organizing literature-defined brain regions. BEAP is an expert-curated neuroinformatics resource that aggregates and spatially indexes literature-defined cortical and subcortical functional regions within a common anatomical reference framework. The project identifies 108 neocortical fields and 18 cerebellar fields defined through an analysis of published figures from 1,453 human studies using functional neuroimaging, intracranial electrophysiology, and cortical stimulation. These regions were manually aligned to standard anatomical templates and associated with parcels of the Human Connectome Project multimodal parcellation (MMP1). Inclusion criteria required convergent functional evidence, lesion support, and boundary-related contrasts. Additionally, 340 allocortical, diencephalic, cerebellar, and brain stem nuclei were delineated through comparison with histological atlases and research articles. The resource is publicly accessible at https://brainatlas.online/3d-brain/, featuring an interactive three-dimensional brain model that interfaces directly with a curated encyclopedia. This platform provides structured entries synthesizing regional functional descriptions, boundary-defining evidence, internal organization, and connectivity annotations. Furthermore, each entry is designed to evolve through community feedback via a dedicated comment section. By providing a unified spatial context at the whole-cortex scale, BEAP enables systematic comparison across studies and facilitates the identification of recurring patterns in cortical organization. It serves as an integrative resource for research and education, supporting the contextualization of neuroimaging findings and the generation of hypotheses regarding large-scale brain organization.

BackgroundThis study examines a competition-based model (C-model) designed to capture the temporal dynamics of successive brain microstates derived from electroencephalography (EEG) recordings during eyes-open conditions. The analyzed data were obtained from a public repository comprising microstate sequences from 60 sessions of a single subject [1]. When applied to microstate dynamics, the C-model posits a stochastic competition among neural circuits underlying the expression of individual microstates. MethodsThe model is formulated at a conceptual level (computational level in Marrs framework) and employs a geometric distribution to account for the long right tail of microstate duration distributions, interpreted as the probability of "failure" of the currently active microstate to persist. To account for the short-lived left tail, the model incorporates a transient increase in the stability of the currently active network, or equivalently, a temporary decrease in the activation probability of competing microstates (refractory period). ResultsThe model provides a good fit to the microstate duration distributions across all 60 sessions. One third of sessions showed microstate identity sequential dependency with respect to the previous microstates. DiscussionThese results suggest that the C-model captures key aspects of microstate temporal structure. Moreover, because microstate probabilities can be modulated by psychophysiological conditions--including the influence of previously active networks--the model may serve as a building block for more comprehensive neurobiological frameworks of neural and behavioral dynamics. In such frameworks, microstate sequences could emerge from structured competition and flow among neural networks supporting microstate expression.

The human cerebellum plays a central role in motor, emotional, and cognitive functions, and is implicated in many brain disorders. To improve the analysis of functional and anatomical imaging from the cerebellum, we introduce SUITPy, an improved and fully revised Python implementation of the widely used SUIT toolbox. For this new version, we developed a U-Net based model to automatically isolate the cerebellum from adjacent cortical tissue, which achieves higher fidelity than existing algorithms. The isolation works robustly without manual corrections for imaging data across the lifespan. We show that isolation and subsequent normalization to a cerebellum-only template lead to a more precise alignment of cerebellar structures across participants compared to normalization using a whole-brain template. We also show the utility of the cerebellar mask to prevent contamination of cerebellar functional data from surrounding cortical structures. The toolbox also provides functionality for visualizing cerebellar data on a flatmap, along with a range of anatomical and functional cerebellar atlases, thereby offering an essential tool that enables accurate cerebellar analysis across the lifespan.

Kilosort 2.0 is a widely adopted spike sorting algorithm recognized for its efficiency and accuracy on planar electrode arrays, such as Neuropixels. To adapt its robust architecture to emerging three-dimensional (3D) neural probes, we present Kilosort 2.0-3D, a modified version that leverages 3D spatial information. Our modification specifically redefines the spatial processing components of Kilosort 2.0 to operate in 3D space while leaving the core template-matching process unchanged. By using synthetic extracellular recording data with ground-truth neuron positions and firing times, we demonstrate that Kilosort 2.0-3D effectively resolves spatial ambiguities and unit misclassifications inherent in 2D spatial assumptions. Our results show that Kilosort 2.0-3D achieves rotational invariance and maintains full backward compatibility with planar arrays. This work establishes a validated, scalable tool for spike sorting of high-density 3D neural electrophysiology data.

Ovarian cancer is one of the gynecological cancer types, which, if metastasized and not detected early, can cause deaths among women. Therefore, there is a need to accurately predict drug responses to ovarian cancer. A gynecological pathologist inspects abnormality in tissues, followed by providing a report about patients; however, such a diagnostic process is (1) hard; (2) requires experience; and (3) time consuming. Moreover, existing tools are far from perfect. Hence, we present a computational pipeline to improve predicting drug response pertaining to ovarian cancer, derived as follows. First, we download digital pathology images pertaining to ovarian bevacizumab response from the cancer imaging archive repository. We employed histogram of oriented gradients to images, constructing feature vectors, provided to Fisher linear discriminant analysis to change the representation through dimensionality reduction. Then, we provide reduced-dimensionality data for regression analysis through support vector regression coupled with various kernels and calculating the area under the ROC curve (AUC). Experimental results against transformer-based models (ViT and Swin) and other deep learning (DL) models (VGG16, ResNet50, InceptionV3, MobileNetV2, and EfficientNetB6) demonstrate that our approach with radial kernel (named SVRD+R) yielded an AUC performance improvements of 17% against the best-performing transformer-based model (ViT) while obtaining an AUC performance improvements of 14.9% when compared against the best DL-based model (MobileNetV2). These results demonstrate the superiority and feasibility of our AI-based pipeline when tackling prediction problems pertaining to gynecologic cancer studies. MSC92B05; 68T09

ObjectivesStage-sensitive markers can aid in early diagnosis of Alzheimers disease (AD) and can improve sensitivity, performance and interpretability. In this study, causal markers from longitudinal imaging data were extracted and integrated with risk factors to improve diagnostic models. Data DescriptionOASIS-3, a longitudinal dataset consisting of 613 controls and 214 cases with very mild to moderate Alzheimers disease is used for this study. A meta model was built using a predisposition model built from risk factors, a stage-sensitization model built from MRI markers at various stages of atrophy and a confirmatory model built using PET markers. The meta model achieved good diagnostic performance (accuracy = 93%, sensitivity = 80%, specificity = 95%). Exclusion of PET data achieved comparable performance (accuracy = 91%, sensitivity = 85%, specificity = 92%). The results demonstrate that integrating causal pathological markers with risk factors improves diagnosis and aids in elucidating stage-specific patterns of AD.

Accurate brain and cerebrospinal fluid (CSF) volume assessment is essential for pediatric hydrocephalus management. Current clinical practice relies on linear measurements that fail to capture complex three-dimensional ventricular morphology, while quantitative volumetric assessment remains limited by laborious processing and lack of clinically optimized automated tools. This study developed a rapid, automated AI-based intracranial segmentation model suitable for clinical workflows. We retrospectively analyzed 167 T2-weighted MRI scans from infants with hydrocephalus, randomly split into training (60%), validation (20%), and hold-out test (20%) sets. All scans were manually segmented into CSF, brain parenchyma, and background. Our model integrates DenseNet and U-Net architectures with feature smoothness regularization to enhance generalizability. Performance was evaluated using Dice scores and absolute relative volume error (ARVE) compared with state-of-the-art methods. The AI model achieved Dice scores of 95.7% for CSF and 96.4% for brain parenchyma on the hold-out test set, significantly outperforming FSL FAST (85.0% and 77.9%) and contemporary deep learning approaches (90.4% and 89.7%). Processing time was 0.8 seconds per scan using GPU acceleration. The model demonstrated consistent performance across different hydrocephalus etiologies and effectively handled challenging scenarios including noise, artifacts, and variable resolution. This study successfully developed a robust MRI segmentation model demonstrating superior accuracy and efficiency compared to existing methods. By incorporating domain-specific enhancements, the model enables rapid, clinically viable brain and CSF volume estimation for pediatric hydrocephalus care.

Brain stimulation is increasingly recognized as an effective and important therapeutic intervention for many brain diseases. Distance between the scalp and other brain regions is a pivotal variable for neurostimulation planning and the development of new techniques, but alterations in the distance between the scalp and other regions in brain diseases are largely unknown. In this study, we developed an automatic pipeline to calculate scalp-to-region distance (SRD) values from T1 MR images and applied it to a total of 1382 participants, including patients with autism spectrum disorder (ASD), Parkinsons disease (PD), Alzheimers disease (AD), and cognitively normal controls (CNs). Cloud points were uniformly sampled on the automatically extracted scalp surface and cortex surface, on which the point-wise distance maps were generated. The brain was then coregistered with the BCI-DNI atlas, and SRD value for each brain region was extracted. Analysis of covariance (ANCOVA) was performed for SRD in each brain region, with age and sex as covariates. Compared with CNs, ASD patients showed widespread SRD decreases across the brain with prominent involvement of the frontal lobe, especially the orbitofrontal cortex and adjacent regions. In contrast, in AD patients, significantly increased SRD values were observed in various regions of the frontal gyrus. No significant SRD alteration was found in PD patients after correction. The automatic SRD calculation pipeline and the different patterns of SRD alterations in these diseases might be helpful for future neurostimulation planning in clinical practice. HighlightsO_LIAutomatic pipeline enables scalp-to-region distance (SRD) measurement, facilitates brain stimulation planning. C_LIO_LIASD patients show widespread SRD decreases, especially in the orbitofrontal cortex and adjacent regions. C_LIO_LIAD patients present increased SRD in the frontal gyrus and decreased SRD in the parahippocampal gyrus. C_LI

IntroductionTo leverage high-frequency oscillations (HFOs) as a biomarker with significant potential, this study compared a large set of detectors on a unified dataset, aiming to evaluate their clinical applicability under realistic conditions. MethodsEleven automatic detectors were applied to a retrospective dataset of intracranial and scalp EEGs from 27 consecutive pediatric patients. Inter-detector agreement was assessed using Spearmans Rho, and the area under the curve (AUC) for seizure onset zone (SOZ) prediction served as a consistent reference standard to enable reliable comparisons across recording modalities. Analyses were conducted separately for HFO and Spike-HFO detections. ResultsThe average age of our cohort was 12.4 years (SD 4.0; range 5-18). AUC values in scalp EEG ranged from 0.61 to 0.67 for HFOs and from 0.53 to 0.63 for Spike-HFO. AUC values in intracranial EEG ranged from 0.48 to 0.66 for HFOs and 0.54 to 0.69 in Spike-HFO. Although only three of the 11 detectors were specifically developed or adapted for scalp EEG, the detectors generally achieved higher AUC values and stronger agreement in scalp EEG ConclusionsWe present the first study comparing intracranial and scalp detectors by testing them beyond the modalities for which they were originally designed. Although the clinical utility of detections was comparable across EEG modalities, it remained lower than reported in original studies assessing the diagnostic value of HFOs. Caution is warranted when applying a publicly available detector to a new dataset, and detector robustness remains a critical issue. Key points- A comprehensive head-to-head comparison of 11 detectors demonstrated significant variability in detector agreement and clinical utility - Clinical utility was not necessarily linked to the EEG recording type the detector was originally designed for - Despite widely accepted use of automatic detections, detector robustness remains a critical issue