Further characterisation of late somatosensory evoked potentials using electroencephalogram and magnetoencephalogram source imaging.

Beside the well-documented involvement of secondary somatosensory area, the cortical network underlying late somatosensory evoked potentials (P60/N60 and P100/N100) is still unknown. Electroencephalogram and magnetoencephalogram source imaging were performed to further investigate the origin of the brain cortical areas involved in late somatosensory evoked potentials, using sensory inputs of different strengths and by testing the correlation between cortical sources. Simultaneous high-density electroencephalograms and magnetoencephalograms were performed in 19 participants, and electrical stimulation was applied to the median nerve (wrist level) at intensity between 1.5 and 9 times the perceptual threshold. Source imaging was undertaken to map the stimulus-induced brain cortical activity according to each individual brain magnetic resonance imaging, during three windows of analysis covering early and late somatosensory evoked potentials. Results for P60/N60 and P100/N100 were compared with those for P20/N20 (early response). According to literature, maximal activity during P20/N20 was found in central sulcus contralateral to stimulation site. During P60/N60 and P100/N100, activity was observed in contralateral primary sensorimotor area, secondary somatosensory area (on both hemispheres) and premotor and multisensory associative cortices. Late responses exhibited similar characteristics but different from P20/N20, and no significant correlation was found between early and late generated activities. Specific clusters of cortical activities were activated with specific input/output relationships underlying early and late somatosensory evoked potentials. Cortical networks, partly common to and distinct from early somatosensory responses, contribute to late responses, all participating in the complex somatosensory brain processing.

Spatial functional reorganizations can serve as potential biomarkers of post facial palsy synkinesis.

Facial palsy can result in a serious complication known as facial synkinesis, causing both physical and psychological harm to the patients. There is growing evidence that patients with facial synkinesis have brain abnormalities, but the brain mechanisms and underlying imaging biomarkers remain unclear. Here, we employed functional magnetic resonance imaging (fMRI) to investigate brain function in 31 unilateral post facial palsy synkinesis patients and 25 healthy controls during different facial expression movements and at rest. Combining surface-based mass-univariate analysis and multivariate pattern analysis, we identified diffused activation and intrinsic connection patterns in the primary motor cortex and the somatosensory cortex on the patient's affected side. Further, we classified post facial palsy synkinesis patients from healthy subjects with favorable accuracy using the support vector machine based on both task-related and resting-state functional magnetic resonance imaging data. Together, these findings indicate the potential of the identified functional reorganizations to serve as neuroimaging biomarkers for facial synkinesis diagnosis.

Functional mapping of the somatosensory cortex using noninvasive fMRI and touch in awake dogs.

Dogs are increasingly used as a model for neuroscience due to their ability to undergo functional MRI fully awake and unrestrained, after extensive behavioral training. Still, we know rather little about dogs' basic functional neuroanatomy, including how basic perceptual and motor functions are localized in their brains. This is a major shortcoming in interpreting activations obtained in dog fMRI. The aim of this preregistered study was to localize areas associated with somatosensory processing. To this end, we touched N = 22 dogs undergoing fMRI scanning on their left and right flanks using a wooden rod. We identified activation in anatomically defined primary and secondary somatosensory areas (SI and SII), lateralized to the contralateral hemisphere depending on the side of touch, and importantly also activation beyond SI and SII, in the cingulate cortex, right cerebellum and vermis, and the sylvian gyri. These activations may partly relate to motor control (cerebellum, cingulate), but also potentially to higher-order cognitive processing of somatosensory stimuli (rostral sylvian gyri), and the affective aspects of the stimulation (cingulate). We also found evidence for individual side biases in a vast majority of dogs in our sample, pointing at functional lateralization of somatosensory processing. These findings not only provide further evidence that fMRI is suited to localize neuro-cognitive processing in dogs, but also expand our understanding of in vivo touch processing in mammals, beyond classically defined primary and secondary somatosensory cortices.

Separation of bimodal fMRI responses in mouse somatosensory areas into V1 and non-V1 contributions.

Multisensory integration is necessary for the animal to survive in the real world. While conventional methods have been extensively used to investigate the multisensory integration process in various brain areas, its long-range interactions remain less explored. In this study, our goal was to investigate interactions between visual and somatosensory networks on a whole-brain scale using 15.2-T BOLD fMRI. We compared unimodal to bimodal BOLD fMRI responses and dissected potential cross-modal pathways with silencing of primary visual cortex (V1) by optogenetic stimulation of local GABAergic neurons. Our data showed that the influence of visual stimulus on whisker activity is higher than the influence of whisker stimulus on visual activity. Optogenetic silencing of V1 revealed that visual information is conveyed to whisker processing via both V1 and non-V1 pathways. The first-order ventral posteromedial thalamic nucleus (VPM) was functionally affected by non-V1 sources, while the higher-order posterior medial thalamic nucleus (POm) was predominantly modulated by V1 but not non-V1 inputs. The primary somatosensory barrel field (S1BF) was influenced by both V1 and non-V1 inputs. These observations provide valuable insights for into the integration of whisker and visual sensory information.

Maternal prenatal depressive symptoms and child brain responses to affective touch at two years of age.

BACKGROUND: Touch is an essential form of mother-child interaction, instigating better social bonding and emotional stability. METHODS: We used diffuse optical tomography to explore the relationship between total haemoglobin (HbT) responses to affective touch in the child's brain at two years of age and maternal self-reported prenatal depressive symptoms (EPDS). Affective touch was implemented via slow brushing of the child's right forearm at 3 cm/s and non-affective touch via fast brushing at 30 cm/s and HbT responses were recorded on the left hemisphere. RESULTS: We discovered a cluster in the postcentral gyrus exhibiting a negative correlation (Pearson's r = -0.84, p = 0.015 corrected for multiple comparisons) between child HbT response to affective touch and EPDS at gestational week 34. Based on region of interest (ROI) analysis, we found negative correlations between child responses to affective touch and maternal prenatal EPDS at gestational week 14 in the precentral gyrus, Rolandic operculum and secondary somatosensory cortex. The responses to non-affective touch did not correlate with EPDS in these regions. LIMITATIONS: The number of mother-child dyads was 16. However, by utilising high-density optode arrangements, individualised anatomical models, and video and accelerometry to monitor movement, we were able to minimize methodological sources of variability in the data. CONCLUSIONS: The results show that maternal depressive symptoms during pregnancy may be associated with reduced child responses to affective touch in the temporoparietal cortex. Responses to affective touch may be considered as potential biomarkers for psychosocial development in children. Early identification of and intervention in maternal depression may be important already during early pregnancy.

The influence of hand posture on tactile processing: Evidence from a 7T functional magnetic resonance imaging study.

Although behavioral evidence has shown that postural changes influence the ability to localize or detect tactile stimuli, little is known regarding the brain areas that modulate these effects. This 7T functional magnetic resonance imaging (fMRI) study explores the effects of touch of the hand as a function of hand location (right or left side of the body) and hand configuration (open or closed). We predicted that changes in hand configuration would be represented in contralateral primary somatosensory cortex (S1) and the anterior intraparietal area (aIPS), whereas change in position of the hand would be associated with alterations in activation in the superior parietal lobule. Multivoxel pattern analysis and a region of interest approach partially supported our predictions. Decoding accuracy for hand location was above chance level in superior parietal lobule (SPL) and in the anterior intraparietal (aIPS) area; above chance classification of hand configuration was observed in SPL and S1. This evidence confirmed the role of the parietal cortex in postural effects on touch and the possible role of S1 in coding the body form representation of the hand.

fMRI, LFP, and anatomical evidence for hierarchical nociceptive routing pathway between somatosensory and insular cortices.

The directional organization of multiple nociceptive regions, particularly within obscure operculoinsular areas, underlying multidimensional pain processing remains elusive. This study aims to establish the fundamental organization between somatosensory and insular cortices in routing nociceptive information. By employing an integrated multimodal approach of high-field fMRI, intracranial electrophysiology, and transsynaptic viral tracing in rats, we observed a hierarchically organized connection of S1/S2 → posterior insula → anterior insula in routing nociceptive information. The directional nociceptive pathway determined by early fMRI responses was consistent with that examined by early evoked LFP, intrinsic effective connectivity, and anatomical projection, suggesting fMRI could provide a valuable facility to discern directional neural circuits in animals and humans non-invasively. Moreover, our knowledge of the nociceptive hierarchical organization of somatosensory and insular cortices and the interface role of the posterior insula may have implications for the development of targeted pain therapies.

Stimulus awareness is associated with secondary somatosensory cortex activation in an inattentional numbness paradigm.

While inattentional blindness and deafness studies have revealed neural correlates of consciousness (NCC) without the confound of task relevance in the visual and auditory modality, comparable studies for the somatosensory modality are lacking. Here, we investigated NCC using functional magnetic resonance imaging (fMRI) in an inattentional numbness paradigm. Participants (N = 44) received weak electrical stimulation on the left hand while solving a demanding visual task. Half of the participants were informed that task-irrelevant weak tactile stimuli above the detection threshold would be applied during the experiment, while the other half expected stimuli below the detection threshold. Unexpected awareness assessments after the experiment revealed that altogether 10 participants did not consciously perceive the somatosensory stimuli during the visual task. Awareness was not significantly modulated by prior information. The fMRI data show that awareness of stimuli led to increased activation in the contralateral secondary somatosensory cortex. We found no significant effects of stimulus awareness in the primary somatosensory cortex or frontoparietal areas. Thus, our results support the hypothesis that somatosensory stimulus awareness is mainly based on activation in higher areas of the somatosensory cortex and does not require strong activation in extended anterior or posterior networks, which is usually seen when perceived stimuli are task-relevant.

Brain Mapping Using a Graphene Electrode Array.

Cortical maps represent the spatial organization of location-dependent neural responses to sensorimotor stimuli in the cerebral cortex, enabling the prediction of physiologically relevant behaviors. Various methods, such as penetrating electrodes, electroencephalography, positron emission tomography, magnetoencephalography, and functional magnetic resonance imaging, have been used to obtain cortical maps. However, these methods are limited by poor spatiotemporal resolution, low signal-to-noise ratio (SNR), high costs, and non-biocompatibility or cause physical damage to the brain. This study proposes a graphene array-based somatosensory mapping method as a feature of electrocorticography that offers superior biocompatibility, high spatiotemporal resolution, desirable SNR, and minimized tissue damage, overcoming the drawbacks of previous methods. This study demonstrated the feasibility of a graphene electrode array for somatosensory mapping in rats. The presented protocol can be applied not only to the somatosensory cortex but also to other cortices such as the auditory, visual, and motor cortex, providing advanced technology for clinical implementation.

Reduced dimension stimulus decoding and column-based modeling reveal architectural differences of primary somatosensory finger maps between younger and older adults.

The primary somatosensory cortex (SI) contains fine-grained tactile representations of the body, arranged in an orderly fashion. The use of ultra-high resolution fMRI data to detect group differences, for example between younger and older adults' SI maps, is challenging, because group alignment often does not preserve the high spatial detail of the data. Here, we use robust-shared response modeling (rSRM) that allows group analyses by mapping individual stimulus-driven responses to a lower dimensional shared feature space, to detect age-related differences in tactile representations between younger and older adults using 7T-fMRI data. Using this method, we show that finger representations are more precise in Brodmann-Area (BA) 3b and BA1 compared to BA2 and motor areas, and that this hierarchical processing is preserved across age groups. By combining rSRM with column-based decoding (C-SRM), we further show that the number of columns that optimally describes finger maps in SI is higher in younger compared to older adults in BA1, indicating a greater columnar size in older adults' SI. Taken together, we conclude that rSRM is suitable for finding fine-grained group differences in ultra-high resolution fMRI data, and we provide first evidence that the columnar architecture in SI changes with increasing age.

Rapid Processing of Observed Touch through Social Perceptual Brain Regions: An EEG-fMRI Fusion Study.

Seeing social touch triggers a strong social-affective response that involves multiple brain networks, including visual, social perceptual, and somatosensory systems. Previous studies have identified the specific functional role of each system, but little is known about the speed and directionality of the information flow. Is this information extracted via the social perceptual system or from simulation from somatosensory cortex? To address this, we examined the spatiotemporal neural processing of observed touch. Twenty-one human participants (seven males) watched 500-ms video clips showing social and nonsocial touch during electroencephalogram (EEG) recording. Visual and social-affective features were rapidly extracted in the brain, beginning at 90 and 150 ms after video onset, respectively. Combining the EEG data with functional magnetic resonance imaging (fMRI) data from our prior study with the same stimuli reveals that neural information first arises in early visual cortex (EVC), then in the temporoparietal junction and posterior superior temporal sulcus (TPJ/pSTS), and finally in the somatosensory cortex. EVC and TPJ/pSTS uniquely explain EEG neural patterns, while somatosensory cortex does not contribute to EEG patterns alone, suggesting that social-affective information may flow from TPJ/pSTS to somatosensory cortex. Together, these findings show that social touch is processed quickly, within the timeframe of feedforward visual processes, and that the social-affective meaning of touch is first extracted by a social perceptual pathway. Such rapid processing of social touch may be vital to its effective use during social interaction.SIGNIFICANCE STATEMENT Seeing physical contact between people evokes a strong social-emotional response. Previous research has identified the brain systems responsible for this response, but little is known about how quickly and in what direction the information flows. We demonstrated that the brain processes the social-emotional meaning of observed touch quickly, starting as early as 150 ms after the stimulus onset. By combining electroencephalogram (EEG) data with functional magnetic resonance imaging (fMRI) data, we show for the first time that the social-affective meaning of touch is first extracted by a social perceptual pathway and followed by the later involvement of somatosensory simulation. This rapid processing of touch through the social perceptual route may play a pivotal role in effective usage of touch in social communication and interaction.

Two-Dimensional Population Receptive Field Mapping of Human Primary Somatosensory Cortex.

Functional magnetic resonance imaging can provide detailed maps of how sensory space is mapped in the human brain. Here, we use a novel 16 stimulator setup (a 4 × 4 grid) to measure two-dimensional sensory maps of between and within-digit (D2-D4) space using high spatial-resolution (1.25 mm isotropic) imaging at 7 Tesla together with population receptive field (pRF) mapping in 10 participants. Using a 2D Gaussian pRF model, we capture maps of the coverage of digits D2-D5 across Brodmann areas and estimate pRF size and shape. In addition, we compare results to previous studies that used fewer stimulators by constraining pRF models to a 1D Gaussian Between Digit or 1D Gaussian Within Digit model. We show that pRFs across somatosensory areas tend to have a strong preference to cover the within-digit axis. We show an increase in pRF size moving from D2-D5. We quantify pRF shapes in Brodmann area (BA) 3b, 3a, 1, 2 and show differences in pRF size in Brodmann areas 3a-2, with larger estimates for BA2. Generally, the 2D Gaussian pRF model better represents pRF coverage maps generated by our data, which itself is produced from a 2D stimulation grid.

Pneumatic tactile stimulus delivery system for studying brain responses evoked by active finger touch with fMRI.

BACKGROUND: Primates use their hands to actively touch objects and collect information. To study tactile information processing, it is important for participants to experience tactile stimuli through active touch while monitoring brain activities. NEW METHOD: Here, we developed a pneumatic tactile stimulus delivery system (pTDS) that delivers various tactile stimuli on a programmed schedule and allows voluntary finger touches during MRI scanning. The pTDS uses a pneumatic actuator to move tactile stimuli and place them in a finger hole. A photosensor detects the time when an index finger touches a tactile stimulus, enabling the analysis of the touch-elicited brain responses. RESULTS: We examined brain responses while the participants actively touched braille objects presented by the pTDS. BOLD responses during tactile perception were significantly stronger in a finger touch area of the contralateral somatosensory cortex compared with that of visual perception. CONCLUSION: The pTDS enables MR studies of brain mechanisms for tactile processes through natural finger touch.

Transcranial focused ultrasound stimulation of cortical and thalamic somatosensory areas in human.

The effects of transcranial focused ultrasound (FUS) stimulation of the primary somatosensory cortex and its thalamic projection (i.e., ventral posterolateral nucleus) on the generation of electroencephalographic (EEG) responses were evaluated in healthy human volunteers. Stimulation of the unilateral somatosensory circuits corresponding to the non-dominant hand generated EEG evoked potentials across all participants; however, not all perceived stimulation-mediated tactile sensations of the hand. These FUS-evoked EEG potentials (FEP) were observed from both brain hemispheres and shared similarities with somatosensory evoked potentials (SSEP) from median nerve stimulation. Use of a 0.5 ms pulse duration (PD) sonication given at 70% duty cycle, compared to the use of 1 and 2 ms PD, elicited more distinctive FEP peak features from the hemisphere ipsilateral to sonication. Although several participants reported hearing tones associated with FUS stimulation, the observed FEP were not likely to be confounded by the auditory sensation based on a separate measurement of auditory evoked potentials (AEP) to tonal stimulation (mimicking the same repetition frequency as the FUS stimulation). Off-line changes in resting-state functional connectivity (FC) associated with thalamic stimulation revealed that the FUS stimulation enhanced connectivity in a network of sensorimotor and sensory integration areas, which lasted for at least more than an hour. Clinical neurological evaluations, EEG, and neuroanatomical MRI did not reveal any adverse or unintended effects of sonication, attesting its safety. These results suggest that FUS stimulation may induce long-term neuroplasticity in humans, indicating its neurotherapeutic potential for various neurological and neuropsychiatric conditions.

Artificial gravity during a spaceflight analog alters brain sensory connectivity.

Spaceflight has numerous untoward effects on human physiology. Various countermeasures are under investigation including artificial gravity (AG). Here, we investigated whether AG alters resting-state brain functional connectivity changes during head-down tilt bed rest (HDBR), a spaceflight analog. Participants underwent 60 days of HDBR. Two groups received daily AG administered either continuously (cAG) or intermittently (iAG). A control group received no AG. We assessed resting-state functional connectivity before, during, and after HDBR. We also measured balance and mobility changes from pre- to post-HDBR. We examined how functional connectivity changes throughout HDBR and whether AG is associated with differential effects. We found differential connectivity changes by group between posterior parietal cortex and multiple somatosensory regions. The control group exhibited increased functional connectivity between these regions throughout HDBR whereas the cAG group showed decreased functional connectivity. This finding suggests that AG alters somatosensory reweighting during HDBR. We also observed brain-behavioral correlations that differed significantly by group. Control group participants who showed increased connectivity between the putamen and somatosensory cortex exhibited greater mobility declines post-HDBR. For the cAG group, increased connectivity between these regions was associated with little to no mobility declines post-HDBR. This suggests that when somatosensory stimulation is provided via AG, functional connectivity increases between the putamen and somatosensory cortex are compensatory in nature, resulting in reduced mobility declines. Given these findings, AG may be an effective countermeasure for the reduced somatosensory stimulation that occurs in both microgravity and HDBR.

Fine-grained topographic organization within somatosensory cortex during resting-state and emotional face-matching task and its association with ASD traits.

Sensory atypicalities are particularly common in autism spectrum disorders (ASD). Nevertheless, our knowledge about the divergent functioning of the underlying somatosensory region and its association with ASD phenotype features is limited. We applied a data-driven approach to map the fine-grained variations in functional connectivity of the primary somatosensory cortex (S1) to the rest of the brain in 240 autistic and 164 neurotypical individuals from the EU-AIMS LEAP dataset, aged between 7 and 30. We estimated the S1 connection topography ('connectopy') at rest and during the emotional face-matching (Hariri) task, an established measure of emotion reactivity, and accessed its association with a set of clinical and behavioral variables. We first demonstrated that the S1 connectopy is organized along a dorsoventral axis, mapping onto the S1 somatotopic organization. We then found that its spatial characteristics were linked to the individuals' adaptive functioning skills, as measured by the Vineland Adaptive Behavior Scales, across the whole sample. Higher functional differentiation characterized the S1 connectopies of individuals with higher daily life adaptive skills. Notably, we detected significant differences between rest and the Hariri task in the S1 connectopies, as well as their projection maps onto the rest of the brain suggesting a task-modulating effect on S1 due to emotion processing. All in all, variation of adaptive skills appears to be reflected in the brain's mesoscale neural circuitry, as shown by the S1 connectivity profile, which is also differentially modulated during rest and emotional processing.

Content Representation of Tactile Mental Imagery in Primary Somatosensory Cortex.

The imagination of tactile stimulation has been shown to activate primary somatosensory cortex (S1) with a somatotopic specificity akin to that seen during the perception of tactile stimuli. Using fMRI and multivariate pattern analysis, we investigate whether this recruitment of sensory regions also reflects content-specific activation (i.e., whether the activation in S1 is specific to the mental content participants imagined). To this end, healthy volunteers (n = 21) either perceived or imagined three types of vibrotactile stimuli (mental content) while fMRI data were acquired. Independent of the content, during tactile mental imagery we found activation of frontoparietal regions, supplemented with activation in the contralateral BA2 subregion of S1, replicating previous reports. While the imagery of the three different stimuli did not reveal univariate activation differences, using multivariate pattern classification, we were able to decode the imagined stimulus type from BA2. Moreover, cross-classification revealed that tactile imagery elicits activation patterns similar to those evoked by the perception of the respective stimuli. These findings promote the idea that mental tactile imagery involves the recruitment of content-specific activation patterns in sensory cortices, namely in S1.

Similar somatotopy for active and passive digit representation in primary somatosensory cortex.

Scientists traditionally use passive stimulation to examine the organisation of primary somatosensory cortex (SI). However, given the close, bidirectional relationship between the somatosensory and motor systems, active paradigms involving free movement may uncover alternative SI representational motifs. Here, we used 7 Tesla functional magnetic resonance imaging to compare hallmark features of SI digit representation between active and passive tasks which were unmatched on task or stimulus properties. The spatial location of digit maps, somatotopic organisation, and inter-digit representational structure were largely consistent between tasks, indicating representational consistency. We also observed some task differences. The active task produced higher univariate activity and multivariate representational information content (inter-digit distances). The passive task showed a trend towards greater selectivity for digits versus their neighbours. Our findings highlight that, while the gross features of SI functional organisation are task invariant, it is important to also consider motor contributions to digit representation.

Talking with hands and feet: Selective somatosensory attention and fMRI enable robust and convenient brain-based communication.

In brain-based communication, voluntarily modulated brain signals (instead of motor output) are utilized to interact with the outside world. The possibility to circumvent the motor system constitutes an important alternative option for severely paralyzed. Most communication brain-computer interface (BCI) paradigms require intact visual capabilities and impose a high cognitive load, but for some patients, these requirements are not given. In these situations, a better-suited, less cognitively demanding information-encoding approach may exploit auditorily-cued selective somatosensory attention to vibrotactile stimulation. Here, we propose, validate and optimize a novel communication-BCI paradigm using differential fMRI activation patterns evoked by selective somatosensory attention to tactile stimulation of the right hand or left foot. Using cytoarchitectonic probability maps and multi-voxel pattern analysis (MVPA), we show that the locus of selective somatosensory attention can be decoded from fMRI-signal patterns in (especially primary) somatosensory cortex with high accuracy and reliability, with the highest classification accuracy (85.93%) achieved when using Brodmann area 2 (SI-BA2) at a probability level of 0.2. Based on this outcome, we developed and validated a novel somatosensory attention-based yes/no communication procedure and demonstrated its high effectiveness even when using only a limited amount of (MVPA) training data. For the BCI user, the paradigm is straightforward, eye-independent, and requires only limited cognitive functioning. In addition, it is BCI-operator friendly given its objective and expertise-independent procedure. For these reasons, our novel communication paradigm has high potential for clinical applications.

Identifying the distinct spectral dynamics of laminar-specific interhemispheric connectivity with bilateral line-scanning fMRI.

Despite extensive efforts to identify interhemispheric functional connectivity (FC) with resting-state (rs-) fMRI, correlated low-frequency rs-fMRI signal fluctuation across homotopic cortices originates from multiple sources. It remains challenging to differentiate circuit-specific FC from global regulation. Here, we developed a bilateral line-scanning fMRI method to detect laminar-specific rs-fMRI signals from homologous forepaw somatosensory cortices with high spatial and temporal resolution in rat brains. Based on spectral coherence analysis, two distinct bilateral fluctuation spectral features were identified: ultra-slow fluctuation (<0.04 Hz) across all cortical laminae versus Layer (L) 2/3-specific evoked BOLD at 0.05 Hz based on 4 s on/16 s off block design and resting-state fluctuations at 0.08-0.1 Hz. Based on the measurements of evoked BOLD signal at corpus callosum (CC), this L2/3-specific 0.05 Hz signal is likely associated with neuronal circuit-specific activity driven by the callosal projection, which dampened ultra-slow oscillation less than 0.04 Hz. Also, the rs-fMRI power variability clustering analysis showed that the appearance of L2/3-specific 0.08-0.1 Hz signal fluctuation is independent of the ultra-slow oscillation across different trials. Thus, distinct laminar-specific bilateral FC patterns at different frequency ranges can be identified by the bilateral line-scanning fMRI method.