Neuroinformatics

Mean-field models are extensively used in large-scale brain simulations because they provide a wieldy description of population dynamics while preserving key features of neural activity. Despite their widespread adoption, no common and reproducible methodology currently exists to systematically derive and validate mean-field models starting from biologically grounded single neuron dynamics. As a result, implementations are often ad hoc, difficult to reproduce and rarely reusable. Here we introduce BRIDGE, a modular, open-source Python pipeline that enables the bottom-up reconstruction, analysis, validation, and simulation of mean-field models from single neurons. The framework integrates single neurons modelling, network simulations, extraction of population statistics, parameters analysis, quantitative comparisons between spiking neural networks and corresponding mean-field representations, and simulation of network of mean-fields. Its flexible architecture allows users to incorporate different neuron models and to generate region-specific or state-dependent mean-field formulations. BRIDGE provides a reproducible foundation for developing biologically informed mean-field models suitable for large-scale and whole-brain simulations, supporting the transition from generic homogeneous population models toward region-specific ones. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=91 SRC="FIGDIR/small/742067v1_ufig1.gif" ALT="Figure 1"> View larger version (28K): org.highwire.dtl.DTLVardef@1124404org.highwire.dtl.DTLVardef@2f8b2aorg.highwire.dtl.DTLVardef@1598f37org.highwire.dtl.DTLVardef@c9814b_HPS_FORMAT_FIGEXP M_FIG C_FIG

A major goal of the BRAIN Initiative Cell Atlas Network (BICAN) is to create a suite of foundational reference cell atlases and associated standards for human and non-human primate brains. Central to this goal is the creation of cross-species harmonized cellular taxonomies and structural parcellations with formal ontologies that can be mapped into 3D reference frameworks bridging neuroimaging and cellular and histological resolutions. We describe here an iterative approach, focused initially on the basal ganglia, to co-create structural and cellular ontologies in human, macaque and marmoset brains, including a Harmonized Ontology of Mammalian Brain Anatomy (HOMBA), and to map and refine structural parcellations into neuroimaging-based common coordinate frameworks. These references provide the framework for documenting and mapping all experimental sampling in BICAN, allowing analyses of cellular and molecular variation as a function of topographic position, and enabling comparisons of cellular, molecular and neuroimaging-based functional variation within and between primate species. HighlightsO_LIA hierarchical Harmonized Ontology of Mammalian Brain Anatomy (HOMBA) covering 2348 structures C_LIO_LIHOMBA-annotated 3D common coordinate frameworks (CCFs) of the basal ganglia across species C_LIO_LIHistologically informed 3D parcellation/atlas of 280 human subcortical structures indexed by HOMBA C_LIO_LIMapping and integration of structural, cellular and functional data with HOMBA and CCFs C_LI

Converting medical images into anatomically detailed, subject-specific finite element (FE) models is a long-standing bottleneck in brain computational modelling. These models are used to predict brain tissue deformation, e.g. in traumatic brain injury, particle diffusion in brain drug delivery, and other biophysical phenomena across neurological disorders. However, existing model creation workflows depend on manual image segmentation, proprietary meshing software, and labour-intensive repair of meningeal and interface structures, limiting reproducibility and cohort analysis. Here we present PARS, a fully automated, open-source pipeline that converts a T1-weighted MRI scan into a simulation-ready FE head model. PARS combines anatomical parcellation with tissue maps and uses iterative neighbourhood-based reclassification, yielding a gap-free whole-head label volume. The volume is directly converted into a hexahedral mesh, augmented with algorithmically reconstructed falx, tentorium, pia and dura mater, and refined by Laplacian smoothing under a node-locking scheme that controls element quality and the explicit-solver stable timestep. We evaluated PARS on 23 subjects spanning cranial volumes of 832 to 1,329 cubic centimeter, at 1.0, 1.5 and 2.0 mm MRI resolutions. At 1 mm, meshes achieved a median Scaled Jacobian of 0.976, and total intracranial volume error of ~0.54; quality remained high at 1.5 mm (SJ of 0.933) and 2 mm (SJ of 0.921). Model creation runtime ranged from 9 to 38 minutes per subject. Models generated by PARS have been validated against cadaveric brain displacement data and demonstrated utility across traumatic brain injury and normal pressure hydrocephalus research. PARS provides an open-access, reproducible resource that substantially lowers the barriers to subject-specific brain modelling.

Glioblastoma (GBM) is the most common and aggressive primary malignant brain tumor in adults with extremely poor prognosis. Complete surgical treatment is practically impossible, as the true extent of GBM infiltration cannot be fully delineated using currently available in vivo neuroimaging methods, leading to frequent recurrences and low overall survival. Consequently, mathematical models are being developed to estimate the GBM expanse beyond the visible tumor mass, providing additional information for treatment planning and patient prognosis. Here, a novel graph-based stochastic mathematical model of GBM invasion using patient-specific structural brain connectome data is proposed. The model is assessed using publicly available UCSF-PDGM dataset to demonstrate GBM invasion dynamics across multiple patients and anatomical locations. Additional scaling using fractional anisotropy (FA) is tested and evaluated. Parameter sensitivity analysis is provided to explore model's behavior under different settings. Ablation testing is performed to suppress model mechanisms utilizing the structural connectome, showing that the tentacle-like extrusions from the tumor core emerge only if the patient-specific connectome is utilized. The model seems to capture GBM micro-infiltration along white matter tracts to a very high degree, making it a potential tool for studying distant recurrences farther from the resection cavity and GBM invasion dynamics in relation to the structural connectome. Full source code is publicly available, ensuring complete transparency of the study.

Advances in high-throughput mesoscale microscopy and machine learning-based image analysis pipelines have made unbiased whole-brain imaging widely accessible. However, translating the resulting atlas-mapped datasets into biologically meaningful results remains a substantial barrier owing to their sheer magnitude and complex hierarchical organization. Consequently, reporting structure and analysis methods vary widely across studies, under-mining rigor and reproducibility. To address this, we developed a user-friendly data reduction workflow, HERO (Hierarchy-aware Expression Region Organization), designed to perform hierarchy-aware selection, refinement, ranking, and visualization of whole-brain cell detec-tion datasets. The workflow is customizable to specific needs, requires minimal coding expe-rience, and outputs transparent, curated results. HERO is designed to function as a seamless plug-in within larger-scale whole-brain cell-detection analysis pipelines, providing efficient, unbiased region selection to streamline subsequent statistical analyses and comparative evaluations. Although HERO is developed with mouse cell-detection datasets, it can, in prin-ciple, be applied to any atlas-mapped dataset that contains hierarchical information. In sum, HERO offers a standardized analysis workflow to reduce whole-brain cell-detection datasets, transforming raw regional cell counts into curated results and advancing the effectiveness, interpretability, and accessibility of whole-brain imaging in neuroscience.

AO_SCPLOWBSTRACTC_SCPLOWThe integration of electrophysiological recordings with multimodal neuroimaging data holds great promise for advancing our understanding of brain function and neurological disorders. To facilitate this endeavour, we present electro-MICA, an open-access Python toolbox designed to project electrophysiological features from scalp and intracranial electroencephalography (EEG) onto cortical and hippocampal surfaces generated by validated multimodal imaging ecosystems. The toolbox comprises two pipelines: one for intracranial EEG (iEEG) recorded with stereo-EEG depth electrodes, and one for scalp EEG source localization. Both pipelines are grounded in numerical solutions to the electromagnetic equations governing electric activity in the brain, solved using the Boundary Element Method. A key methodological contribution is the use of a current density double layer model for neural generators, which avoids the mathematical singularities introduced by conventional dipole-based models when electrodes are near the cortical surface, a situation that can arise in iEEG. Electrode contacts are additionally modeled with non-zero length, improving physical realism. Scalp EEG source localization is performed using eLORETA on a subject-specific three-layer head model derived from the anatomical input. Validation against empirical gamma-band iEEG data from 32 subjects demonstrates that the distributed generator model outperforms both distance-based and dipole-based alternatives. An illustrative clinical example demonstrates the toolboxs capacity to reveal associations between intracranial spike rates, cortical thickness, and anatomical connectivity in an epilepsy patient. Electro-MICA requires no parameter selection from the user, facilitating straightforward multimodal analyses in both research and clinical settings. The toolbox is available at github.com/MICA-MNI/electromica with extensive online documentation at electromica.readthedocs.io.

Early detection of schizophrenia (SZ) remains challenging due to the subtlety of early-stage brain alterations and reliance on subjective clinical assessment. We propose a frequency-aware 3D convolutional neural network (CNN) pipeline that integrates NeuroMark-HiFi high-pass spatial filtering with a modified VGGNet3D architecture featuring 3D Laplacian kernel initialization and dilated convolutions. Using the FBIRN dataset (N=311; 150 healthy controls, 161 SZ) with all 53 intrinsic connectivity networks (ICNs) per subject, we evaluate four experimental conditions across two hyperparameter configurations to isolate the contributions of enhanced input representations and frequency-aware model design. Under the optimized configuration, Condition 3 (HiFi + Laplacian initialization) achieved the best mean test accuracy of 75.54% with a peak single-fold accuracy of 87.10%, representing a 5.44% absolute gain over the optimized baseline. These results demonstrate that high-frequency spatial features are more discriminative for SZ classification than raw intensities, and that aligning Laplacian-initialized kernels with HiFi-filtered input creates a beneficial inductive bias--even with a compact model of approximately 1.4M parameters.

EegFun.jl is an open-source package for electroencephalography (EEG) analysis implemented in the Julia programming language. EegFun.jl provides a flexible framework for EEG research, covering data import from standard file formats, filtering and re-referencing, Independent Component Analysis (ICA) for artifact detection/correction, epoch extraction, and ERP averaging and visualisation. The Julia language provides the readability of a high-level scripting environment together with execution speeds comparable to compiled code. EegFun.jl combines interactive data visualization with high-performance execution, making large-scale analyses both efficient and easy. Here, we provide a brief overview and introductory tutorial of the core stages of the EEG analysis workflow to illustrate the packages capabilities. The package is freely available under the MIT license.

Understanding vascular contributions to disease is a major unmet need. White matter lesions (WML) are an accepted imaging marker of cerebral small vessel disease, giving insights into its related pathologies. A unified approach for WML analyses in large multi-site data is lacking despite the need for pooling of data to overcome the limitations of often small heterogenous MRI studies which make subtyping and identifying patterns within disease groups difficult. Our ENIGMA-PD-WML pipeline is an open-source containerized pipeline containing all the code and packages required for pre-processing, processing and post-processing of T1-weighted and FLAIR data, outputting accurate and reproducible binary WML maps using a UNet approach. The pipeline provides a standardized image analysis approach for WML and outputs data in both native and MNI space to allow for sharing and pooling of data from multiple sites for large-data analysis. In addition to a reliable standardized approach for WML segmentation, key priorities when developing the pipeline included: usability, i.e., requiring minimal manual input and technical expertise to use, and suitability to run on various MRI scanners and acquisition parameters as is common in multi-site data. This paper describes the pipeline in detail, with rationale for each step, providing transparency and facilitating its usage to overcome reproducibility issues in large-scale WML analyses.

The corpus callosum is a major white matter bundle responsible for connecting both hemispheres. In mammals, due to a variety of causes, the development of the corpus callosum can be impaired - this brain malformation is known as corpus callosum dysgenesis (CCD). The clinical presentation of CCD varies, with patients exhibiting three morphological phenotypes: agenesis, partial dysgenesis, and hypoplasia. Although the first two presentations are easily detectable on MRI scans, the latter is more challenging, as the structure is fully formed but has a reduced area. In this study, we develop (1) a pipeline to generate synthetic MRI scans with apparent anatomical variation and (2) train a U-Net-based tool to automatically segment the corpus callosum of marmosets in both healthy and disease contexts. Methodologically, a custom script was devised to apply rotation and translation to T1-weighted MRI scans at the volume level. Because the slicing grid remains unchanged, these rigid transformations translate into apparent anatomical variations at the slice level. We compared corpus callosum measurements obtained from automatically segmented masks with those from manually delineated masks. The average Dice score was above 0.90, and the Hausdorff distance was below 0.4 mm. We also stratified our cohort according to phenotype (healthy controls and hypoplastic animals). The magnitude of the effect and the significance level observed between the voxel counts of healthy and hypoplastic animals using manually delineated masks were comparable to those obtained via automatic segmentations. These results show that our pipeline can generate a sufficiently varied training pool to build an accurate U-Net segmentation model with high diagnostic capability.

In patients with epilepsy, seizures are associated with pathological neural synchronization. However, the preictal period preceding a seizure often exhibits reduced spatial synchronization compared to normal cognition. This observation aligns with the concept of the brain as a complex dynamical system, where a reduction in dimensionality and resilience can precede a phase transition. The Critical Brain Hypothesis suggests a connection between the loss of healthy scale-free behavior and various disorders, including epilepsy. Our study investigates preictal changes by utilizing network features, such as mean node degree and mean clustering coefficient, derived from thresholded correlation matrices of patient intracranial electrocorticographic electrode data. We observed a suppression of intermittent high-synchronization periods within the feature space during the minutes leading up to seizure onset. This constriction of the explored hypervolume in the preictal state indicates a breakdown in the brains ability to maintain normal coherence. We use these preictal changes to predict the probability of seizure onset using a Support Vector Machine algorithm. These discrete predictions can then be combined into real-time continuous seizure risk forecasts via Bayesian updating. This innovative and computationally lightweight approach has the potential to significantly improve upon static predictions, providing opportunities for more adaptable, quantitative, and interpretable tools for managing seizures.

Mapping hierarchical brain networks within traditional Euclidean space causes significant structural distortion, undermining neuroimaging diagnostic frameworks. While hyperbolic models like the Poincare ball preserve these nested topologies, they demand heavy computational overhead due to intricate Mobius operations and curved geodesics. This paper introduces a highly efficient non-Euclidean framework for analyzing neurocognitive decline utilizing the Beltrami-Klein ball model. By projecting hyperbolic geodesics as Euclidean straight lines, this approach converts complex distance calculations into simple dot products, radically reducing processing demands. We validated our methodology against state-of-the-art Poincare and Lorentz baselines using datasets for Schizophrenia, Parkinsons Disease, and Alzheimers Disease. The Klein-based framework demonstrates superior performance, delivering both higher diagnostic precision and accelerated processing velocities across all three neurocognitive disorders.

Traumatic brain injury (TBI) triggers widespread biomarker activation, including astrocytic markers such as glial fibrillary acidic protein (GFAP) and microglia markers such as ionized calcium-binding adapter molecule 1 (IBA1). Quantifying and analyzing these biomarkers are critical for understanding injury impact; however, current methods are labor-intensive and time-consuming. In this study, we propose an automated deep learning framework for dual-biomarker segmentation and TBI classification using GFAP and IBA1 immunofluorescent images. Four U-Net variants: Baseline U-Net, U-Net++, MANet, and LinkNet were trained for segmentation. Three classification models, ResNet50, Swin_T, and MaxViT, were trained to distinguish TBI from control images under single- and dual-biomarker conditions. The baseline U-Net achieved the highest segmentation Dice score for GFAP (0.9259), while the U-Net++ achieved the highest Dice score for IBA1 (0.9676). Trained segmentation models demonstrated significantly better performance compared to QuPath alternatives. While GFAP alone supported high classification accuracy, IBA1 alone was less effective. Multimodal fusion of GFAP and IBA1 significantly improved classification performance across all models, with Swin_T achieving the highest overall accuracy (0.9489), and ResNet50 achieving the highest F1-score (0.9499). These findings demonstrate that integrating complementary biomarkers enhances automated TBI classification, and deep learning offers a robust alternative to manual analysis for immunofluorescent brain injury imaging. This framework is scalable to additional biomarkers and injury models, offering a reproducible approach to accelerate biomarker research.

Accurate alignment of hippocampal anatomy across individuals remains challenging due to complex and highly variable folding patterns that are not well captured by conventional volumetric approaches. HippUnfold introduced a surface-based representation of the hippocampus, but key components--including coordinate estimation and inter-subject correspondence--were defined in the volumetric domain, making them susceptible to topological errors and interpolation artifacts. Here, we introduce a surface-intrinsic formulation of hippocampal unfolding in which geometry, intrinsic coordinates, and correspondence are defined directly on subject-specific surface manifolds. Intrinsic anterior-posterior and proximal-distal coordinates are computed by solving Laplace equations on the surface, and correspondence is established through surface-based resampling in unfolded space, replacing inverse volumetric warping. Relative to the original HippUnfold approach, this formulation improves test-retest consistency, subject identifiability, and mesh quality, while better preserving subject-specific gyral and sulcal morphology. Surface representations show reduced distortion between folded and unfolded spaces and eliminate misplaced or outlier vertices associated with volumetric warping. These improvements translate to enhanced sensitivity in a clinical application, improving lateralization of temporal lobe epilepsy. These results demonstrate that a surface-intrinsic formulation provides a principled and robust foundation for hippocampal unfolding, enabling topology-preserving alignment and more accurate characterization of inter-individual variability in health and disease.

We propose a study protocol for routine clinical electroencephalograms (EEGs) from public hospitals, which represents a vast resource for neuroscience research. These non-invasive measures of brain function, paired with rich clinical annotations from large and diverse patient populations, are critical for developing robust artificial intelligence (AI) models and conducting population-level studies. This protocol presents a scalable methodology for curating and harmonizing extensive clinical EEG datasets, encompassing over 40,000 individual studies, to facilitate research applications. Key steps include: (i) integration of raw EEG recordings with corresponding clinical records, including neurological reports, diagnostic codes, and potentially medication data; and (ii) spatial standardization of EEG signals by mapping them to a common brain space defined by functional and anatomical landmarks. The resulting harmonized datasets enable the development of large-scale EEG foundation models, the discovery of novel EEG waveform representations, and the creation of normative "brain charts" for electrophysiological assessment across the lifespan. By enabling standardised, large-scale analyses of real-world clinical EEG data, this protocol supports data-intensive solutions for EEG applications and addresses the challenge of generalising AI models. Our approach promotes the translation of AI tools from research to diverse patient populations, advancing population neuroscience.

Wave recurrent neural networks (wRNNs) are biologically inspired recurrent architectures that use traveling-wave dynamics to support sequence learning and memory. However, their input-to-hidden pathway remains relatively simple compared with biological neurons, where dendrites perform nonlinear input integration. In this study, we introduce the Dendritic Wave Recurrent Neural Network (DWRNN), which augments the input pathway of the wRNN with nonlinear basal dendritic branches while preserving the original recurrent wave dynamics. We evaluate DW-RNN on a simple copy task, sequential MNIST (sMNIST), permuted sequential MNIST (psMNIST), and noisy sequential CIFAR-10 (nsCIFAR-10). On the copy task, DW-RNN shows learning behavior comparable to the standard wRNN, suggesting that dendritic input integration does not disrupt the recurrent wave-based memory mechanism. On the three sequential image-classification benchmarks, DW-RNN outperforms the standard wRNN, improving accuracy from 97.27 {+/-} 0.15% to 97.82 {+/-} 0.12% on sMNIST, from 96.74 {+/-} 0.17% to 96.92 {+/-} 0.10% on psMNIST, and from 54.30 {+/-} 0.79% to 55.65 {+/-} 0.55% on nsCIFAR-10. In addition to improving mean accuracy, DW-RNN exhibits lower across-seed variability on all three classification benchmarks, suggesting that dendritic input integration may improve the stability of wRNN training. Hidden-activity visualizations further show that DW-RNN preserves the characteristic traveling-wave patterns of the original wRNN. These results suggest that dendritic computation and traveling-wave recurrent dynamics provide complementary mechanisms for biologically inspired sequence learning.

Micro-electrode array (MEA) recordings are widely used to characterize functional connectivity in neural cultures and have gained traction for the analysis of human brain slices. However, the impact of graph construction methodology on the resulting network topology has not been systematically quantified. Here, we benchmark three methods - shared spiking activity, Pearson cross-correlation, and the spike time tiling coefficient (STTC) - across 37 recordings from human cortical slice cultures classified into low, moderate, and high activity groups. We show that method choice alone produces large topological differences (Cohens d = 0.86-1.14 for clustering coefficient, d > 1.0 for node count), while higher-order features such as modularity remain stable. Each method exhibits a distinct sensitivity profile: shared spiking detects activity-dependent changes primarily through network size, correlation uniquely captures clustering differences, and STTC combines strong biological sensitivity with negligible parameter dependence across lag windows (all d < 0.1). Within shared spiking, z-score normalization dominates all other parameter choices (d > 1.0 versus bin size effects of d < 0.23), functioning as an implicit analytical null model that fundamentally reshapes the edge set rather than merely rescaling weights. Inter-method edge overlap is low (Jaccard index 0.08-0.45) and activity dependent, demonstrating that these methods identify substantially different connections from identical data. Our results reveal that methodological choices including construction method, threshold, and normalization introduce hidden degrees of freedom with effect sizes comparable to the biological signals being measured. We provide practical recommendations for parameter selection, reporting, and cross-method validation in MEA-based network neuroscience. Author SummaryWhen we record electrical activity from brain tissue using grids of electrodes, we can ask how different sites influence one another and map the tissue as a network of connections. Thanks to novel culturing methods, this approach is increasingly used to study human brain slices. However, deciding what is "connected" is not well defined. Researchers use several different methods, and it has never been clear how much this choice shapes the network they end up describing. Here we compared three widely used methods on 37 recordings from human cortical slices spanning a range of activity levels. We found that the method alone can change the apparent structure of the network as much as real biological differences do. The methods frequently disagreed about which connections exist and some technical choices, including normalization techniques, had surprisingly large effects. Because these hidden choices can rival the biological signal, we provide this benchmarking work with practical recommendations for selecting, reporting, and cross-checking methods, so that network studies of brain tissue become more transparent, comparable, and reproducible.

Over the last decade, pre-clinical research has witnessed the advancement of human induced pluripotent stem-cell derived 3D brain organoid models and their differentiation into specific brain regions. In the realm of Parkinsons disease research, development of midbrain-specific organoids has enabled studies of this neurodegenerative disorder in patient derived 3D organoid models that attempt to recapitulate the human brain complexity. Within this line of research, neural functionality of the organoid models is established through electrophysiology. As a novel methodological approach, this study aimed to establish whether clustering of electrophysiological activity originating from midbrain organoids would aid in identifying different types of action-potential waveforms exhibited by neurons within the organoid model. Long-term extracellular electrophysiological recordings were conducted by use of a multi-electrode array device. The local field potential signal was spike-sorted, and the extracted putative neuron units were clustered into groups of spike waveform profiles. After establishing this methodological analysis pipeline, the clusters of waveform types were further analyzed in terms of electrophysiology. Results revealed that the clustering approach was successful at identifying three types of spike waveforms categories. Furthermore, it was proposed that one spike waveform profile potentially originated from dopaminergic neurons, which were one on the neural cells populating the organoid models. Overall, this study has successfully established a new methodological clustering approach to analyze electrophysiological data recoded from 3D organoid models in the context of Parkinsons disease modelling and organoid model development research.

As multisite and multi-study data aggregation becomes increasingly common for improving statistical power and sample diversity, robust harmonization methods are needed to address biases introduced by batch variation, particularly in neuroimaging research. Although a variety of harmonization approaches are available, the lack of systematic guidance for diagnosing batch effects and selecting appropriate methods remains a major challenge. To address this gap, we introduce ComBatFamQC, a comprehensive R package designed to streamline batch-effect diagnosis, harmonization, and post-harmonization analysis. ComBatFamQC integrates a user-friendly Shiny app for interactive batch-effect diagnostics, state-of-the-art harmonization methods from the ComBat family, including ComBat, longitudinal ComBat, ComBat-GAM, and CovBat, and tools for downstream analysis after harmonization. The package provides qualitative visualizations, statistical tests for batch-effect assessment, and a consistent interface that supports both in-sample and out-of-sample harmonization through the Shiny app, the R console, or the command line. In addition, it includes functions for post-harmonization analyses to facilitate downstream modeling. Its modular design also supports the systematic incorporation of future harmonization methods and expanded downstream analysis capabilities.

GapAutomated spike sorting algorithms have revolutionized the way neuronal activity is extracted from extracellular recordings, yet they remain imperfect. Specifically, inaccurate acceptance of noise-based units not only leaves researchers with clusters that require extensive manual curation, an essential but time-consuming process, that also leads to significant subjectivity in the selection of units. In an era of high-density probes like Neuropixels, where an hour of data can exceed 80 GB, manual curation is no longer scalable, automation of standard criteria can speed data curation and ensure quality of datasets. Here, we developed a semi-automated curation pipeline to label the quality of units after automated curation by Kilosort. ApproachOur algorithm standardizes criteria for labeling of Noise, Multi-Unit Activity (MUA), and Good Units using a combination of spike rate, spike timing metrics (from autocorrelogram), and waveform-based physiological features such as peak amplitude, slopes, half-width, and inter-channel correlation. Based on these features, clusters are assigned standardized labels (good, noise, multi-unit activity) that can be imported directly into Phy, where they serve as curation aids rather than absolute classifications, supporting but not replacing expert judgment. Heuristically, "noise" units are those unlikely to be neuronal in origin; "MUA" includes units with significant neural contribution (i.e., neuronal waveform) but with some degree of clear imperfection to be further cleaned, and "good" units are those without any clear deviation from ideal unit criteria. By ensuring accurate selection of acceptable units, we enable robust downstream analyses such as neural decoding and longitudinal tracking of neuron identity. Thresholds for all metrics were chosen to maximize the matching of algorithm output to that of 2 expert manual curators. Of note, users may alter thresholds either based on their own judgment or using an included tool to semi-automatically find thresholds that optimize SpikeCleaner with their own expert curation. Results: To benchmark, we compared the outputs of our algorithm to expert-labels curated in Phy by two expert users across three recordings. SpikeCleaner achieved an average of 97% accuracy vs. experts & 92% F1 score in classifying Single Units. It achieved an accuracy of 97% & 92% F1 score in full-category agreement (SU, MUA, Noise), and 97% accuracy & 95% F1 score in distinguishing Neuronal vs. Non-Neuronal units.