3.
<ul><li>In our previous study, we found that the amplitude of the event-related potential (ERP) component N140 evoked by the tactile stimulus was increased when the subject expected a coming visual stimulus that had been paired with the tactile stimulus in comparison to this component evoked by the same tactile stimulus without crossmodal expectation . It has been suggested that the somatosensory N140 is generated by sources in multiple cortical areas, including frontal cortex and SII cortex. </li></ul>

4.
<ul><li>By applying independent component analysis (ICA) in the present study to the EEG (electroencephalogram) data recorded in the unimodal and crossmodal tasks, we explored independent components (ICs) that represented neural activities in cortical areas. We found that the crossmodal modulation of the N140 represented the neural activities in somatosensory (SI, and possibly SII as well) and frontal cortical areas that cooperated with each other in crossmodal association in the tasks. </li></ul>Linear noiseless ICA The components xi of the observed random vector X =(X1,….,Xm)T are generated as a sum of the independent components sk K=1,….,n Xi= ai,1s1 +…..+ ai,ksk The same generative model can be written in vectorial form as X= sigma(1,n) aksk

5.
<ul><li>Participants </li></ul><ul><li>Twelve paid normal adult volunteers were recruited for the present study (10 men, 2 women, aged 19–47 years). </li></ul><ul><li>Stimuli and EEG recording </li></ul><ul><li>Experiments were carried out in a quiet, dimly lit room. Participants sat in a comfortable chair, facing a light-emitting diode (LED) about 1.5 m away in the center of the visual field at eye level. Visual stimuli (LED) comprised green or red light 100 ms in duration. Tactile stimuli were generated by a mechanical vibrator with frequency 80 Hz and 150 Hz , and delivered on the subject's left index finger-tip. During performance of tasks, participants placed their left hand on a supporter in their usual position to receive the tactile stimulus , and their right hand on another supporter to press two buttons as the response with their fingers, the index finger for left button; the middle finger for right button. </li></ul><ul><li>Electroencephalograms (EEG) were recorded by an EEG recording system (SynAmp, Neuroscan, Ltd Corp). Thirty-two Ag-AgCl scalp electrodes (Quick-Caps, Neuroscan) were arranged in standard locations (Guideline thirteen, American Clinical Neurophysiology Society, 2003). EEG signals from all of these electrodes were referenced to linked earlobes. The impedance of each electrode was kept below 5 kΩ. The electro-oculogram (EOG) was recorded for horizontal eye movements and for vertical eye movements. Signals of EEG (30 electrodes) and EOG (2 electrodes) were filtered (0.1–100 Hz band-pass), amplified, digitized (500 Hz sampling rate), and saved for off-line analysis. </li></ul>

6.
Independent component analysis Analysis of EEG data recorded from 30 electrodes was performed by using Matlab 7.0 (Math Works, Natick, MA) and EEGLAB4.51(Swartz Center for Computational Neurosciences, La Jolla, CA; http://www.sccn.ucsd.edu/eeglab ), a freely available open source software toolbox. BESA2000 (MEGIS, Graefeling, Germany) was also used to localize dipoles of independent components (ICs). The filtered-EEG data (2~30 Hz) that preserved theta, alpha, and beta band information were used for the ICA study. The onset of S-1 in each trial was used as the task-event marker to separate a trial into a period before the onset and a period after the onset. In each trial, filtered-EEG of 1,500 ms (500 ms before the onset of S-1 and 1,000 ms after it) were extracted from the continuous EEG to form a data epoch. The mean value of EEG amplitude in the first 100 ms of the epoch was calculated from all trials of each task in each individual subject . To obtain the EEG data epoch for further processing, this mean value was then subtracted from each corresponding data epoch to reduce the influence of EEG variance across trials. All data epochs were put together and submitted to infomax ICA that comes from the ICA families performing blind source separation. A 30×30 unmixing square matrix was found by using Infomax ICA. When this matrix was multiplied by the EEG data epochs, maximally, temporally independent activities were obtained. In this calculation of the independent activities, a weight change of 10e-6 together with iterations <800 were set as the stop criterion .

7.
Independent component analysis Let X denote the EEG data and M denote the unmixing square matrix. Then independent activities (S) are: S = MX. We can change the formula into X = M−1S. In this formula, one row of the matrix S represents the temporal activity of one IC, and the corresponding column of the matrix M−1 represents this IC's spatial pattern at the scalp electrodes. The back-projection of an IC at one electrode is obtained by multiplying the temporal activity of this IC with its coefficient of the corresponding spatial pattern at this electrode. EEG at one electrode can be considered as the sum of back-projections of all ICs at this electrode. The temporal independent activity and its corresponding spatial pattern together characterize an IC that may correlate with the activity of a neuronal clique. In the present study, we screened activities of ICs to determine potentially common temporal patterns of those ICs across all trials of each task and subjects, and we also visually screened topographies of ICs to assess their potentially common spatial patterns across tasks and subjects. ICs showing event (onset of S-1)-related activities consistently across trials of each task and subjects, and spatial topographies consistently across tasks and subjects, were selected in the screening. The spatial topographies may reflect the dipole activity, presumably caused by partially synchronous activities within certain cortical source patches that produce far-field potentials through volume conduction. The above process of selection resulted in identification of ICs for each subject and each task, which presumably reflected activities of neuronal cliques in certain cortical areas.

8.
Grand average back-projections of the selected ICs to the scalp electrode sites were compared with the original scalp ERPs. The contribution of each IC to an original ERP component was assessed by calculating the proportion of the back-projection of the IC to that ERP component. Latency and amplitude of the components of both the original ERPs and the back-projections of ICs were statistically analyzed by performing a four-way repeated measures MANOVA with Modality, Type (ERPs, and IC's back-projected activity), AP, and LR as within-subject factors. The amplitude of the components of IC projections and ERPs was normalized by subtracting its corresponding baseline (200 ms preceding the onset of S-1) mean values before the statistical analysis. The grand average topography across subjects and tasks of a selected IC was submitted to BESA2000 that uses a standard four-shell spherical head model (i.e., brain, cerebrospinal fluid, bone, and scalp) to find the location of the IC-related dipole (source model) in the brain. The dipole was derived in BESA2000 by fitting it iteratively to the averaged IC topography parameters until minimal residual variance was reached. In the present study, the values of residual variance lower than 10% were used as the threshold [74] .

9.
Time-frequency representations (TFRs) and coherence TFRs of ICs and ERPs from the electrodes that showed the largest IC back-projections were computed on single trials in the frequency range of 2~30 Hz by using Hanning windowed short time Fourier transformation . The window had a fixed length of 250 ms, moving across every time stamp . The mean value of the windowed-period was taken away to avoid the variation of direct current. Zeros were then added after each windowed-period to make TFRs smoother across the Frequency axis. The ratio of the zero-pad to the windowed-period was 32. The TFRs were then normalized for each frequency by subtracting the baseline (200 ms before the onset of S-1) mean value, and dividing by the baseline standard deviation The coherence spectrum of the two independent components (ICs) was calculated by using: C 12 =[S 12 ] 2 /S 1 S 2 S 1 , S 2 are the power spectra of the two ICs respectively . S 12 is the cross spectrum between the two ICs . The value of coherence (C 12 ) at frequency f ranges from 1, indicating maximum interdependence between the two ICs, down to 0, indicating no interdependence. Trials were shuffled 200 times to examine the significance of the coherence values (p<0.05). The window length used in the coherence calculation was the same as those in calculation of the power spectrum. The ratio of the zero-pad to the windowed-period was 8. A three-way repeated-measures MANOVA was applied to compare mean coherence values among tasks, time durations, and frequency bands with Modality (crossmodal or unimodal), Duration (−100~0 ms, 0~100 ms, 100~200 ms, 200~300 ms, with the onset of S-1 as time 0) and Frequency-Band (Theta Band: 2~8 Hz; Alpha Band: 8~14 Hz; Early Beta Band: 14~20 Hz; Late Beta Band: 20~30 Hz) as the within-subject factors.

11.
Figure 1. Tasks and scalp electrode distributions. Upper-left: Schematic description of delayed matching-to-sample tasks. In the unimodal matching task, stimulus-1 (S-1) is a tactile vibration (150 Hz or 80 Hz) delivered on the subject's left index fingertip. Stimulus-2 (S-2) is also a tactile vibration. In the crossmodal matching task, S-2 is a light (red or green) from a light-emitting diode (LED) presented in front of the subject at eye level. The green light matches high frequency and the red light, low frequency. Upper-right: a top view of scalp electrode distributions. Nose and ears are shown in the diagram. Ag-AgCl electrodes are in a standard arrangement for locations. Lower: Grand average ERPs recorded in performance of the matching tasks. ERP components P45, P100, and N140 are indicated by arrows. The ERPs are time-locked to the onset of stimulus-1 (S-1).

12.
Figure 2. Topographic maps of an independent component (IC-F) located in frontal areas. Color-scale shows the value of the projection coefficient of the component. The topography of the IC-F is consistent across subjects (n = 10, indicated by numbers) and between tasks, unimodal (U) and crossmodal (C). Figure 3. Topographic maps of an independent component (IC-RS) located in right somatosensory areas.

13.
Figure 4. The average topography of IC-F (upper-left) and IC-RS (lower-left), and the corresponding BESA fitting dipole positions. The grand mean of the topographic maps is from 10 subjects across the tasks. Dipoles indicating the source of the components are located in medial prefrontal areas (IC-F, upper-right) and somatosensory areas (IC-RS, lower-right) respectively. Image views of the brain for each component are (clockwise from the top-left): sagittal (Sag), coronal (Cor), horizontal (Hori), and three-dimensional (3D). A: anterior; P: posterior; L: left; R: right.

14.
Figure 5. Comparisons of ERP components and IC projections. Left: Comparisons of latency and amplitude of ERP component P100 and IC-F projections. Middle: Grand averages of the original ERPs recorded from those midline electrodes, and grand average back-projections of the IC-F component to those electrodes. Right: Comparisons of latency and amplitude of ERP component N140 and IC-F projections. The percentage number indicates the proportion of potential that the IC-F contributes to the original ERP N140. C: crossmodal. U: unimodal. In bar graphs, the range of the ordinate for latency of P100 and N140 is 0–200 ms; the range of the ordinate for absolute values of amplitude of P100 and N140 is 0–7 uV. Error bars represent SEMs in this figure and other figures.

15.
Figure 6. Grand averages of the original ERPs recorded at the electrodes (C4, CP4, and P4) contralateral to the tactile stimulus, and also at those (C3, CP3, and P3) ipsilateral to it. ERP component P45 is shown at those contralateral electrodes. Grand average back-projections of the IC-RS component to those electrodes are also shown, where the projections have the largest peaks. Note those ERP P45 peaks and IC-RS back-projections are similar in both latency and amplitude. Figure 7. Grand average back-projections of the IC-RS. The significant difference in the projections at the electrodes between unimodal and crossmodal tasks are labeled with asterisks in three durations: 30~70 ms (yellow), 70~100 ms (gray), 100~160 ms (green).

16.
Figure 8. Time-frequency representation (TFR) for IC-F and the original ERPs at FCz (upper), and for IC-RS and the original ERPs at C4 (lower). Results are the average of all trials over 10 subjects and displayed in units of standard deviation of the baseline. Time zero is the onset of stimulus-1 in the tasks, crossmodal and unimodal. The peak frequency is indicated by a white square in each corresponding representation.

17.
Figure 9. Average coherence between IC-F and IC-RS across subjects. Upper : Results in the crossmodal task; Lower : Results in the unimodal task. Time-frequency representation of coherence index is shown on the left side for both tasks. The coherence index across different frequency bands during different time durations is shown on the right side for the tasks. Post hoc (Tukey HSD) test shows that the theta-band coherence during 100~200 ms is significantly different from the baseline (−100~0 ms) in the crossmodal task.

18.
Figure 10. Granger causality analysis between IC-F and IC-RS for the crossmodal task. Post hoc (Tukey HSD) test shows that in the theta band the bottom-up connectivity in the period of 100~200 ms is significantly stronger than that in the baseline (−100~0 ms). Top-down : Granger causality from IC-F to IC-RS; Bottom-up : Granger causality from IC-RS to IC-F.

19.
Discussion <ul><li>In the present study, the use of the ICA technique enabled us to find from the original EEG data two ICs (IC-F and IC-RS) that represented neural activities correlated with tactile working memory tasks, unimodal or crossmodal. </li></ul><ul><li>This finding strongly suggests that cortical locations of those two ICs, medial prefrontal cortex and SI are involved in perception of the tactile stimulus and crossmodal associations in the task, and it may therefore provide us with a better understanding of the neural mechanism underlying the crossmodal working memory. </li></ul>

20.
Discussion <ul><li>The theta oscillation in ACC may play an important role in the attentional modulation Our coherence analysis indicated that in the crossmodal task, compared with the baseline period the coherence in theta range during the period of 100~200 ms after the onset of the tactile stimulus was significantly increased between IC-F and IC-RS, showing that the activity in SI cortex may be synchronized with the activity in ACC in crossmodal association. </li></ul>

21.
Abstract Our previous studies on scalp-recorded event-related potentials (ERPs) showed that somatosensory N140 evoked by a tactile vibration in working memory tasks was enhanced when human subjects expected a coming visual stimulus that had been paired with the tactile stimulus. The results suggested that such enhancement represented the cortical activities involved in tactile-visual crossmodal association. In the present study, we further hypothesized that the enhancement represented the neural activities in somatosensory and frontal cortices in the crossmodal association. By applying independent component analysis (ICA) to the ERP data, we found independent components (ICs) located in the medial prefrontal cortex (around the anterior cingulate cortex, ACC) and the primary somatosensory cortex (SI). The activity represented by the IC in SI cortex showed enhancement in expectation of the visual stimulus. Such differential activity thus suggested the participation of SI cortex in the task-related crossmodal association. Further, the coherence analysis and the Granger causality spectral analysis of the ICs showed that SI cortex appeared to cooperate with ACC in attention and perception of the tactile stimulus in crossmodal association. The results of our study support with new evidence an important idea in cortical neurophysiology: higher cognitive operations develop from the modality-specific sensory cortices (in the present study, SI cortex) that are involved in sensation and perception of various stimuli.