U-fiber analysis: a toolbox for automated quantification of U-fibers and white matter hyperintensities

Background

Whether white matter hyperintensities (WMHs) involve U-fibers is of great value in understanding the different etiologies of cerebral white matter (WM) lesions. However, clinical practice currently relies only on the naked eye to determine whether WMHs are in the vicinity of U-fibers, and there is a lack of good neuroimaging tools to quantify WMHs and U-fibers.

Methods

Here, we developed a multimodal neuroimaging toolbox named U-fiber analysis (UFA) that can automatically extract WMHs and quantitatively characterize the volume and number of WMHs in different brain regions. In addition, we proposed an anatomically constrained U-fiber tracking scheme and quantitatively characterized the microstructure diffusion properties, fiber length, and number of U-fibers in different brain regions to help clinicians to quantitatively determine whether WMHs in the proximal cortex disrupt the microstructure of U-fibers. To validate the utility of the UFA toolbox, we analyzed the neuroimaging data from 246 patients with cerebral small vessel disease (cSVD) enrolled at Zhongshan Hospital between March 2018 and November 2019 in a cross-sectional study.

Results

According to the manual judgment of the clinician, the patients with cSVD were divided into a WMHs involved U-fiber group (U-fiber-involved group, 51 cases) and WMHs not involved U-fiber group (U-fiber-spared group, 163 cases). There were no significant differences between the U-fiber-spared group and the U-fiber-involved group in terms of age (P=0.143), gender (P=0.462), education (P=0.151), Mini-Mental State Examination (MMSE) scores (P=0.151), and Montreal Cognitive Assessment (MoCA) scores (P=0.411). However, patients in the U-fiber-involved group had higher Fazekas scores (P<0.001) and significantly higher whole brain WMHs (P=0.046) and deep WMH volumes (P<0.001) compared to patients in the U-fiber-spared group. Moreover, the U-fiber-involved group had higher WMH volumes in the bilateral frontal [P(left) <0.001, P(right) <0.001] and parietal lobes [P(left) <0.001, P(right) <0.001]. On the other hand, patients in the U-fiber-involved group had higher mean diffusivity (MD) and axial diffusivity (AD) in the bilateral parietal [P(left, MD) =0.048, P(right, MD) =0.045, P(left, AD) =0.015, P(right, AD) =0.015] and right frontal-parietal regions [P(MD) =0.048, P(AD) =0.027], and had significantly reduced mean fiber length and number in the right parietal [P(length) =0.013, P(number) =0.028] and right frontal-parietal regions [P(length) =0.048] compared to patients in the U-fiber-spared group.

Conclusions

Our results suggest that WMHs in the proximal cortex may disrupt the microstructure of U-fibers. Our tool may provide new insights into the understanding of WM lesions of different etiologies in the brain.

A MATLAB toolbox to fit and forecast growth trajectories using phenomenological growth models: Application to epidemic outbreaks

Background

Simple dynamic modeling tools can be useful for generating real-time short-term forecasts with quantified uncertainty of the trajectory of diverse growth processes unfolding in nature and society, including disease outbreaks. An easy-to-use and flexible toolbox for this purpose is lacking.

Results

In this tutorial-based primer, we introduce and illustrate a user-friendly MATLAB toolbox for fitting and forecasting time-series trajectories using phenomenological dynamic growth models based on ordinary differential equations. This toolbox is accessible to various audiences, including students training in time-series forecasting, dynamic growth modeling, parameter estimation, parameter uncertainty and identifiability, model comparison, performance metrics, and forecast evaluation, as well as researchers and policymakers who need to conduct short-term forecasts in real-time. The models included in the toolbox capture exponential and sub-exponential growth patterns that typically follow a rising pattern followed by a decline phase, a common feature of contagion processes. Models include the 2-parameter generalized-growth model, which has proved useful to characterize and forecast the ascending phase of epidemic outbreaks, and the Gompertz model as well as the 3-parameter generalized logistic-growth model and the Richards model, which have demonstrated competitive performance in forecasting single peak outbreaks.The toolbox provides a tutorial for forecasting time-series trajectories that include the full uncertainty distribution, derived through parametric bootstrapping, which is needed to construct prediction intervals and evaluate their accuracy. Functions are available to assess forecasting performance across different models, estimation methods, error structures in the data, and forecasting horizons. The toolbox also includes functions to quantify forecasting performance using metrics that evaluate point and distributional forecasts, including the weighted interval score.

Conclusions

We have developed the first comprehensive toolbox to characterize and forecast time-series data using simple phenomenological growth models. As a contagion process takes off, the tools presented in this tutorial can facilitate policymaking to guide the implementation of control strategies and assess the impact of interventions. The toolbox functionality is demonstrated through various examples, including a tutorial video, and is illustrated using weekly data on the monkeypox epidemic in the USA.

A MATLAB toolbox to fit and forecast growth trajectories using phenomenological growth models: Application to epidemic outbreaks

BACKGROUND

Simple dynamic modeling tools can be useful for generating real-time short-term forecasts with quantified uncertainty of the trajectory of diverse growth processes unfolding in nature and society, including disease outbreaks. An easy-to-use and flexible toolbox for this purpose is lacking.

RESULTS

In this tutorial-based primer, we introduce and illustrate a user-friendly MATLAB toolbox for fitting and forecasting time-series trajectories using phenomenological dynamic growth models based on ordinary differential equations. This toolbox is accessible to various audiences, including students training in time-series forecasting, dynamic growth modeling, parameter estimation, parameter uncertainty and identifiability, model comparison, performance metrics, and forecast evaluation, as well as researchers and policymakers who need to conduct short-term forecasts in real-time. The models included in the toolbox capture exponential and sub-exponential growth patterns that typically follow a rising pattern followed by a decline phase, a common feature of contagion processes. Models include the 2-parameter generalized-growth model, which has proved useful to characterize and forecast the ascending phase of epidemic outbreaks, and the Gompertz model as well as the 3-parameter generalized logistic-growth model and the Richards model, which have demonstrated competitive performance in forecasting single peak outbreaks.The toolbox provides a tutorial for forecasting time-series trajectories that include the full uncertainty distribution, derived through parametric bootstrapping, which is needed to construct prediction intervals and evaluate their accuracy. Functions are available to assess forecasting performance across different models, estimation methods, error structures in the data, and forecasting horizons. The toolbox also includes functions to quantify forecasting performance using metrics that evaluate point and distributional forecasts, including the weighted interval score.

CONCLUSIONS

We have developed the first comprehensive toolbox to characterize and forecast time-series data using simple phenomenological growth models. As a contagion process takes off, the tools presented in this tutorial can facilitate policymaking to guide the implementation of control strategies and assess the impact of interventions. The toolbox functionality is demonstrated through various examples, including a tutorial video, and is illustrated using weekly data on the monkeypox epidemic in the USA.

A user-friendly automatic toolbox for hand kinematic analysis, clinical assessment and postural synergies extraction

The clinical assessment of the human hand is typically conducted through questionnaires or tests that include objective (e.g., time) and subjective (e.g., grasp quality) outcome measures. However, there are other important indicators that should be considered to quantify grasp and movement quality in addition to the time needed by a subject to execute a task, and this is essential for human and artificial hands that attempt to replicate the human hand properties. The correct estimation of hand kinematics is fundamental for computing these indicators with high fidelity, and a technical background is typically required to perform this analysis. In addition, to understand human motor control strategies as well as to replicate them on artificial devices, postural synergies were widely explored in recent years. Synergies should be analyzed not only to investigate possible modifications due to musculoskeletal and/or neuromuscular disorders, but also to test biomimetic hands. The aim of this work is to present an open source toolbox to perform all-in-one kinematic analysis and clinical assessment of the hand, as well as to perform postural synergies extraction. In the example provided in this work, the tool takes as input the position of 28 retroreflective markers with a diameter of 6 mm, positioned on specific anatomical landmarks of the hand and recorded with an optoelectronic motion capture system, and automatically performs 1) hand kinematic analysis (i.e., computation of 23 joint angles); 2) clinical assessment, by computing indicators that allow quantifying movement efficiency (Peak Grip Aperture), smoothness (Normalized Dimensionless Jerk Grasp Aperture) and speed (Peak Velocity of Grasp Aperture), planning capabilities (Time to Peak Grip Aperture), spatial posture (Wrist and Finger Joint Angles) and grasp stability (Posture of Hand Finger Joints), and 3) postural synergies extraction and analysis through the Pareto, Scree and Loadings plots. Two examples are described to demonstrate the applicability of the toolbox: the first one aiming at performing a clinical assessment of a volunteer and the second one aiming at extracting and analyzing the volunteer's postural synergies. The tool allows calculating joint angles with high accuracy (reconstruction errors below 4 mm and 3.2 mm for the fingers and wrist respectively) and automatically performing clinical assessment and postural synergies extraction. Results can be visually inspected, and data can be saved for any desired post processing analysis. Custom-made protocols to extract joint angles, based on different markersets, could be also integrated in the toolbox. The tool can be easily exploitable in clinical contexts, as it does not require any particular technical knowledge to be used, as confirmed by the usability evaluation conducted (perceived usability = 94.2 ± 5.4). In addition, it can be integrated with the SynGrasp toolbox to perform grasp analysis of underactuated virtual hands based on postural synergies.

Software for non-parametric image registration of 2-photon imaging data

Functional 2-photon microscopy is a key technology for imaging neuronal activity. The recorded image sequences, however, can contain non-rigid movement artifacts which requires high-accuracy movement correction. Variational optical flow (OF) estimation is a group of methods for motion analysis with established performance in many computer vision areas. However, it has yet to be adapted to the statistics of 2-photon neuroimaging data. In this work, we present the motion compensation method Flow-Registration that outperforms previous alignment tools and allows to align and reconstruct even low signal-to-noise ratio 2-photon imaging data and is able to compensate high-divergence displacements during local drug injections. The method is based on statistics of such data and integrates previous advances in variational OF estimation. Our method is available as an easy-to-use ImageJ/FIJI plugin as well as a MATLAB toolbox with modular, object oriented file IO, native multi-channel support and compatibility with existing 2-photon imaging suites.

NIRS-ICA: A MATLAB Toolbox for Independent Component Analysis Applied in fNIRS Studies

Independent component analysis (ICA) is a multivariate approach that has been widely used in analyzing brain imaging data. In the field of functional near-infrared spectroscopy (fNIRS), its promising effectiveness has been shown in both removing noise and extracting neuronal activity-related sources. The application of ICA remains challenging due to its complexity in usage, and an easy-to-use toolbox dedicated to ICA processing is still lacking in the fNIRS community. In this study, we propose NIRS-ICA, an open-source MATLAB toolbox to ease the difficulty of ICA application for fNIRS studies. NIRS-ICA incorporates commonly used ICA algorithms for source separation, user-friendly GUI, and quantitative evaluation metrics assisting source selection, which facilitate both removing noise and extracting neuronal activity-related sources. The options used in the processing can also be reported easily, which promotes using ICA in a more reproducible way. The proposed toolbox is validated and demonstrated based on both simulative and real fNIRS datasets. We expect the release of the toolbox will extent the application for ICA in the fNIRS community.

NICA: A Novel Toolbox for Near-Infrared Spectroscopy Calculations and Analyses

Functional near-infrared spectroscopy (fNIRS) measures the functional activity of the cerebral cortex. The concentration changes of oxygenated (oxy-Hb) and deoxygenated hemoglobin (deoxy-Hb) can be detected and associated with activation of the cortex in the investigated area (neurovascular coupling). Recorded signals of hemodynamic responses may contain influences from physiological signals (systemic influences, physiological artifacts) which do not originate from the cerebral cortex activity. The physiological artifacts contain the blood pressure (BP), respiratory patterns, and the pulsation of the heart. In order to perform a comprehensive analysis of recorded fNIRS data, a proper correction of these physiological artifacts is necessary. This article introduces NICA – a novel toolbox for near-infrared spectroscopy calculations and analyses based on MATLAB. With NICA it is possible to process and visualize fNIRS data, including different signal processing methods for physiological artifact correction. The artifact correction methods used in this toolbox are common average reference (CAR), independent component analysis (ICA), and transfer function (TF) models. A practical example provides results from a study, where NICA was used for analyzing the measurement data, in order to demonstrate the signal processing steps and the physiological artifact correction. The toolbox was developed for fNIRS data recorded with the NIRScout 1624 measurement device and the corresponding recording software NIRStar.

Network-based Observability and Controllability Analysis of Dynamical Systems: the NOCAD toolbox

The network science-based determination of driver nodes and sensor placement has become increasingly popular in the field of dynamical systems over the last decade. In this paper, the applicability of the methodology in the field of life sciences is introduced through the analysis of the neural network of Caenorhabditis elegans. Simultaneously, an Octave and MATLAB-compatible NOCAD toolbox is proposed that provides a set of methods to automatically generate the relevant structural controllability and observability associated measures for linear or linearised systems and compare the different sensor placement methods.

A MATLAB toolbox for multivariate analysis of brain networks

Complex brain networks formed via structural and functional interactions among brain regions are believed to underlie information processing and cognitive function. A growing number of studies indicate that altered brain network topology is associated with physiological, behavioral, and cognitive abnormalities. Graph theory is showing promise as a method for evaluating and explaining brain networks. However, multivariate frameworks that provide statistical inferences about how such networks relate to covariates of interest, such as disease phenotypes, in different study populations are yet to be developed. We have developed a freely available MATLAB toolbox with a graphical user interface that bridges this important gap between brain network analyses and statistical inference. The modeling framework implemented in this toolbox utilizes a mixed-effects multivariate regression framework that allows assessing brain network differences between study populations as well as assessing the effects of covariates of interest such as age, disease phenotype, and risk factors on the density and strength of brain connections in global (i.e., whole-brain) and local (i.e., subnetworks) brain networks. Confounding variables, such as sex, are controlled for through the implemented framework. A variety of neuroimaging data such as fMRI, EEG, and DTI can be analyzed with this toolbox, which makes it useful for a wide range of studies examining the structure and function of brain networks. The toolbox uses SAS, R, or Python (depending on software availability) to perform the statistical modeling. We also provide a clustering-based data reduction method that helps with model convergence and substantially reduces modeling time for large data sets.

Neuroscience Information Toolbox: An Open Source Toolbox for EEG-fMRI Multimodal Fusion Analysis

Recently, scalp electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) multimodal fusion has been pursued in an effort to study human brain function and dysfunction to obtain more comprehensive information on brain activity in which the spatial and temporal resolutions are both satisfactory. However, a more flexible and easy-to-use toolbox for EEG–fMRI multimodal fusion is still lacking. In this study, we therefore developed a freely available and open-source MATLAB graphical user interface toolbox, known as the Neuroscience Information Toolbox (NIT), for EEG–fMRI multimodal fusion analysis. The NIT consists of three modules: (1) the fMRI module, which has batch fMRI preprocessing, nuisance signal removal, bandpass filtering, and calculation of resting-state measures; (2) the EEG module, which includes artifact removal, extracting EEG features (event onset, power, and amplitude), and marking interesting events; and (3) the fusion module, in which fMRI-informed EEG analysis and EEG-informed fMRI analysis are included. The NIT was designed to provide a convenient and easy-to-use toolbox for researchers, especially for novice users. The NIT can be downloaded for free at http://www.neuro.uestc.edu.cn/NIT.html, and detailed information, including the introduction of NIT, user’s manual and example data sets, can also be observed on this website. We hope that the NIT is a promising toolbox for exploring brain information in various EEG and fMRI studies.

MATLAB Toolboxes for Reference Electrode Standardization Technique (REST) of Scalp EEG

Reference electrode standardization technique (REST) has been increasingly acknowledged and applied as a re-reference technique to transform an actual multi-channels recordings to approximately zero reference ones in electroencephalography/event-related potentials (EEG/ERPs) community around the world in recent years. However, a more easy-to-use toolbox for re-referencing scalp EEG data to zero reference is still lacking. Here, we have therefore developed two open-source MATLAB toolboxes for REST of scalp EEG. One version of REST is closely integrated into EEGLAB, which is a popular MATLAB toolbox for processing the EEG data; and another is a batch version to make it more convenient and efficient for experienced users. Both of them are designed to provide an easy-to-use for novice researchers and flexibility for experienced researchers. All versions of the REST toolboxes can be freely downloaded at http://www.neuro.uestc.edu.cn/rest/Down.html, and the detailed information including publications, comments and documents on REST can also be found from this website. An example of usage is given with comparative results of REST and average reference. We hope these user-friendly REST toolboxes could make the relatively novel technique of REST easier to study, especially for applications in various EEG studies.

COMETS2: An advanced MATLAB toolbox for the numerical analysis of electric fields generated by transcranial direct current stimulation

BACKGROUND

Since there is no way to measure electric current generated by transcranial direct current stimulation (tDCS) inside the human head through in vivo experiments, numerical analysis based on the finite element method has been widely used to estimate the electric field inside the head. In 2013, we released a MATLAB toolbox named COMETS, which has been used by a number of groups and has helped researchers to gain insight into the electric field distribution during stimulation. The aim of this study was to develop an advanced MATLAB toolbox, named COMETS2, for the numerical analysis of the electric field generated by tDCS.

NEW METHOD

COMETS2 can generate any sizes of rectangular pad electrodes on any positions on the scalp surface. To reduce the large computational burden when repeatedly testing multiple electrode locations and sizes, a new technique to decompose the global stiffness matrix was proposed.

RESULTS

As examples of potential applications, we observed the effects of sizes and displacements of electrodes on the results of electric field analysis. The proposed mesh decomposition method significantly enhanced the overall computational efficiency.

COMPARISON WITH EXISTING METHODS

We implemented an automatic electrode modeler for the first time, and proposed a new technique to enhance the computational efficiency.

CONCLUSIONS

In this paper, an efficient toolbox for tDCS analysis is introduced (freely available at http://www.cometstool.com). It is expected that COMETS2 will be a useful toolbox for researchers who want to benefit from the numerical analysis of electric fields generated by tDCS.

TMSEEG: A MATLAB-Based Graphical User Interface for Processing Electrophysiological Signals during Transcranial Magnetic Stimulation

Concurrent recording of electroencephalography (EEG) during transcranial magnetic stimulation (TMS) is an emerging and powerful tool for studying brain health and function. Despite a growing interest in adaptation of TMS-EEG across neuroscience disciplines, its widespread utility is limited by signal processing challenges. These challenges arise due to the nature of TMS and the sensitivity of EEG to artifacts that often mask TMS-evoked potentials (TEP)s. With an increase in the complexity of data processing methods and a growing interest in multi-site data integration, analysis of TMS-EEG data requires the development of a standardized method to recover TEPs from various sources of artifacts. This article introduces TMSEEG, an open-source MATLAB application comprised of multiple algorithms organized to facilitate a step-by-step procedure for TMS-EEG signal processing. Using a modular design and interactive graphical user interface (GUI), this toolbox aims to streamline TMS-EEG signal processing for both novice and experienced users. Specifically, TMSEEG provides: (i) targeted removal of TMS-induced and general EEG artifacts; (ii) a step-by-step modular workflow with flexibility to modify existing algorithms and add customized algorithms; (iii) a comprehensive display and quantification of artifacts; (iv) quality control check points with visual feedback of TEPs throughout the data processing workflow; and (v) capability to label and store a database of artifacts. In addition to these features, the software architecture of TMSEEG ensures minimal user effort in initial setup and configuration of parameters for each processing step. This is partly accomplished through a close integration with EEGLAB, a widely used open-source toolbox for EEG signal processing. In this article, we introduce TMSEEG, validate its features and demonstrate its application in extracting TEPs across several single- and multi-pulse TMS protocols. As the first open-source GUI-based pipeline for TMS-EEG signal processing, this toolbox intends to promote the widespread utility and standardization of an emerging technology in brain research.

moviEEG: An animation toolbox for visualization of intracranial electroencephalography synchronization dynamics

OBJECTIVE

We introduce and describe the functions of moviEEG (Multiple Overlay Visualizations for Intracranial ElectroEncephaloGraphy), a novel MATLAB-based toolbox for spatiotemporal mapping of network synchronization dynamics in intracranial electroencephalography (iEEG) data.

METHODS

The toolbox integrates visualizations of inter-electrode phase-locking relationships in peri-ictal epileptogenic networks with signal spectral properties and graph-theoretical network measures overlaid upon operating room images of the electrode grid. Functional connectivity between every electrode pair is evaluated over a sliding window indexed by phase synchrony.

RESULTS

Two case studies are presented to provide preliminary evidence for the application of the toolbox to guide network-based mapping of epileptogenic cortex and to distinguish these regions from eloquent brain networks. In both cases, epileptogenic cortex was visually distinct.

CONCLUSION

We introduce moviEEG, a novel toolbox for animation of oscillatory network dynamics in iEEG data, and provide two case studies showing preliminary evidence for utility of the toolbox in delineating the epileptogenic zone.

SIGNIFICANCE

Despite evidence that atypical network synchronization has shown to be altered in epileptogenic brain regions, network based techniques have yet to be incorporated into clinical pre-surgical mapping. moviEEG provides a set of functions to enable easy visualization with network based techniques.