The Brain/MINDS Marmoset Connectivity Resource: An open-access platform for cellular-level tracing and tractography in the primate brain.

The primate brain has unique anatomical characteristics, which translate into advanced cognitive, sensory, and motor abilities. Thus, it is important that we gain insight on its structure to provide a solid basis for models that will clarify function. Here, we report on the implementation and features of the Brain/MINDS Marmoset Connectivity Resource (BMCR), a new open-access platform that provides access to high-resolution anterograde neuronal tracer data in the marmoset brain, integrated to retrograde tracer and tractography data. Unlike other existing image explorers, the BMCR allows visualization of data from different individuals and modalities in a common reference space. This feature, allied to an unprecedented high resolution, enables analyses of features such as reciprocity, directionality, and spatial segregation of connections. The present release of the BMCR focuses on the prefrontal cortex (PFC), a uniquely developed region of the primate brain that is linked to advanced cognition, including the results of 52 anterograde and 164 retrograde tracer injections in the cortex of the marmoset. Moreover, the inclusion of tractography data from diffusion MRI allows systematic analyses of this noninvasive modality against gold-standard cellular connectivity data, enabling detection of false positives and negatives, which provide a basis for future development of tractography. This paper introduces the BMCR image preprocessing pipeline and resources, which include new tools for exploring and reviewing the data.

Dopamine D2 receptors in discrimination learning and spine enlargement.

Dopamine D2 receptors (D2Rs) are densely expressed in the striatum and have been linked to neuropsychiatric disorders such as schizophrenia1,2. High-affinity binding of dopamine suggests that D2Rs detect transient reductions in dopamine concentration (the dopamine dip) during punishment learning3-5. However, the nature and cellular basis of D2R-dependent behaviour are unclear. Here we show that tone reward conditioning induces marked stimulus generalization in a manner that depends on dopamine D1 receptors (D1Rs) in the nucleus accumbens (NAc) of mice, and that discrimination learning refines the conditioning using a dopamine dip. In NAc slices, a narrow dopamine dip (as short as 0.4 s) was detected by D2Rs to disinhibit adenosine A2A receptor (A2AR)-mediated enlargement of dendritic spines in D2R-expressing spiny projection neurons (D2-SPNs). Plasticity-related signalling by Ca2+/calmodulin-dependent protein kinase II and A2ARs in the NAc was required for discrimination learning. By contrast, extinction learning did not involve dopamine dips or D2-SPNs. Treatment with methamphetamine, which dysregulates dopamine signalling, impaired discrimination learning and spine enlargement, and these impairments were reversed by a D2R antagonist. Our data show that D2Rs refine the generalized reward learning mediated by D1Rs.

Cellular-resolution gene expression profiling in the neonatal marmoset brain reveals dynamic species- and region-specific differences.

Precise spatiotemporal control of gene expression in the developing brain is critical for neural circuit formation, and comprehensive expression mapping in the developing primate brain is crucial to understand brain function in health and disease. Here, we developed an unbiased, automated, large-scale, cellular-resolution in situ hybridization (ISH)-based gene expression profiling system (GePS) and companion analysis to reveal gene expression patterns in the neonatal New World marmoset cortex, thalamus, and striatum that are distinct from those in mice. Gene-ontology analysis of marmoset-specific genes revealed associations with catalytic activity in the visual cortex and neuropsychiatric disorders in the thalamus. Cortically expressed genes with clear area boundaries were used in a three-dimensional cortical surface mapping algorithm to delineate higher-order cortical areas not evident in two-dimensional ISH data. GePS provides a powerful platform to elucidate the molecular mechanisms underlying primate neurobiology and developmental psychiatric and neurological disorders.

Impact of CT scan parameters on deformable image registration accuracy using deformable thorax phantom.

The purpose of this study was to evaluate the deformable image registration (DIR) accuracy using various CT scan parameters with deformable thorax phantom. Our developed deformable thorax phantom (Dephan, Chiyoda Technol Corp, Tokyo, Japan) was used. The phantom consists of a base phantom, an inner phantom, and a motor-derived piston. The base phantom is an acrylic cylinder phantom with a diameter of 180 mm, which simulates the chest wall. The inner phantom consists of deformable, 20 mm thick disk-shaped sponges with 48 Lucite beads and 48 nylon cross-wires which simulate the vascular and bronchial bifurcations of the lung. Peak-exhale and peak-inhale images of the deformable phantom were acquired using a CT scanner (Aquilion LB, TOSHIBA). To evaluate the impact of CT scan parameters on DIR accuracy, we used the four tube voltages (80, 100, 120, and 135 kV) and six reconstruction algorithms (FC11, FC13, FC15, FC41, FC44, and FC52). Intensity-based DIR was performed between the two images using MIM Maestro (MIM software, Cleveland, USA). Fiducial markers (beads and cross-wires) based target registration error (TRE) was used for quantitative evaluation of DIR. In case with different tube voltages, the range of average TRE were 4.44-5.69 mm (reconstruction algorithm: FC13). In case with different reconstruction algorithms, the range of average TRE were 4.26-4.59 mm (tube voltage: 120 kV). The TRE were differed by up to 3.0 mm (3.96-6.96 mm) depending on the combination of tube voltage and reconstruction algorithm. Our result indicated that CT scan parameters had moderate impact of TRE, especially for reconstruction algorithms for the deformable thorax phantom.

The critical balance between dopamine D2 receptor and RGS for the sensitive detection of a transient decay in dopamine signal.

In behavioral learning, reward-related events are encoded into phasic dopamine (DA) signals in the brain. In particular, unexpected reward omission leads to a phasic decrease in DA (DA dip) in the striatum, which triggers long-term potentiation (LTP) in DA D2 receptor (D2R)-expressing spiny-projection neurons (D2 SPNs). While this LTP is required for reward discrimination, it is unclear how such a short DA-dip signal (0.5-2 s) is transferred through intracellular signaling to the coincidence detector, adenylate cyclase (AC). In the present study, we built a computational model of D2 signaling to determine conditions for the DA-dip detection. The DA dip can be detected only if the basal DA signal sufficiently inhibits AC, and the DA-dip signal sufficiently disinhibits AC. We found that those two requirements were simultaneously satisfied only if two key molecules, D2R and regulators of G protein signaling (RGS) were balanced within a certain range; this balance has indeed been observed in experimental studies. We also found that high level of RGS was required for the detection of a 0.5-s short DA dip, and the analytical solutions for these requirements confirmed their universality. The imbalance between D2R and RGS is associated with schizophrenia and DYT1 dystonia, both of which are accompanied by abnormal striatal LTP. Our simulations suggest that D2 SPNs in patients with schizophrenia and DYT1 dystonia cannot detect short DA dips. We finally discussed that such psychiatric and movement disorders can be understood in terms of the imbalance between D2R and RGS.

Signaling models for dopamine-dependent temporal contiguity in striatal synaptic plasticity.

Animals remember temporal links between their actions and subsequent rewards. We previously discovered a synaptic mechanism underlying such reward learning in D1 receptor (D1R)-expressing spiny projection neurons (D1 SPN) of the striatum. Dopamine (DA) bursts promote dendritic spine enlargement in a time window of only a few seconds after paired pre- and post-synaptic spiking (pre-post pairing), which is termed as reinforcement plasticity (RP). The previous study has also identified underlying signaling pathways; however, it still remains unclear how the signaling dynamics results in RP. In the present study, we first developed a computational model of signaling dynamics of D1 SPNs. The D1 RP model successfully reproduced experimentally observed protein kinase A (PKA) activity, including its critical time window. In this model, adenylate cyclase type 1 (AC1) in the spines/thin dendrites played a pivotal role as a coincidence detector against pre-post pairing and DA burst. In particular, pre-post pairing (Ca2+ signal) stimulated AC1 with a delay, and the Ca2+-stimulated AC1 was activated by the DA burst for the asymmetric time window. Moreover, the smallness of the spines/thin dendrites is crucial to the short time window for the PKA activity. We then developed a RP model for D2 SPNs, which also predicted the critical time window for RP that depended on the timing of pre-post pairing and phasic DA dip. AC1 worked for the coincidence detector in the D2 RP model as well. We further simulated the signaling pathway leading to Ca2+/calmodulin-dependent protein kinase II (CaMKII) activation and clarified the role of the downstream molecules of AC1 as the integrators that turn transient input signals into persistent spine enlargement. Finally, we discuss how such timing windows guide animals' reward learning.

Noise-resistant developmental reproducibility in vertebrate somite formation.

The reproducibility of embryonic development is remarkable, although molecular processes are intrinsically stochastic at the single-cell level. How the multicellular system resists the inevitable noise to acquire developmental reproducibility constitutes a fundamental question in developmental biology. Toward this end, we focused on vertebrate somitogenesis as a representative system, because somites are repeatedly reproduced within a single embryo whereas such reproducibility is lost in segmentation clock gene-deficient embryos. However, the effect of noise on developmental reproducibility has not been fully investigated, because of the technical difficulty in manipulating the noise intensity in experiments. In this study, we developed a computational model of ERK-mediated somitogenesis, in which bistable ERK activity is regulated by an FGF gradient, cell-cell communication, and the segmentation clock, subject to the intrinsic noise. The model simulation generated our previous in vivo observation that the ERK activity was distributed in a step-like gradient in the presomitic mesoderm, and its boundary was posteriorly shifted by the clock in a stepwise manner, leading to regular somite formation. Here, we showed that this somite regularity was robustly maintained against the noise. Removing the clock from the model predicted that the stepwise shift of the ERK activity occurs at irregular timing with irregular distance owing to the noise, resulting in somite size variation. This model prediction was recently confirmed by live imaging of ERK activity in zebrafish embryos. Through theoretical analysis, we presented a mechanism by which the clock reduces the inherent somite irregularity observed in clock-deficient embryos. Therefore, this study indicates a novel role of the segmentation clock in noise-resistant developmental reproducibility.

System level analysis of motor-related neural activities in larval Drosophila.

The way in which the central nervous system (CNS) governs animal movement is complex and difficult to solve solely by the analyses of muscle movement patterns. We tackle this problem by observing the activity of a large population of neurons in the CNS of larval Drosophila. We focused on two major behaviors of the larvae - forward and backward locomotion - and analyzed the neuronal activity related to these behaviors during the fictive locomotion that occurs spontaneously in the isolated CNS. We expressed a genetically-encoded calcium indicator, GCaMP and a nuclear marker in all neurons and then used digitally scanned light-sheet microscopy to record (at a fast frame rate) neural activities in the entire ventral nerve cord (VNC). We developed image processing tools that automatically detected the cell position based on the nuclear staining and allocate the activity signals to each detected cell. We also applied a machine learning-based method that we recently developed to assign motor status in each time frame. Our experimental procedures and computational pipeline enabled systematic identification of neurons that showed characteristic motor activities in larval Drosophila. We found cells whose activity was biased toward forward locomotion and others biased toward backward locomotion. In particular, we identified neurons near the boundary of the subesophageal zone (SEZ) and thoracic neuromeres, which were strongly active during an early phase of backward but not forward fictive locomotion.

Inverse tissue mechanics of cell monolayer expansion.

Living tissues undergo deformation during morphogenesis. In this process, cells generate mechanical forces that drive the coordinated cell motion and shape changes. Recent advances in experimental and theoretical techniques have enabled in situ measurement of the mechanical forces, but the characterization of mechanical properties that determine how these forces quantitatively affect tissue deformation remains challenging, and this represents a major obstacle for the complete understanding of morphogenesis. Here, we proposed a non-invasive reverse-engineering approach for the estimation of the mechanical properties, by combining tissue mechanics modeling and statistical machine learning. Our strategy is to model the tissue as a continuum mechanical system and to use passive observations of spontaneous tissue deformation and force fields to statistically estimate the model parameters. This method was applied to the analysis of the collective migration of Madin-Darby canine kidney cells, and the tissue flow and force were simultaneously observed by the phase contrast imaging and traction force microscopy. We found that our monolayer elastic model, whose elastic moduli were reverse-engineered, enabled a long-term forecast of the traction force fields when given the tissue flow fields, indicating that the elasticity contributes to the evolution of the tissue stress. Furthermore, we investigated the tissues in which myosin was inhibited by blebbistatin treatment, and observed a several-fold reduction in the elastic moduli. The obtained results validate our framework, which paves the way to the estimation of mechanical properties of living tissues during morphogenesis.

Identification of animal behavioral strategies by inverse reinforcement learning.

Animals are able to reach a desired state in an environment by controlling various behavioral patterns. Identification of the behavioral strategy used for this control is important for understanding animals' decision-making and is fundamental to dissect information processing done by the nervous system. However, methods for quantifying such behavioral strategies have not been fully established. In this study, we developed an inverse reinforcement-learning (IRL) framework to identify an animal's behavioral strategy from behavioral time-series data. We applied this framework to C. elegans thermotactic behavior; after cultivation at a constant temperature with or without food, fed worms prefer, while starved worms avoid the cultivation temperature on a thermal gradient. Our IRL approach revealed that the fed worms used both the absolute temperature and its temporal derivative and that their behavior involved two strategies: directed migration (DM) and isothermal migration (IM). With DM, worms efficiently reached specific temperatures, which explains their thermotactic behavior when fed. With IM, worms moved along a constant temperature, which reflects isothermal tracking, well-observed in previous studies. In contrast to fed animals, starved worms escaped the cultivation temperature using only the absolute, but not the temporal derivative of temperature. We also investigated the neural basis underlying these strategies, by applying our method to thermosensory neuron-deficient worms. Thus, our IRL-based approach is useful in identifying animal strategies from behavioral time-series data and could be applied to a wide range of behavioral studies, including decision-making, in other organisms.

Reconstruction of Spatial Thermal Gradient Encoded in Thermosensory Neuron AFD in Caenorhabditis elegans.

During navigation, animals process temporal sequences of sensory inputs to evaluate the surrounding environment. Thermotaxis of Caenorhabditis elegans is a favorable sensory behavior to elucidate how navigating animals process sensory signals from the environment. Sensation and storage of temperature information by a bilaterally symmetric pair of thermosensory neurons, AFD, is essential for the animals to migrate toward the memorized temperature on a thermal gradient. However, the encoding mechanisms of the spatial environment with the temporal AFD activity during navigation remain to be elucidated. Here, we show how the AFD neuron encodes sequences of sensory inputs to perceive spatial thermal environment. We used simultaneous calcium imaging and tracking system for a freely moving animal and characterized the response property of AFD to the thermal stimulus during thermotaxis. We show that AFD neurons respond to shallow temperature increases with intermittent calcium pulses and detect temperature differences with a critical time window of 20 s, which is similar to the timescale of behavioral elements of C. elegans, such as turning. Convolution of a thermal stimulus and the identified response property successfully reconstructs AFD activity. Conversely, deconvolution of the identified response kernel and AFD activity reconstructs the shallow thermal gradient with migration trajectory, indicating that AFD activity and the migration trajectory are sufficient as the encoded signals for thermal environment. Our study demonstrates bidirectional transformation between environmental thermal information and encoded neural activity. SIGNIFICANCE STATEMENT: Deciphering how information is encoded in the nervous system is an important challenge for understanding the principles of information processing in neural circuits. During navigation behavior, animals transform spatial information to temporal patterns of neural activity. To elucidate how a sensory system achieves this transformation, we focused on a thermosensory neuron in Caenorhabditis elegans called AFD, which plays a major role in a sensory behavior. Using tracking and calcium imaging system for freely moving animals, we identified the response property of the AFD. The identified response property enabled us to reconstruct both neural activity from a temperature stimulus and a spatial thermal environment from neural activity. These results shed light on how a sensory system encodes the environment.

Uncertainty-Dependent Extinction of Fear Memory in an Amygdala-mPFC Neural Circuit Model.

Uncertainty of fear conditioning is crucial for the acquisition and extinction of fear memory. Fear memory acquired through partial pairings of a conditioned stimulus (CS) and an unconditioned stimulus (US) is more resistant to extinction than that acquired through full pairings; this effect is known as the partial reinforcement extinction effect (PREE). Although the PREE has been explained by psychological theories, the neural mechanisms underlying the PREE remain largely unclear. Here, we developed a neural circuit model based on three distinct types of neurons (fear, persistent and extinction neurons) in the amygdala and medial prefrontal cortex (mPFC). In the model, the fear, persistent and extinction neurons encode predictions of net severity, of unconditioned stimulus (US) intensity, and of net safety, respectively. Our simulation successfully reproduces the PREE. We revealed that unpredictability of the US during extinction was represented by the combined responses of the three types of neurons, which are critical for the PREE. In addition, we extended the model to include amygdala subregions and the mPFC to address a recent finding that the ventral mPFC (vmPFC) is required for consolidating extinction memory but not for memory retrieval. Furthermore, model simulations led us to propose a novel procedure to enhance extinction learning through re-conditioning with a stronger US; strengthened fear memory up-regulates the extinction neuron, which, in turn, further inhibits the fear neuron during re-extinction. Thus, our models increased the understanding of the functional roles of the amygdala and vmPFC in the processing of uncertainty in fear conditioning and extinction.

Empirical Bayesian significance measure of neuronal spike response.

BACKGROUND: Functional connectivity analyses of multiple neurons provide a powerful bottom-up approach to reveal functions of local neuronal circuits by using simultaneous recording of neuronal activity. A statistical methodology, generalized linear modeling (GLM) of the spike response function, is one of the most promising methodologies to reduce false link discoveries arising from pseudo-correlation based on common inputs. Although recent advancement of fluorescent imaging techniques has increased the number of simultaneously recoded neurons up to the hundreds or thousands, the amount of information per pair of neurons has not correspondingly increased, partly because of the instruments' limitations, and partly because the number of neuron pairs increase in a quadratic manner. Consequently, the estimation of GLM suffers from large statistical uncertainty caused by the shortage in effective information. RESULTS: In this study, we propose a new combination of GLM and empirical Bayesian testing for the estimation of spike response functions that enables both conservative false discovery control and powerful functional connectivity detection. We compared our proposed method's performance with those of sparse estimation of GLM and classical Granger causality testing. Our method achieved high detection performance of functional connectivity with conservative estimation of false discovery rate and q values in case of information shortage due to short observation time. We also showed that empirical Bayesian testing on arbitrary statistics in place of likelihood-ratio statistics reduce the computational cost without decreasing the detection performance. When our proposed method was applied to a functional multi-neuron calcium imaging dataset from the rat hippocampal region, we found significant functional connections that are possibly mediated by AMPA and NMDA receptors. CONCLUSIONS: The proposed empirical Bayesian testing framework with GLM is promising especially when the amount of information per a neuron pair is small because of growing size of observed network.

Learning a common dictionary for subject-transfer decoding with resting calibration.

Brain signals measured over a series of experiments have inherent variability because of different physical and mental conditions among multiple subjects and sessions. Such variability complicates the analysis of data from multiple subjects and sessions in a consistent way, and degrades the performance of subject-transfer decoding in a brain-machine interface (BMI). To accommodate the variability in brain signals, we propose 1) a method for extracting spatial bases (or a dictionary) shared by multiple subjects, by employing a signal-processing technique of dictionary learning modified to compensate for variations between subjects and sessions, and 2) an approach to subject-transfer decoding that uses the resting-state activity of a previously unseen target subject as calibration data for compensating for variations, eliminating the need for a standard calibration based on task sessions. Applying our methodology to a dataset of electroencephalography (EEG) recordings during a selective visual-spatial attention task from multiple subjects and sessions, where the variability compensation was essential for reducing the redundancy of the dictionary, we found that the extracted common brain activities were reasonable in the light of neuroscience knowledge. The applicability to subject-transfer decoding was confirmed by improved performance over existing decoding methods. These results suggest that analyzing multisubject brain activities on common bases by the proposed method enables information sharing across subjects with low-burden resting calibration, and is effective for practical use of BMI in variable environments.

A statistical method of identifying interactions in neuron-glia systems based on functional multicell Ca2+ imaging.

Crosstalk between neurons and glia may constitute a significant part of information processing in the brain. We present a novel method of statistically identifying interactions in a neuron-glia network. We attempted to identify neuron-glia interactions from neuronal and glial activities via maximum-a-posteriori (MAP)-based parameter estimation by developing a generalized linear model (GLM) of a neuron-glia network. The interactions in our interest included functional connectivity and response functions. We evaluated the cross-validated likelihood of GLMs that resulted from the addition or removal of connections to confirm the existence of specific neuron-to-glia or glia-to-neuron connections. We only accepted addition or removal when the modification improved the cross-validated likelihood. We applied the method to a high-throughput, multicellular in vitro Ca2+ imaging dataset obtained from the CA3 region of a rat hippocampus, and then evaluated the reliability of connectivity estimates using a statistical test based on a surrogate method. Our findings based on the estimated connectivity were in good agreement with currently available physiological knowledge, suggesting our method can elucidate undiscovered functions of neuron-glia systems.

In vitro reconstitution of a CaMKII memory switch by an NMDA receptor-derived peptide.

Ca(2+)/Calmodulin-dependent protein kinase II (CaMKII) has been shown to play a major role in establishing memories through complex molecular interactions including phosphorylation of multiple synaptic targets. However, it is still controversial whether CaMKII itself serves as a molecular memory because of a lack of direct evidence. Here, we show that a single holoenzyme of CaMKII per se serves as an erasable molecular memory switch. We reconstituted Ca(2+)/Calmodulin-dependent CaMKII autophosphorylation in the presence of protein phosphatase 1 in vitro, and found that CaMKII phosphorylation shows a switch-like response with history dependence (hysteresis) only in the presence of an N-methyl-D-aspartate receptor-derived peptide. This hysteresis is Ca(2+) and protein phosphatase 1 concentration-dependent, indicating that the CaMKII memory switch is not simply caused by an N-methyl-D-aspartate receptor-derived peptide lock of CaMKII in an active conformation. Mutation of a phosphorylation site of the peptide shifted the Ca(2+) range of hysteresis. These functions may be crucial for induction and maintenance of long-term synaptic plasticity at hippocampal synapses.

Decoding spatial attention by using cortical currents estimated from electroencephalography with near-infrared spectroscopy prior information.

For practical brain-machine interfaces (BMIs), electroencephalography (EEG) and near-infrared spectroscopy (NIRS) are the only current methods that are non-invasive and available in non-laboratory environments. However, the use of EEG and NIRS involves certain inherent problems. EEG signals are generally a mixture of neural activity from broad areas, some of which may not be related to the task targeted by BMI, hence impairing BMI performance. NIRS has an inherent time delay as it measures blood flow, which therefore detracts from practical real-time BMI utility. To try to improve real environment EEG-NIRS-based BMIs, we propose here a novel methodology in which the subjects' mental states are decoded from cortical currents estimated from EEG, with the help of information from NIRS. Using a Variational Bayesian Multimodal EncephaloGraphy (VBMEG) methodology, we incorporated a novel form of NIRS-based prior to capture event related desynchronization from isolated current sources on the cortical surface. Then, we applied a Bayesian logistic regression technique to decode subjects' mental states from further sparsified current sources. Applying our methodology to a spatial attention task, we found our EEG-NIRS-based decoder exhibited significant performance improvement over decoding methods based on EEG sensor signals alone. The advancement of our methodology, decoding from current sources sparsely isolated on the cortex, was also supported by neuroscientific considerations; intraparietal sulcus, a region known to be involved in spatial attention, was a key responsible region in our task. These results suggest that our methodology is not only a practical option for EEG-NIRS-based BMI applications, but also a potential tool to investigate brain activity in non-laboratory and naturalistic environments.

Belief inference for hierarchical hidden states in spatial navigation.

Uncertainty abounds in the real world, and in environments with multiple layers of unobservable hidden states, decision-making requires resolving uncertainties based on mutual inference. Focusing on a spatial navigation problem, we develop a Tiger maze task that involved simultaneously inferring the local hidden state and the global hidden state from probabilistically uncertain observation. We adopt a Bayesian computational approach by proposing a hierarchical inference model. Applying this to human task behaviour, alongside functional magnetic resonance brain imaging, allows us to separate the neural correlates associated with reinforcement and reassessment of belief in hidden states. The imaging results also suggest that different layers of uncertainty differentially involve the basal ganglia and dorsomedial prefrontal cortex, and that the regions responsible are organised along the rostral axis of these areas according to the type of inference and the level of abstraction of the hidden state, i.e. higher-order state inference involves more anterior parts.

Shredding artifacts: extracting brain activity in EEG from extreme artifacts during skateboarding using ASR and ICA.

Introduction: To understand brain function in natural real-world settings, it is crucial to acquire brain activity data in noisy environments with diverse artifacts. Electroencephalography (EEG), while susceptible to environmental and physiological artifacts, can be cleaned using advanced signal processing techniques like Artifact Subspace Reconstruction (ASR) and Independent Component Analysis (ICA). This study aims to demonstrate that ASR and ICA can effectively extract brain activity from the substantial artifacts occurring while skateboarding on a half-pipe ramp. Methods: A dual-task paradigm was used, where subjects were presented with auditory stimuli during skateboarding and rest conditions. The effectiveness of ASR and ICA in cleaning artifacts was evaluated using a support vector machine to classify the presence or absence of a sound stimulus in single-trial EEG data. The study evaluated the effectiveness of ASR and ICA in artifact cleaning using five different pipelines: (1) Minimal cleaning (bandpass filtering), (2) ASR only, (3) ICA only, (4) ICA followed by ASR (ICAASR), and (5) ASR preceding ICA (ASRICA). Three skateboarders participated in the experiment. Results: Results showed that all ICA-containing pipelines, especially ASRICA (69%, 68%, 63%), outperformed minimal cleaning (55%, 52%, 50%) in single-trial classification during skateboarding. The ASRICA pipeline performed significantly better than other pipelines containing ICA for two of the three subjects, with no other pipeline performing better than ASRICA. The superior performance of ASRICA likely results from ASR removing non-stationary artifacts, enhancing ICA decomposition. Evidenced by ASRICA identifying more brain components via ICLabel than ICA alone or ICAASR for all subjects. For the rest condition, with fewer artifacts, the ASRICA pipeline (71%, 82%, 75%) showed slight improvement over minimal cleaning (73%, 70%, 72%), performing significantly better for two subjects. Discussion: This study demonstrates that ASRICA can effectively clean artifacts to extract single-trial brain activity during skateboarding. These findings affirm the feasibility of recording brain activity during physically demanding tasks involving substantial body movement, laying the groundwork for future research into the neural processes governing complex and coordinated body movements.

Expression profiling using a tumor-specific cDNA microarray predicts the prognosis of intermediate risk neuroblastomas.

To predict the prognosis of neuroblastoma patients and choose a better therapeutic protocol, we developed a cDNA microarray carrying 5340 genes obtained from primary neuroblastomas and examined 136 tumor samples. We made a probabilistic output statistical classifier that provided a high accuracy in prognosis prediction (89% at 5 years) and a highly reliable method to validate it. Kaplan-Meier analysis indicated that the patients in an intermediate group defined by existing markers are divided by microarray into two further groups with 5 year survivals for 36% and 89% of patients (p < 10(-4)), i.e., with unfavorably and favorably predicted neuroblastomas, respectively. According to these results, we developed a gene subset chip for a clinical tool, for which our classifier exhibited 88% prediction accuracy.