Neural specificity of acupuncture stimulation from support vector machine classification analysis

Acupoint specificity, as a crucial issue in acupuncture neuroimaging studies, is still a controversial topic. Previous studies have generally adopted a block-based general linear model (GLM) approach, which predicts the temporal changes in the blood oxygenation level-dependent signal conforming to the “on-off” specifications. However, this method might become impractical since the precise timing and duration of acupuncture actions cannot be specified a priori. In the current study, we applied a data-driven multivariate classification approach, namely, support vector machine (SVM), to explore the neural specificity of acupuncture at gall bladder 40 (GB40) using kidney 3 (KI3) as a control condition (belonging to different meridians but the same nerve segment). In addition, to verify whether the typical GLM approach is sensitive enough in exploring the neural response patterns evoked by acupuncture, we also employed the GLM method to the same data sets. The SVM analysis detected distinct neural response patterns between GB40 and KI3 — positive predominantly for the GB40, while negative following the KI3. By contrast, group analysis from the GLM showed that acupuncture at these different acupoints can both evoke similar widespread signal decreases in multiple brain regions, and most of these regions were spatially overlapped, mainly distributing in the limbic and subcortical structures. Our findings may provide additional evidence to support the specificity of acupuncture, relevant to its clinical efficacy. Moreover, we also proved that GLM analysis is prone to be susceptible to errors and is not appropriate for detecting neural response patterns evoked by acupuncture stimulation.     

Exploring vision-related acupuncture point specificity with multivoxel pattern analysis

Acupoint specificity is one of the central issues of functional magnetic resonance imaging (fMRI) studies of acupuncture and has been under discussed. However, strong and consistent proof has not been provided for the existence of acupoint specificity, and unsuitable analysis approach applied could be the reason. We observed that previous researches of acupoint specificity were mostly based on model-based methods which were limited to make exploration of acupoint specificity because of the inaccurate specified prior. Here we applied multi-voxel pattern analysis (MVPA) to investigate the specificity of brain activation patterns induced by acupuncture stimulations at a vision-related acupoint (GB37) and a nearby nonacupoint (NAP). Results showed that multiple brain areas could differentiate the central neural response patterns induced by acupuncture stimulation at these two sites with higher accuracy above the chance level. These regions included occipital cortex, limbic-cerebellar areas and somatosensory cortex. Our results support that the characteristic neural response patterns of brain cortex to the acupuncture stimulation at GB37 and a nearby NAP could differ from each other effectively with the application of MVPA approach.     

FMRI connectivity analysis of acupuncture effects on the whole brain network in mild cognitive impairment patients

The increased risk for the elderly with mild cognitive impairment (MCI) to progress to Alzheimer's disease makes it an appropriate condition for investigation. While the use of acupuncture as a complementary therapeutic method for treating MCI is popular in certain parts of the world, the underlying mechanism is still elusive. We sought to investigate the acupuncture effects on the functional connectivity throughout the entire brain in MCI patients compared to healthy controls (HC). The functional magnetic resonance imaging experiment was performed with two different paradigms, namely, deep acupuncture (DA) and superficial acupuncture (SA), at acupoint KI3. We first identified regions showing abnormal functional connectivity in the MCI group compared to HC during the resting state and subsequently tested whether these regions could be modulated by acupuncture. Then, we made the comparison of MCI vs. HC to test whether there were any specific modulatory patterns in the poststimulus resting brain between the two groups. Finally, we made the comparisons of DA vs. SA in each group to test the effect of acupuncture with different needling depths. We found the temporal regions (hippocampus, thalamus, fusiform gyrus) showing abnormal functional connectivity during the resting state. These regions are implicated in memory encoding and retrieving. Furthermore, we found significant changes in functional connectivity related with the abnormal regions in MCI patients following acupuncture. Compared to HC, the correlations related with the temporal regions were enhanced in the poststimulus resting brain in MCI patients. Compared to SA, significantly increased correlations related with the temporal regions were found for the DA condition. The enhanced correlations in the memory-related brain regions following acupuncture may be related to the purported therapeutically beneficial effects of acupuncture for the treatment of MCI. The heterogeneous modulatory patterns between DA and SA may suggest that deep muscle insertion of acupuncture is necessary to achieve the appreciable clinical effect.     

A serial functional connectivity MRI study in healthy individuals assessing the variability of connectivity measures: reduced interhemispheric connectivity in the motor network during continuous performance

To date, little data is available on the reproducibility of functional connectivity MRI (fcMRI) studies. Here, we tested the variability and reproducibility of both the functional connectivity itself and different statistical methods to analyze this phenomenon. In the main part of our study, we repeatedly examined two healthy subjects in 10 sessions over 6 months with fcMRI. Cortical areas involved in motor function were examined under two different cognitive states: during continuous performance (CP) of a flexion/extension task of the fingers of the right hand and while subjects were at rest. Connectivity to left primary motor cortex (lSM1) was calculated by correlation analysis. The resulting correlation coefficients were transformed to z-scores of the standard normal distribution. For each subject, multisession statistical analyses were carried out with the z-score maps of the resting state (RS) and the CP experiments. First, voxel based t tests between the two groups of fcMRI experiments were performed. Second, ROI analyses were carried out for contralateral right SM1 and for supplementary motor area (SMA). For both ROI, mean and maximum z-score were calculated for each experiment. Also, the fraction of significantly (P<.05) correlated voxels (FCV) in each ROI was calculated. To evaluate the differences between the RS and the CP condition, paired t tests were performed for the mean and maximum z-scores, and Wilcoxon signed ranks tests for matched pairs were carried out for the FCV. All statistical methods and connectivity measures under investigation yielded a distinct loss in left–right SM1 connectivity under the CP condition. For SMA, interindividual differences were apparent. We therefore repeated the fcMRI experiments and the ROI analyses in a group of seven healthy subjects (including the two subjects of the main study). In this substudy, we were able to verify the reduction of left–right SM1 connectivity during unilateral performance. Still, the direction of SMA to lSM1 connectivity change during the CP condition remained undefined as four subjects showed a connectivity increase and three showed a decrease. In summary, we were able to demonstrate a distinct reduction in left–right SM1 synchrony in the CP condition compared to the RS both in the longitudinal and in the multisubject study. This effect was reproducible with all statistical methods and all measures of connectivity under investigation. We conclude that despite intra- and interindividual variability, serial and cross-sectional assessment of functional connectivity reveals stable and reliable results.     

Post-stimulus response in hemodynamics observed by functional magnetic resonance imaging--difference between the primary sensorimotor area and the supplementary motor area

The blood oxygen level dependency (BOLD) contrast is a useful tool for functional neuroimaging based on the hemodynamic response to neuronal activation. We observed different hemodynamic responses in the BOLD signal between the primary sensorimotor area (SM1) and the supplementary motor area (SMA) in the sequential finger movement task. In the SMA, a stronger initial overshoot and a post-stimulus overshoot were observed. It was hypothesized from the time course analysis that the stronger initial overshoot reflected the activation of the SMA for motor control programming in the initial phase. Although the post-stimulus overshoot may be partially explained by cerebral blood flow (CBF) cerebral blood volume (CBV) uncoupling, its mechanism remained unknown. In the SM1, only the initial overshoot was observed and the level of BOLD signal was almost constant after the initial overshoot during the task period. These observations suggested that the BOLD signal is characterized by both CBF-CBV uncoupling and the neuronal activation characteristics in each region.     

高级脑活动的功能性核磁共振成像

功能性核磁共振成像技术可以显示大脑各个区域内静脉毛细血管中血液氧合状态所起的磁共振信号的小变化。使用ｆＭＲＩ的方法，可以在正常的活体上无损伤地实现大脑活动的功能定位，时空分辨率可分别达到秒和毫米数量级，尽管目前还面临一系列技术的困难，ｆＭＲＩ已经日益成为观察大脑活动，研究有脑的拓扑结构，进而揭示脑和思维关系的一种重要方法。     

Surface-based functional magnetic resonance imaging analysis of partial brain echo planar imaging data at 1.5 T

Surface-based functional magnetic resonance imaging (fMRI) analysis is more sensitive and accurate than volume-based analysis for detecting neural activation. However, these advantages are less important in practical fMRI experiments with commonly used 1.5-T magnetic resonance devices because of the resolution gap between the echo planar imaging data and the cortical surface models. We expected high-resolution segmented partial brain echo planar imaging (EPI) data to overcome this problem, and the activation patterns of the high-resolution data could be different from the low-resolution data. For the practical applications of surface-based fMRI analysis using segmented EPI techniques, the effects of some important factors (e.g., activation patterns, registration and local distortions) should be intensively evaluated because the results of surface-based fMRI analyses could be influenced by them. In this study, we demonstrated the difference between activations detected from low-resolution EPI data, which were covering whole brain, and high-resolution segmented EPI data covering partial brain by volume- and surface-based analysis methods. First, we compared the activation maps of low- and high-resolution EPI datasets detected by volume- and surface-based analyses, with the spatial patterns of activation clusters, and analyzed the distributions of activations in occipital lobes. We also analyzed the high-resolution EPI data covering motor areas and fusiform gyri of human brain, and presented the differences of activations detected by volume- and surface-based methods.     

Assessing functional connectivity in the human brain by fMRI

Functional magnetic resonance imaging (fMRI) is widely used to detect and delineate regions of the brain that change their level of activation in response to specific stimuli and tasks. Simple activation maps depict only the average level of engagement of different regions within distributed systems. FMRI potentially can reveal additional information about the degree to which components of large-scale neural systems are functionally coupled together to achieve specific tasks. In order to better understand how brain regions contribute to functionally connected circuits, it is necessary to record activation maps either as a function of different conditions, at different times or in different subjects. Data obtained under different conditions may then be analyzed by a variety of techniques to infer correlations and couplings between nodes in networks. Several multivariate statistical methods have been adapted and applied to analyze variations within such data. An approach of particular interest that is suited to studies of connectivity within single subjects makes use of acquisitions of runs of MRI images obtained while the brain is in a so-called steady state, either at rest (i.e., without any specific stimulus or task) or in a condition of continuous activation. Interregional correlations between fluctuations of MRI signal potentially reveal functional connectivity. Recent studies have established that interregional correlations between different components of circuits in each of the visual, language, motor and working memory systems can be detected in the resting state. Correlations at baseline are changed during the performance of a continuous task. In this review, various methods available for assessing connectivity are described and evaluated.     

Quantification of the power changes in BOLD signals using Welch spectrum method during different single-hand motor imageries

Motor imagery is an experimental paradigm implemented in cognitive neuroscience and cognitive psychology. To investigate the asymmetry of the strength of cortical functional activity due to different single-hand motor imageries, functional magnetic resonance imaging (fMRI) data from right handed normal subjects were recorded and analyzed during both left-hand and right-hand motor imagery processes. Then the average power of blood oxygenation level-dependent (BOLD) signals in temporal domain was calculated using the developed tool that combines Welch power spectrum and the integral of power spectrum approach of BOLD signal changes during motor imagery. Power change analysis results indicated that cortical activity exhibited a stronger power in the precentral gyrus and medial frontal gyrus with left-hand motor imagery tasks compared with that from right-hand motor imagery tasks. These observations suggest that right handed normal subjects mobilize more cortical nerve cells for left-hand motor imagery. Our findings also suggest that the approach based on power differences of BOLD signals is a suitable quantitative analysis tool for quantification of asymmetry of brain activity intensity during motor imagery tasks.     

Quantification of venous blood signal contribution to BOLD functional activation in the auditory cortex at 3 T

Most modern techniques for functional magnetic resonance imaging (fMRI) rely on blood-oxygen-level-dependent (BOLD) contrast as the basic principle for detecting neuronal activation. However, the measured BOLD effect depends on a transfer function related to neurophysiological changes accompanying electrical neural activation. The spatial accuracy and extension of the region of interest are determined by vascular effect, which introduces incertitude on real neuronal activation maps. Our efforts have been directed towards the development of a new methodology that is capable of combining morphological, vascular and functional information; obtaining new insight regarding foci of activation; and distinguishing the nature of activation on a pixel-by-pixel basis. Six healthy volunteers were studied in a parametric auditory functional experiment at 3 T; activation maps were overlaid on a high-resolution brain venography obtained through a novel technique. The BOLD signal intensities of vascular and nonvascular activated voxels were analyzed and compared: it was shown that nonvascular active voxels have lower values for signal peak (P<10(-7)) and area (P<10(-8)) with respect to vascular voxels. The analysis showed how venous blood influenced the measured BOLD signals, supplying a technique to filter possible venous artifacts that potentially can lead to misinterpretation of fMRI results. This methodology, although validated in the auditory cortex activation, maintains a general applicability to any cortical fMRI study, as the basic concepts on which it relies on are not limited to this cortical region. The results obtained in this study can represent the basis for new methodologies and tools that are capable of adding further characterization to the BOLD signal properties.     

Performance evaluation of nonnegative matrix factorization algorithms to estimate task-related neuronal activities from fMRI data

Nonnegative matrix factorization (NMF) is a blind source separation (BSS) algorithm which is based on the distinct constraint of nonnegativity of the estimated parameters as well as on the measured data. In this study, according to the potential feasibility of NMF for fMRI data, the four most popular NMF algorithms, corresponding to the following two types of (1) least-squares based update [i.e., alternating least-squares NMF (ALSNMF) and projected gradient descent NMF] and (2) multiplicative update (i.e., NMF based on Euclidean distance and NMF based on divergence cost function), were investigated by using them to estimate task-related neuronal activities. These algorithms were applied firstly to individual data from a single subject and, subsequently, to group data sets from multiple subjects. On the single-subject level, although all four algorithms detected task-related activation from simulated data, the performance of multiplicative update NMFs was significantly deteriorated when evaluated using visuomotor task fMRI data, for which they failed in estimating any task-related neuronal activities. In group-level analysis on both simulated data and real fMRI data, ALSNMF outperformed the other three algorithms. The presented findings may suggest that ALSNMF appears to be the most promising option among the tested NMF algorithms to extract task-related neuronal activities from fMRI data.     

BOLD functional MRI in disease and pharmacological studies: room for improvement?

In the past decade the use of blood oxygen level-dependent (BOLD) fMRI to investigate the effect of diseases and pharmacological agents on brain activity has increased greatly. BOLD fMRI does not measure neural activity directly, but relies on a cascade of physiological events linking neural activity to the generation of MRI signal. However, most of the disease and pharmacological studies performed so far have interpreted changes in BOLD fMRI as "brain activation," ignoring the potential confounds that can arise through drug- or disease-induced modulation of events downstream of the neural activity. This issue is especially serious in diseases (like multiple sclerosis, brain tumours and stroke) and drugs (like anaesthetics or those with a vascular action) that are known to influence these physiological events. Here we provide evidence that, to extract meaningful information on brain activity in patient and pharmacological BOLD fMRI studies, it is important to identify, characterise and possibly correct these influences that potentially confound the results. We suggest a series of experimental measures to improve the interpretability of BOLD fMRI studies. We have ranked these according to their potential information and current practical feasibility. First-line, necessary improvements consist of (1) the inclusion of one or more control tasks, and (2) the recording of physiological parameters during scanning and subsequent correction of possible between-group differences. Second-line, highly recommended important aim to make the results of a patient or drug BOLD study more interpretable and include the assessment of (1) baseline brain perfusion, (2) vascular reactivity, (3) the inclusion of stimulus-related perfusion fMRI and (4) the recording of electrophysiological responses to the stimulus of interest. Finally, third-line, desirable improvements consist of the inclusion of (1) simultaneous EEG-fMRI, (2) cerebral blood volume and (3) rate of metabolic oxygen consumption measurements and, when relevant, (4) animal studies investigating signalling between neural cells and blood vessels.     

Nonlinear blood oxygen level-dependent responses for transient activations and deactivations in V1 - insights into the hemodynamic response function with the balloon model

We report studies of the nonlinear nature of blood oxygen level-dependent (BOLD) responses to short transient deactivations in human visual cortex. Both functional magnetic resonance imaging (fMRI) and near-infrared spectroscopy (NIRS) have been used to compare and contrast the hemodynamic response functions (HRFs) associated with transient activation and deactivation in primary visual cortex. We show that signal decreases for short duration deactivations are smaller than corresponding signal increases in activation studies. Moreover, the standard balloon model of BOLD effects may be modified to account for the observed nonlinear nature of deactivations by appropriate changes to simple hemodynamic parameters without recourse to new assumptions about the nature of the coupling between activity and oxygen use.     

Independent component model of the default-mode brain function: combining individual-level and population-level analyses in resting-state fMRI

Resting-state functional magnetic resonance imaging (RS-fMRI) is a technique used to investigate the spontaneous correlations of blood-oxygen-level-dependent signals across different regions of the brain. Using functional connectivity tools, it is possible to investigate a specific RS-fMRI network, referred to as "default-mode" (DM) network, that involves cortical regions deactivated in fMRI experiments with cognitive tasks. Previous works have reported a significant effect of aging on DM regions activity. Independent component analysis (ICA) is often used for generating spatially distributed DM functional connectivity patterns from RS-fMRI data without the need for a reference region. This aspect and the relatively easy setup of an RS-fMRI experiment even in clinical trials have boosted the combined use of RS-fMRI and ICA-based DM analysis for noninvasive research of brain disorders. In this work, we considered different strategies for combining ICA results from individual-level and population-level analyses and used them to evaluate and predict the effect of aging on the DM component. Using RS-fMRI data from 20 normal subjects and a previously developed group-level ICA methodology, we generated group DM maps and showed that the overall ICA-DM connectivity is negatively correlated with age. A negative correlation of the ICA voxel weights with age existed in all DM regions at a variable degree. As an alternative approach, we generated a distributed DM spatial template and evaluated the correlation of each individual DM component fit to this template with age. Using a "leave-one-out" procedure, we discuss the importance of removing the bias from the DM template-generation process.     

Fractional amplitude of low-frequency fluctuation changes in functional dyspepsia: A resting-state fMRI study

Recently, there is an increasing interest in the study of the role of brain dysfunction in the pathogenesis of symptoms of functional dyspepsia (FD). More specifically, abnormal brain activities in patients with FD during the resting state have been proven by several positron emission tomography (PET) studies. Resting-state functional magnetic resonance imaging (fMRI) is also a valuable tool in investigating spontaneous brain activity abnormalities in pathological conditions. In the present study, we examined the amplitude of low-frequency fluctuations (ALFF) and fractional (f)ALFF changes in patients with FD by using fMRI. Twenty-nine patients with FD and sixteen healthy controls participated in this study. Between-group differences in ALFF/fALFF were examined using a permutation-based nonparametric test after accounting for the gender and age effects. The results revealed a significant between-group difference in fALFF but not in ALFF in multiple brain regions including the right insula, brainstem and cerebellum. Seed-based resting-state functional connectivity analysis revealed that FD patients have increased correlations between the right cerebellum and multiple brain regions including the bilateral brainstem, bilateral cerebellum, bilateral thalamus, left para-/hippocampus, left pallidum and left putamen. Furthermore, fLAFF values in the right insula were positively correlated with the severity of the disease. These findings have provided further evidence of spontaneous brain activity abnormalities in FD patients which might contribute to our understanding of the pathophysiology of the disease.     

Hadamard-encoded sub-slice fMRI for reduced signal dropout

T2?-weighted blood oxygen level-dependent functional magnetic resonance imaging is adversely affected by susceptibility-induced field gradients in brain regions adjoining air interfaces that cause image distortion and signal dropout. Reducing slice thickness diminishes signal dropout but is accompanied by reduced signal-to-noise ratio (SNR). This study proposes simultaneous excitation of subslices with total width equal to the desired slice thickness, employing alternating Hadamard-encoded radiofrequency pulses coupled with incoherent addition of the subslices to achieve reduction of through-plane dephasing with minimal SNR loss but at the expense of a reduction in temporal resolution. Using a sensory task and hypercapnic challenge with breathholding (BH), results with two subslices per slice and a twofold reduction in temporal resolution show improved activation relative to a conventional acquisition. Average (eight subjects) T-scores in the BH task increased by 16% (P<.0003), and activation extent increased by 12% (not significant). In frontal brain regions, significant improvements in BH activation extent (11.4%, P<.05) and T-scores (18%, P<.0002) were demonstrated. Higher temporal resolution can be achieved by tradeoff of SNR.     

Depression recognition using resting-state and event-related fMRI signals

Purpose

This paper aimed to develop a method for depression detection using blood-oxygen-level-dependent (BOLD) response estimated from event-related signals and resting-state functional magnetic resonance imaging (fMRI) signals together.

Materials and Methods

Thirteen patients with unipolar depression and matched healthy subjects were recruited. Resting state data of each subject were collected. Thereafter, event-related paradigm was undertaken using sad facial stimuli. The resting-state fMRI signal was deemed as the baseline of each subject's activity. Coefficient marks were designed to sort and select temporal independent components of event-related signals. Thereafter, stimulus-evoked BOLD response components inside event-related signal were extracted and taken as features to discriminate depressive patients from healthy controls.

Results

Accuracy rate for depression recognition was 77.27% with P value of .017 for whole-brain analysis and 81.82% with P value of .009 for region-of-interest analysis. The effectiveness and the superiority of the proposed method for disease recognition were demonstrated via the performance comparison with three other typical methods.

Conclusions

The proposed model was effective in depression recognition.     

MRI gradient fields increase brain mannitol space

Following nephrectomy and intravenous injection of tritiated mannitol, adult male rats were exposed to magnetic resonance imaging (MRI) procedures at 1.5 T, 0.5 T, and 0.3 T. Compared to rats similarly handled but not exposed to MRI procedures, brain mannitol concentration, expressed as a percentage of mean body concentration, was significantly increased at 0.3 T and 0.5 T but not at 1.5 T. At 0.3 T, exposure to gradient-field fluctuations used for imaging increased brain mannitol concentration, but exposures to static main field and pulsed radiofrequency energies did not. Increased brain mannitol associated with gradient-field flux may reflect increased blood-brain barrier permeability or blood volume in brain. MRI effects on brain mannitol space are of uncertain clinical significance, but are consistent with prior evidence of an MRI-induced increase of brain capillary endothelial cell transport observed with horseradish peroxidase. Further studies are needed to confirm these findings and to explore the processes underlying changes in mannitol distribution related to MRI.     

核磁共振成像教学实验

利用超小型核磁共振成像仪，完成了核磁共振空间编码成像的一系列实验．研究了超小型核磁共振成像仪的性能，并对2种观赏植物、动物骨骼、蛔虫、蚯蚓等生物样品进行了核磁共振成像以及T1加权和氢核密度加权成像．