A High-Resolution LED Stimulator for Steady-State Visual Stimulation: Customizable, Affordable, and Open Source.

Visually evoked steady-state potentials (SSVEPs) are neural responses elicited by visual stimuli oscillating at specific frequencies. In this study, we introduce a novel LED stimulator system explicitly designed for steady-state visual stimulation, offering precise control over visual stimulus parameters, including frequency resolution, luminance, and the ability to control the phase at the end of the stimulation. The LED stimulator provides a personalized, modular, and affordable option for experimental setups. Based on the Teensy 3.2 board, the stimulator utilizes direct digital synthesis and pulse width modulation techniques to control the LEDs. We validated its performance through four experiments: the first two measured LED light intensities directly, while the last two assessed the stimulator's impact on EEG recordings. The results demonstrate that the stimulator can deliver a stimulus suitable for generating SSVEPs with the desired frequency and phase resolution. As an open source resource, we provide comprehensive documentation, including all necessary codes and electrical diagrams, which facilitates the system's replication and adaptation for specific experimental requirements, enhancing its potential for widespread use in the field of neuroscience setups.

Effects of Tau and Sampling Frequency on the Regularity Analysis of ECG and EEG Signals Using ApEn and SampEn Entropy Estimators.

Electrocardiography (ECG) and electroencephalography (EEG) signals provide clinical information relevant to determine a patient's health status. The nonlinear analysis of ECG and EEG signals allows for discovering characteristics that could not be found with traditional methods based on amplitude and frequency. Approximate entropy (ApEn) and sampling entropy (SampEn) are nonlinear data analysis algorithms that measure the data's regularity, and these are used to classify different electrophysiological signals as normal or pathological. Entropy calculation requires setting the parameters r (tolerance threshold), m (immersion dimension), and τ (time delay), with the last one being related to how the time series is downsampled. In this study, we showed the dependence of ApEn and SampEn on different values of τ, for ECG and EEG signals with different sampling frequencies (Fs), extracted from a digital repository. We considered four values of Fs (128, 256, 384, and 512 Hz for the ECG signals, and 160, 320, 480, and 640 Hz for the EEG signals) and five values of τ (from 1 to 5). We performed parametric and nonparametric statistical tests to confirm that the groups of normal and pathological ECG and EEG signals were significantly different (p < 0.05) for each F and τ value. The separation between the entropy values of regular and irregular signals was variable, demonstrating the dependence of ApEn and SampEn with Fs and τ. For ECG signals, the separation between the conditions was more robust when using SampEn, the lowest value of Fs, and τ larger than 1. For EEG signals, the separation between the conditions was more robust when using SampEn with large values of Fs and τ larger than 1. Therefore, adjusting τ may be convenient for signals that were acquired with different Fs to ensure a reliable clinical classification. Furthermore, it is useful to set τ to values larger than 1 to reduce the computational cost.

Elevating understanding: Linking high-altitude hypoxia to brain aging through EEG functional connectivity and spectral analyses.

High-altitude hypoxia triggers brain function changes reminiscent of those in healthy aging and Alzheimer's disease, compromising cognition and executive functions. Our study sought to validate high-altitude hypoxia as a model for assessing brain activity disruptions akin to aging. We collected EEG data from 16 healthy volunteers during acute high-altitude hypoxia (at 4,000 masl) and at sea level, focusing on relative changes in power and aperiodic slope of the EEG spectrum due to hypoxia. Additionally, we examined functional connectivity using wPLI, and functional segregation and integration using graph theory tools. High altitude led to slower brain oscillations, that is, increased δ and reduced α power, and flattened the 1/f aperiodic slope, indicating higher electrophysiological noise, akin to healthy aging. Notably, functional integration strengthened in the θ band, exhibiting unique topographical patterns at the subnetwork level, including increased frontocentral and reduced occipitoparietal integration. Moreover, we discovered significant correlations between subjects' age, 1/f slope, θ band integration, and observed robust effects of hypoxia after adjusting for age. Our findings shed light on how reduced oxygen levels at high altitudes influence brain activity patterns resembling those in neurodegenerative disorders and aging, making high-altitude hypoxia a promising model for comprehending the brain in health and disease.

Exposure to high-altitude hypoxia, with reduced oxygen levels, can replicate brain function changes akin to aging and Alzheimer's disease. In our work, we propose high-altitude hypoxia as a possible reversible model of human brain aging. We gathered EEG data at high altitude and sea level, investigating the impact of hypoxia on brainwave patterns and connectivity. Our findings revealed that high-altitude exposure led to slower and noisier brain oscillations and produced altered brain connectivity, resembling some remarkable changes seen in the aging process. Intriguingly, these changes were linked to age, even when hypoxia's effects were considered. Our research unveils how high-altitude conditions emulate brain patterns associated with aging and neurodegenerative conditions, providing valuable insights into the understanding of both normal and impaired brain function.

Temporal differences in the myoelectric activity of lower limb muscles during rearfoot and forefoot running: A statistical parametric mapping approach.

Forefoot (FF) and rearfoot (RF) running techniques can induce different lower-limb muscle activation patterns. However, few studies have evaluated temporal changes in the electromyographic activity (EMG) of lower limb muscles during running. The aim of this study was to compare temporal changes in EMG amplitude between RF and FF running techniques. Eleven recreational runners ran on a treadmill at a self-selected speed, once using a RF strike pattern and once using a FF strike pattern (randomized order). The EMG of five lower limb muscles [rectus femoris (RFe), biceps femoris (BF), tibialis anterior (TA), medial and lateral gastrocnemius (MG and LG)] was evaluated, using bipolar electrodes. EMG data from the RF and FF running techniques was then processed and compared with statistical parametric mapping (SPM), dividing the analysis of the running cycle into stance and swing phases. The MG and LG muscles showed higher activation during FF running at the beginning of the stance phase and at the end of the swing phase. During the end of the swing phase, the TA muscle's EMG amplitude was higher, when the RF running technique was used. A higher level of co-activation between the gastrocnemius and TA muscles was observed in both stance and swing phases using RF. The myoelectric behaviour of the RFe and BF muscles was similar during both running techniques. The current findings highlight that the two running techniques predominately reflect adjustments of the shank and not the thigh muscles, in both phases of the running cycle.HighlightsStatistical parametric mapping (SPM) can reveal temporal differences in muscle activity between running techniques.The medial and lateral gastrocnemius muscles were more active at specific time-instants of the initial stance and late swing phases during forefoot (FF) running compared to rearfoot (RF) running.Higher activation was observed for the tibialis anterior muscle at the end of the swing phase during RF runningContrary to the muscle activity differences observed in the leg muscles, the muscle activity of the thigh muscles was similar during RF and FF running.

DIVA Meets EEG: Model Validation Using Formant-Shift Reflex.

The neurocomputational model 'Directions into Velocities of Articulators' (DIVA) was developed to account for various aspects of normal and disordered speech production and acquisition. The neural substrates of DIVA were established through functional magnetic resonance imaging (fMRI), providing physiological validation of the model. This study introduces DIVA_EEG an extension of DIVA that utilizes electroencephalography (EEG) to leverage the high temporal resolution and broad availability of EEG over fMRI. For the development of DIVA_EEG, EEG-like signals were derived from original equations describing the activity of the different DIVA maps. Synthetic EEG associated with the utterance of syllables was generated when both unperturbed and perturbed auditory feedback (first formant perturbations) were simulated. The cortical activation maps derived from synthetic EEG closely resembled those of the original DIVA model. To validate DIVA_EEG, the EEG of individuals with typical voices (N = 30) was acquired during an altered auditory feedback paradigm. The resulting empirical brain activity maps significantly overlapped with those predicted by DIVA_EEG. In conjunction with other recent model extensions, DIVA_EEG lays the foundations for constructing a complete neurocomputational framework to tackle vocal and speech disorders, which can guide model-driven personalized interventions.

Biased instantaneous regional muscle activation maps: Embedded fuzzy topology and image feature analysis.

The instantaneous spatial representation of electrical propagation produced by muscle contraction may introduce bias in surface electromyographical (sEMG) activation maps. Here, we described the effect of instantaneous spatial representation (sEMG segmentation) on embedded fuzzy topological polyhedrons and image features extracted from sEMG activation maps. We analyzed 73,008 topographic sEMG activation maps from seven healthy participants (age 21.4 ± 1.5 years and body mass 74.5 ± 8.5 kg) who performed submaximal isometric plantar flexions with 64 surface electrodes placed over the medial gastrocnemius muscle. Window lengths of 50, 100, 150, 250, 500, and 1,000 ms and overlap of 0, 25, 50, 75, and 90% to change sEMG map generation were tested in a factorial design (grid search). The Shannon entropy and volume of global embedded tri-dimensional geometries (polyhedron projections), and the Shannon entropy, location of the center (LoC), and image moments of maps were analyzed. The polyhedron volume increased when the overlap was <25% and >75%. Entropy decreased when the overlap was <25% and >75% and when the window length was <100 ms and >500 ms. The LoC in the x-axis, entropy, and the histogram moments of maps showed effects for overlap (p < 0.001), while the LoC in the y-axis and entropy showed effects for both overlap and window length (p < 0.001). In conclusion, the instantaneous sEMG maps are first affected by outer parameters of the overlap, followed by the length of the window. Thus, choosing the window length and overlap parameters can introduce bias in sEMG activation maps, resulting in distorted regional muscle activation.

Kalman Filter Implementation of Subglottal Impedance-Based Inverse Filtering to Estimate Glottal Airflow during Phonation.

Subglottal Impedance-Based Inverse Filtering (IBIF) allows for the continuous, non-invasive estimation of glottal airflow from a surface accelerometer placed over the anterior neck skin below the larynx. It has been shown to be advantageous for the ambulatory monitoring of vocal function, specifically in the use of high-order statistics to understand long-term vocal behavior. However, during long-term ambulatory recordings over several days, conditions may drift from the laboratory environment where the IBIF parameters were initially estimated due to sensor positioning, skin attachment, or temperature, among other factors. Observation uncertainties and model mismatch may result in significant deviations in the glottal airflow estimates; unfortunately, they are very difficult to quantify in ambulatory conditions due to a lack of a reference signal. To address this issue, we propose a Kalman filter implementation of the IBIF filter, which allows for both estimating the model uncertainty and adapting the airflow estimates to correct for signal deviations. One-way analysis of variance (ANOVA) results from laboratory experiments using the Rainbow Passage indicate an improvement using the modified Kalman filter on amplitude-based measures for phonotraumatic vocal hyperfunction (PVH) subjects compared to the standard IBIF; the latter showing a statistically difference (p-value = 0.02, F = 4.1) with respect to a reference glottal volume velocity signal estimated from a single notch filter used here as ground-truth in this work. In contrast, maximum flow declination rates from subjects with vocal phonotrauma exhibit a small but statistically difference between the ground-truth signal and the modified Kalman filter when using one-way ANOVA (p-value = 0.04, F = 3.3). Other measures did not have significant differences with either the modified Kalman filter or IBIF compared to ground-truth, with the exception of H1-H2, whose performance deteriorates for both methods. Overall, both methods (modified Kalman filter and IBIF) show similar glottal airflow measures, with the advantage of the modified Kalman filter to improve amplitude estimation. Moreover, Kalman filter deviations from the IBIF output airflow might suggest a better representation of some fine details in the ground-truth glottal airflow signal. Other applications may take more advantage from the adaptation offered by the modified Kalman filter implementation.

The Transcranial Light Therapy Improves Synaptic Plasticity in the Alzheimer's Disease Mouse Model.

Alzheimer's disease (AD) is the main cause of dementia worldwide. Emerging non-invasive treatments such as photobiomodulation target the mitochondria to minimize brain damage, improving cognitive functions. In this work, an experimental design was carried out to evaluate the effect of transcranial light therapy (TLTC) on synaptic plasticity (SP) and cognitive functions in an AD animal model. Twenty-three mice were separated into two general groups: an APP/PS1 (ALZ) transgenic group and a wild-type (WT) group. Each group was randomly subdivided into two subgroups: mice with and without TLTC, depending on whether they would undergo treatment with TLTC. Cognitive function, measured through an object recognition task, showed non-significant improvement after TLTC. SP, on the other hand, was evaluated using four electrophysiological parameters from the Schaffer-CA1 collateral hippocampal synapses: excitatory field potentials (fEPSP), paired pulse facilitation (PPF), long-term depression (LTD), and long-term potentiation (LTP). An improvement was observed in subjects treated with TLTC, showing higher levels of LTP than those transgenic mice that were not exposed to the treatment. Therefore, the results obtained in this work showed that TLTC could be an efficient non-invasive treatment for AD-associated SP deficits.

Understanding the effect of window length and overlap for assessing sEMG in dynamic fatiguing contractions: A non-linear dimensionality reduction and clustering.

The Short-Time Fourier transform (STFT) is a helpful tool to identify muscle fatigue with clinical and sports applications. However, the choice of STFT parameters may affect the estimation of myoelectrical manifestations of fatigue. Here, we determine the effect of window length and overlap selections on the frequency slope and the coefficient of variation from EMG spectrum features in fatiguing contractions. We also determine whether STFT parameters affect the relationship between frequency slopes and task failure. Eighty-eight healthy adult men performed one-leg heel-rise until exhaustion. A factorial design with a window length of 50, 100, 250, 500, and 1000 ms with 0, 25, 50, 75, and 90% of overlap was used. The frequency slope was non-linearly fitted as a task failure function, followed by a dimensionality reduction and clustering analysis. The STFT parameters elicited five patterns. A small window length produced a higher slope frequency for the peak frequency (p < 0.001). The contrary was found for the mean and median frequency (p < 0.001). A larger window length elicited a higher slope frequency for the mean and peak frequencies. The largest frequency slope and dispersion was found for a window length of 50 ms without overlap using peak frequency. A combination of 250 ms with 50% of overlap reduced the dispersion both for peak, median, and mean frequency, but decreased the slope frequency. Therefore, the selection of STFT parameters during dynamic contractions should be accompanied by a mechanical measure of the task failure, and its parameters should be adjusted according to the experiment's requirements.

Persistence of EEG Alpha Entrainment Depends on Stimulus Phase at Offset.

Neural entrainment is the synchronization of neural activity to the frequency of repetitive external stimuli, which can be observed as an increase in the electroencephalogram (EEG) power spectrum at the driving frequency, -also known as the steady-state response. Although it has been systematically reported that the entrained EEG oscillation persists for approximately three cycles after stimulus offset, the neural mechanisms underpinning it remain unknown. Focusing on alpha oscillations, we adopt the dynamical excitation/inhibition framework, which suggests that phases of entrained EEG signals correspond to alternating excitatory/inhibitory states of the neural circuitry. We hypothesize that the duration of the persistence of entrainment is determined by the specific functional state of the entrained neural network at the time the stimulus ends. Steady-state visually evoked potentials (SSVEP) were elicited in 19 healthy volunteers at the participants' individual alpha peaks. Visual stimulation consisted of a sinusoidally-varying light terminating at one of four phases: 0, π/2, π, and 3π/2. The persistence duration of the oscillatory activity was analyzed as a function of the terminating phase of the stimulus. Phases of the SSVEP at the stimulus termination were distributed within a constant range of values relative to the phase of the stimulus. Longer persistence durations were obtained when visual stimulation terminated towards the troughs of the alpha oscillations, while shorter persistence durations occurred when stimuli terminated near the peaks. Source localization analysis suggests that the persistence of entrainment reflects the functioning of fronto-occipital neuronal circuits, which might prime the sensory representation of incoming visual stimuli based on predictions about stimulus rhythmicity. Consequently, different states of the network at the end of the stimulation, corresponding to different states of intrinsic neuronal coupling, may determine the time windows over which coding of incoming sensory stimulation is modulated by the preceding oscillatory activity.

Estimation of auditory steady-state responses based on the averaging of independent EEG epochs.

The amplitude of auditory steady-state responses (ASSRs) generated in the brainstem of rats exponentially decreases over the sequential averaging of EEG epochs. This behavior is partially due to the adaptation of the ASSR induced by the continuous and monotonous stimulation. In this study, we analyzed the potential clinical relevance of the ASSR adaptation. ASSR were elicited in eight anesthetized adult rats by 8-kHz tones, modulated in amplitude at 115 Hz. We called independent epochs to those EEG epochs acquired with sufficiently long inter-stimulus interval, so the ASSR contained in any given epoch is not affected by the previous stimulation. We tested whether the detection of ASSRs is improved when the response is computed by averaging independent EEG epochs, containing only unadapted auditory responses. The improvements in the ASSR detection obtained with standard, weighted and sorted averaging were compared. In the absence of artifacts, when the ASSR was elicited by continuous acoustic stimulation, the computation of the ASSR amplitude relied upon the averaging method. While the adaptive behavior of the ASSR was still evident after the weighting of epochs, the sorted averaging resulted in under-estimations of the ASSR amplitude. In the absence of artifacts, the ASSR amplitudes computed by averaging independent epochs did not depend on the averaging procedure. Averaging independent epochs resulted in higher ASSR amplitudes and halved the number of EEG epochs needed to be acquired to achieve the maximum detection rate of the ASSR. Acquisition protocols based on averaging independent EEG epochs, in combination with appropriate averaging methods for artifact reduction might contribute to develop more accurate hearing assessments based on ASSRs.

Ant Colony Clustering for ROI Identification in Functional Magnetic Resonance Imaging.

Brain network analysis using functional magnetic resonance imaging (fMRI) is a widely used technique. The first step of brain network analysis in fMRI is to detect regions of interest (ROIs). The signals from these ROIs are then used to evaluate neural networks and quantify neuronal dynamics. The two main methods to identify ROIs are based on brain atlas registration and clustering. This work proposes a bioinspired method that combines both paradigms. The method, dubbed HAnt, consists of an anatomical clustering of the signal followed by an ant clustering step. The method is evaluated empirically in both in silico and in vivo experiments. The results show a significantly better performance of the proposed approach compared to other brain parcellations obtained using purely clustering-based strategies or atlas-based parcellations.

Use of accelerometers for automatic regional chest movement recognition during tidal breathing in healthy subjects.

Recognition of breathing patterns helps clinicians to understand acute and chronic adaptations during exercise and pathological conditions. Wearable technologies combined with a proper data analysis provide a low cost option to monitor chest and abdominal wall movements. Here we set out to determine the feasibility of using accelerometry and machine learning to detect chest-abdominal wall movement patterns during tidal breathing. Furthermore, we determined the accelerometer positions included in the clusters, considering principal component domains. Eleven healthy participants (age: 21 ±â€¯0.2 y, BMI: 23.4 ±â€¯0.7 kg/m2, FEV1: 4.1 ±â€¯0.3 L, VO2: 4.6 ±â€¯0.2 mL/min kg) were included in this cross-sectional study. Spirometry and ergospirometry assessments were performed with participants seated with 13 accelerometers placed over the thorax. Data collection lasted 10  min. Following signal pre-processing, principal components and clustering analyses were performed. The Euclidean distances in respect to centroids were compared between the clusters (p < 0.05), identifying two clusters (p < 0.001). The first cluster included sensors located at the right and left second rib midline, body of sternum, left fourth rib midline, right and left second thoracic vertebra midline, and fifth thoracic vertebra. The second cluster included sensors at the fourth right rib midline, right and left seventh ribs, abdomen at linea alba, and right and left tenth thoracic vertebra midline. Costal-superior and costal-abdominal patterns were also recognized. We conclude that accelerometers placed on the chest and abdominal wall permit the identification of two clusters of movements regarding respiration biomechanics.

A Method for Tracking the Time Evolution of Steady-State Evoked Potentials.

Neural entrainment refers to the synchronization of neural activity to the periodicity of sensory stimuli. This synchronization defines the generation of steady-state evoked responses (i.e., oscillations in the electroencephalogram phase-locked to the driving stimuli). The classic interpretation of the amplitude of the steady-state evoked responses assumes a stereotypical time-invariant neural response plus random background fluctuations, such that averaging over repeated presentations of the stimulus recovers the stereotypical response. This approach ignores the dynamics of the steady-state, as in the case of the adaptation elicited by prolonged exposures to the stimulus. To analyze the dynamics of steady-state responses, it can be assumed that the time evolution of the response amplitude is the same in different stimulation runs separated by sufficiently long breaks. Based on this assumption, a method to characterize the time evolution of steady-state responses is presented. A sufficiently large number of recordings are acquired in response to the same experimental condition. Experimental runs (recordings) are column-wise averaged (i.e., runs are averaged but epoch within recordings are not averaged with the preceding segments). The column-wise averaging allows analysis of steady-state responses in recordings with remarkably high signal-to-noise ratios. Therefore, the averaged signal provides an accurate representation of the time evolution of the steady-state response, which can be analyzed in both the time and frequency domains. In this study, a detailed description of the method is provided, using steady-state visually evoked potentials as an example of a response. Advantages and caveats are evaluated based on a comparison with single-trial methods designed to analyze neural entrainment.

Towards a conceptual model for the use of home healthcare medical devices: The multi-parameter monitor case.

In the last decade there has been an increase in the use of medical devices in the home environment. These devices are commonly the same as those used in hospitals by healthcare professionals. The use of these these devices by lay users outside of a clinical environment may become unsafe. This study presents a methodology that allows decision makers to identify potential risk situations that may arise when lay users operate healthcare medical devices at home. Through a usability study based on the Grounded Theory methodology, we create a conceptual model in which we identified problems and errors related to the use of a multi-parameter monitor in a home environment by a group of lay users. The conceptual model is reified as a graphical representation, which allows stakeholders to identify (i) the weaknesses of the device, (ii) unsafe operation modes, and (iii) the most suitable device for a specific user.

Coactivación Muscular. Una Propuesta para Calcular su Índice utilizando Python / Muscle coactivation. A proposal to calculate its index utilising Python

Introducción. Fenómenos neurofisiológicos, como la coactivación muscular, se han utilizado para identificar tareas motoras que requieren una mayor estabilidad articular en personas sanas o con trastornos del movimiento. Sin embargo, existen varias formas de calcular el índice de coactivación (IC) muscular. Objetivo. El objetivo de este artículo fue crear una propuesta de procesamiento para calcular el IC muscular mediante el diseño de dos funciones utilizando el lenguaje Python. La primera función calcula el IC utilizando la fórmula planteada por Falconer y Winter, definida como "coactivation index". Se requiere introducir dos señales de músculos antagonistas con una misma longitud de datos y frecuencia de muestreo. Estas señales son previamente normalizadas a la contracción voluntaria máxima utilizando valores promedios rectificados. La segunda función definida como "plot_coactivacion" despliega una figura con los cambios de amplitud de ambos músculos y su área común. Estas funciones fueron creadas con un lenguaje de libre acceso (Python), destacando su clara sintaxis y la amplia gama de librerías en el procesamiento de señales biomédicas.

Introduction. Neurophysiological phenomena, such as muscle coactivation, have been used to identify motor tasks requiring greater joint stability in healthy people or with movement disorders. Nonetheless, there are many ways to calculate the coactivation index (CI). This article aimed to create a processing pipeline to calculate the muscular CI by designing two functions with the Python language. The first function calculates the CI utilising the formula proposed by Falconer and Winter, defined as "coactivation_index". It is required to introduce two signals of antagonist muscles with the same data long and sample frequency. These signals were previously normalised to the maximum voluntary contraction using the averaged rectified values. The second function was defined as "plot_coactivacion", which unfolds a figure that describes the amplitude changes for both muscles and their common area. These functions were designed with a freely accessible language (Python), highlighting its clear syntax and the number of libraries associated with biomedical signal processing.