Anterior insular cortex plays a critical role in interoceptive attention

Accumulating evidence indicates that the anterior insular cortex (AIC) mediates interoceptive attention which refers to attention towards physiological signals arising from the body. However, the necessity of the AIC in this process has not been demonstrated. Using a novel task that directs attention toward breathing rhythm, we assessed the involvement of the AIC in interoceptive attention in healthy participants using functional magnetic resonance imaging and examined the necessity of the AIC in interoceptive attention in patients with AIC lesions. Results showed that interoceptive attention was associated with increased AIC activation, as well as enhanced coupling between the AIC and somatosensory areas along with reduced coupling between the AIC and visual sensory areas. In addition, AIC activation was predictive of individual differences in interoceptive accuracy. Importantly, AIC lesion patients showed disrupted interoceptive discrimination accuracy and sensitivity. These results provide compelling evidence that the AIC plays a critical role in interoceptive attention.


Introduction
Our brain consistently receives physiological signals arising from sensory inputs and our body. Although attention toward inputs from the external environment (i.e., exteroceptive attention) has been extensively investigated, the attentional mechanism of the awareness and conscious focus on bodily somatic and visceral signals or responses (i.e., interoceptive attention) has been less studied because of difficulties in its measurement (Brener and Ring, 2016;Ring et al., 2015;Craig, 2002;Craig, 2003;Craig, 2010;Critchley, 2005;Critchley, 2004;Farb et al., 2013a). Previous theories argue that subjective emotions arise from these bodily reactions and visceral experiences (Cannon, 1987;Critchley and Harrison, 2013;Damasio, 1996;Dolan, 2002;Tranel and Damasio, 1991) in which interoceptive awareness plays a critical role. Appropriate attention to bodily states and accurate perception of interoceptive information are essential in emotional awareness and in the maintenance of normal physiological conditions (Craig, 2002;Craig, 2003;Craig, 2010;Critchley, 2005;Wiens, 2005). The link between deficits in interoceptive attention and psychiatric symptoms may also be explained by the James-Lange theory of emotion (Cannon, 1987), the somatic marker hypothesis (Damasio, 1996;Damasio et al., 1991) for the embodied mind mediated by interoception (Garfinkel and Critchley, 2013), and the embodied predictive processing model (Allen et al., 2016;Allen and Friston, 2018;Barrett and Simmons, 2015;Seth, 2013;Seth and Critchley, 2013;Seth et al., 2011).
Recent human studies have emphasized the role of the insula in interoceptive representations (Daubenmier et al., 2013;Farb et al., 2013b;Ronchi et al., 2015). Neuroanatomical evidence, consistent with neuroimaging findings, suggests that the anterior insular cortex (AIC) is an important structure for encoding and representing interoceptive information (Craig, 2002;Craig, 2003;Craig, 2009;Critchley et al., 2004;Stephani et al., 2011). Although the AIC has been recognized as an interoceptive cortex (Craig, 2003;Critchley et al., 2004;Ernst et al., 2014;Singer et al., 2009;Terasawa et al., 2013), these findings remain equivocal because AIC activation seems ubiquitous across a wide range of tasks involving cognition, emotion, and other cognitive processes in addition to interoceptive attention (Allen et al., 2013;Seeley et al., 2007;Uddin et al., 2014;Yarkoni et al., 2011). Therefore, a task that selectively and reliably engages interoceptive attention needs to be employed. In addition, the correlational AIC activation found in functional neuroimaging studies alone does not provide causal evidence for its role in interoceptive attention, leaving the question of whether the AIC is critical in interoceptive attention unanswered. Studying patients with focal lesions in the AIC (Gu et al., 2012;Gu et al., 2015;Ronchi et al., 2015;Starr et al., 2009;Wang et al., 2014;Wu et al., 2019) would thus provide a unique opportunity to examine the necessity of the AIC in this fundamental process.
One challenge to the study of interoceptive attention is the vague nature of interoceptive awareness. According to the classic definition of attention by James (1890), only the contents that are clearly perceived and represented by the mind can be the target of attention. However, most existing tasks measuring interoceptive attention fail to meet this criterion (Ring et al., 2015). In contrast to exteroceptive attention toward external sensory inputs, precise measurements of interoceptive attention are difficult to obtain experimentally because of the imprecise perception of visceral changes, such as heart rate (Brener and Ring, 2016;Paulus and Stein, 2010;Ring et al., 2015;Windmann et al., 1999). Multiple sources of physical information contribute to bodily signals, and most of these sources of somatic feedback cannot be described accurately by mindful introspection in normal physiological states (Ring et al., 2015). This limitation impedes accurate measurement of interoceptive attention and examination of the neural mechanisms underlying this process. To overcome this barrier, a perceivable visceral channel needs to be used.
Breathing is an essential activity for maintaining human life and, more importantly, is an easily perceivable bodily signal. As an autonomous vital movement, breathing can be measured and actively controlled in humans (Daubenmier et al., 2013;Davenport et al., 2007). The unique physiological characteristics of respiration render breath detection an ideal method for measuring interoceptive accuracy and sensitivity (Garfinkel et al., 2015) and for exploring the neural activity underlying interoceptive attention. Thus, we designed a breath detection task to engage interoceptive attention (attention to bodily signals), in which participants were required to indicate whether a presented breathing curve is delayed or not relative to their own breathing rhythm (breath detection task, BDT), in contrast to engaging exteroceptive attention (attention to visual signals), in which participants were required to indicate whether a visual dot stimulus is flashed on the breathing curve (dot flash detection task, DDT). This design enabled us to examine the involvement of the AIC in interoceptive processing in healthy participants and the necessity of the AIC in this processing in patients with AIC lesions.
Basing from previous evidence (e.g., Critchley, 2004), we hypothesized that the AIC is critical for interoceptive attention to reach subjective awareness by integrating information from an individual's homeostatic state and the external environment. We first conducted functional magnetic resonance imaging (fMRI) studies with two samples to map the neural substrates underlying interoceptive attention to internal bodily signals in contrast to exteroceptive attention to external visual signals in healthy participants while they performed the tasks. We then investigated the necessity of the AIC in interoceptive attention by assessing interoceptive attention in patients with focal AIC lesions (AIC group) in comparison to brain-damaged controls (BDC group, patients with lesions in areas other than insular-or somatosensory-related cortices) and matched neurologically intact normal controls (NC group). We predicted that the AIC is involved in interoceptive attention and that patients with AIC lesions would show deficits in performance on the interoceptive but not exteroceptive attention task.

Results
Behavioral results of the fMRI studies Performance accuracy (%) and discrimination sensitivity (d') in the BDT were 82.1 ± 14.7% and 2.2 ± 1.1 (mean ± SD) for the first sample, and 74.9 ± 9.6% and 1.6 ± 0.6 (mean ± SD) for the second sample, respectively, which were significantly above the chance levels (50% and 0 for accuracy and d', respectively; For the first sample: t(43) = 14.51, p<0.001, Cohen's d = 2.18 for accuracy and t(43) = 13.09, p<0.001, Cohen's d = 2.0 for d', respectively; For the second sample: t(27) = 13.77, p<0.001, Cohen's d = 2.59 for accuracy and t(27) = 12.89, p<0.001, Cohen's d = 2.67 for d', respectively), but lower than the DDT accuracy of 87.3 ± 9.8% and d' of 2.6 ± 0.8 for the first sample (t (43) Figure 1-figure supplements 1 and 2 for details of the behavior results for the first and second samples, respectively, in accuracy, RT, d', and b. Data were plotted in R using 'raincloud' script (Allen et al., 2018a;Allen et al., 2018b); See Table 1 for the statistics of behavioral results for the first and second samples). The split-half reliability of the BDT and DDT were 0.86 and 0.85 for the first sample, and 0.68 and 0.89 for the second sample, respectively.
For the first sample, the relative interoceptive accuracy was negatively correlated with the subjectively scored difficulty of the BDT relative to the DDT (Pearson r = À0.43, corrected p=0.02, Bayes Factor (BF) = 10.38), but not significantly correlated with the 'awareness of bodily processes' subtest of the body perception questionnaire (BPQ) after correction for multiple comparisons (Pearson r = 0.27, corrected p=0.38, BF = 0.86). No significant correlations were observed between relative interoceptive accuracy and subjective emotion experiences, including trait positive affective experience (measured by Positive and Negative Affect Schedule, PANAS) (Watson, 1988) (Pearson r = 0.31, corrected p=0.20, BF = 1.38), anxiety (Pearson r = À0.006, p>0.9, BF = 0.19) or depression score (Pearson r = À0.002, p>0.9, BF = 0.19). For the second sample, however, we did not find significant correlations between relative interoceptive accuracy and scores of questionnaires (awareness of bodily processes: Pearson r = À0.17, corrected p>0.9, BF = 0.33; trait positive affective experience: Pearson r = 0.12, corrected p>0.9, BF = 0.27; anxiety: Pearson r = 0.29, corrected p=0.56, BF = 0.69; depression: Pearson r = 0.03, corrected p>0.9; note that we did not collect subjective rating of task difficulty in the second sample). In addition, we also calculated correlation coefficients between task performance and questionnaires by pooling the two samples (See Table 2 for Pearson correlation strength and Bayesian tests between all behavioral measures in the first sample, the second sample, and across the two samples).
Imaging results of the whole brain analysis of the first fMRI study Main effects of interoceptive attention and feedback delay, and their interaction The main effect of interoceptive attention, compared to exteroceptive attention (BDT vs. DDT), was associated with enhanced activity in the cognitive control network (Fan, 2014;Wu et al., 2015;Xuan et al., 2016), including the AIC, the dorsal anterior cingulate cortex (ACC) and the supplementary motor area (SMA), and the superior frontal and the parietal cortices (the frontal eye field, FEF; and the areas near/along the intraparietal sulcus; Figure 2a, Table 3). In addition, this contrast revealed significantly less activation, or deactivation, in the core regions of the default mode network (Raichle et al., 2001), including the the ventral medial prefrontal cortex, the middle temporal gyrus (MTG), and the posterior cingulate cortex.
Activation in the AIC, the middle frontal gyrus, the SMA, and the temporal parietal junction was associated with the main effect of feedback delay ( Figure 2b, Table 4). The regions showing the main effect of feedback delay also showed the interaction effect between attentional focus (interoceptively in the BDT and exteroceptively in the DDT) and feedback (with and without delay) ( Figure 2c, Table 5). The task-induced responses extracted from the bilateral AIC, defined by the attention by feedback interaction map, shows the activation pattern under different task conditions (Figure 2d). In specific, (1) the bilateral AIC demonstrated greater activation during interoceptive processing (BDT) than during exteroceptive processing (DDT), irrespective of the feedback type; (2) the delayed trials induced greater activation in the bilateral AIC in comparison to the non-delayed trials only during interoceptive processing (BDT). The evidence of this interaction effect in the AIC suggests that the AIC was actively engaged in interoceptive processing.  Figure 1. Experimental setup, trial structure of the tasks, and stimulus conditions. (a) The respiratory effort is converted into electronic signal changes using a respiratory transducer, amplified by BIOPAC, digitized using an A/D converter, and sent to the test computer for the final visual display as a dynamic breath curve, with or without a 400 ms delay. (b) This panel shows two trials for the breath detection task (BDT) and flash dot detection task (DDT) runs, respectively. Each trial begins with a 3 s blank display, followed by a 12 s display of respiratory curve presented with or without a 400 ms delay and with or without a 30 ms red dot flashed at a random position on the curve, and ends with a 3 s response window during which participants make a forced-choice button-press response to two alternative choices depending on the block type (BDT or DDT) to indicate whether the feedback curve is synchronous or delayed (for the BDT run) or whether a dot has appeared (for the DDT run). (c) The task represents a 2 Â 2 Â 2 factorial design with the factors of attention to breath or dot (block design), presence or absence of breath curve delay, and presence or absence of a dot flashed. The dashed line represents the actual breath curve, while the solid line represents the feedback breath curve displayed on the screen. DOI: https://doi.org/10.7554/eLife.42265.002 The following source data and figure supplements are available for figure 1:

Correlation between interoceptive accuracy and AIC activation
Voxel-wise regression analysis revealed the relationship between the interoceptive task-induced activation strength (map of the interaction contrast) and participants' interoceptive accuracy (performance accuracy in the BDT), with exteroceptive accuracy (performance accuracy in the DDT) controlled as a covariate. Higher interoceptive accuracy was associated with greater interaction effect of the bilateral AIC (and MTG) across participants ( Figure 3a, Table 6). The AIC activation during the interoceptive processing involved attending to physiological signals and matching bodily signals to external visual input, which predicted individual differences in interoceptive attention (see Figure 3b for the illustration).
Functional and effective connectivity of the AIC Psychophysiological interaction (PPI) analysis showed augmented connectivity between the right AIC (as the seed) and the SMA/ACC, the FEF, the inferior frontal gyrus (IFG), and the postcentral gyrus (PoCG) during interoceptive (versus exteroceptive) attention (BDT vs. DDT) in contrast to the reduced connectivity between the right AIC and visual cortices (VCs) modulated by interoceptive attention (Figure 4a, Table 7). This result indicates that an increase in activation in the right AIC was associated with a greater increase in activation in the FEF, the IFG, and the PoCG and a greater decrease in activation in the VCs under interoceptive attention compared with exteroceptive Here we showed the left AIC for the visualization of the seed for the ROI analysis in the second fMRI sample, although the cluster with 210 voxels did not survive GRF correction. Red color represents an increased activation; Blue color represents a decreased activation. (d) Activation of the left and the right AIC under the four task conditions, and the pattern of the interaction. DOI: https://doi.org/10.7554/eLife.42265. 010 The following source data is available for figure 2: Source data 1. CSV file containing data for Figure 2d. DOI: https://doi.org/10.7554/eLife.42265.011 Table 3. Activation and deactivation of the brain regions involved in interoceptive attention (interoception -exteroception Table 3 continued on next page attention ( Figure 4b). Similar PPI results were obtained when the left AIC was used as the seed ( On the basis of the PPI results, VCs of the right V2/3 (x = 14, y = À90, z = 28 as indicated by negative PPI) and the right PoCG (x = 58, y = À16, z = 32 as indicated by positive PPI) were included in the dynamic causal modeling (DCM) model. Data from one participant were excluded because significant activation in the V2/3 region of interest could not be identified. For model comparison, random-effects (RFX) Bayesian model selection (BMS) indicated that the winning model (with an exceedance probability of 29.84%) was the one with the modulatory effects of interoceptive and  Table 4. Activation and deactivation of the brain regions involved in feedback delay (delay -non-delay). exteroceptive attention (BDT and DDT) exerting on the connection from the AIC to the PoCG and from the AIC to V2/3 ( Figure 4c and Figure 4-figure supplement 2). The BMS indicated that interoceptive and exteroceptive attention were achieved through modulating the top-down connectivity from the AIC to these two sensory cortices. We performed parameter inference by using Bayesian model averaging (BMA), which considers uncertainty by pooling information across all models in a weighted fashion . For BMA (Figure 4d), the modulatory effect of interoceptive attention (BDT) was significant on the connection from the AIC to the PoCG (t(42) = 4.85, Bonferroni corrected p<0.001). The modulatory effect of exteroceptive attention (DDT) on the connection from the AIC to the V2/3 was significant without correction (t(42) = 2.25, uncorrected p=0.03). The BMA results were consistent with the winning model selected by model comparison and the PPI results: the modulatory effect from the AIC to the PoCG was driven by interoceptive attention (BDT), whereas the modulatory effect from the AIC to the V2/3 was driven by exteroceptive attention (DDT). In addition, the BMA results highlighted the importance of the intrinsic efferent connection from the AIC to the PoCG in the network (t(42) = 3.61, Bonferroni corrected p=0.01).

Region-of-interest (ROI) analysis results of the second fMRI study
The interaction between attentional focus (interoceptively in BDT and exteroceptively in DDT) and feedback (with and without delay) was significant in both left and right AICs (left: F(1, 27)=6.12, p=0.020; right: F(1,27) = 5.88, p=0.022; Figure 5a), which confirmed the interaction effect in the bilateral AIC revealed by whole brain analyses of the first sample. The main effect of attentional focus (BDT vs. DDT) was significant in the right AIC with greater activation during the BDT than during the DDT (F(1, 27)=4.20, p=0.05) but not significant in the left AIC (F(1, 27)<1, p=0.51). The main effect of the feedback was not significant in either left or right AIC (left: F<1; right: F<1). In addition, similar to the results of the first sample, we found a significant correlation between the interaction effect of both left and right AICs and relative interoceptive accuracy (left: Pearson r = 0.32, p=0.050, one-tailed; right: Pearson r = 0.42, p=0.014, one-tailed; Figure 5b). In addition, we examined the pattern of the respiratory volume under BDT and DDT (see Figure 5-figure supplement 1). Despite the difference in the respiratory volume between interoceptive and exteroceptive conditions (BDT and DDT was significant, F(1,27) = 15.88, p<0.001), this difference was canceled out for the interaction effect (F < 1). These results further illustrated that the interaction effect in the AIC is not subject to the confounding of breathing effort difference between the two tasks.
The whole brain analysis of the second fMRI sample showed significant overlap between the activations without and with physiological correction for the main and the interaction effects (see Figure 5-figure supplement 2). We further checked how much physiological noise impacted AIC activation by comparing the contrast maps without and with physiological correction at an extremely permissive threshold (p<0.05 uncorrected). The difference in signals of the AIC between the analyses with and without physiological corrections was only evident for the main effect of interoceptive vs. exteroceptive attention (BDT vs. DDT) but not for the interaction contrast, confirming that the . Altogether, these ROI results from the second sample confirmed that the AIC was actively engaged in interoceptive processing.
Lesion study results: the necessity of the AIC in interoceptive attention Figure 6 shows the insular lesion overlap for the AIC patient group. The area with the most overlap was identified as the AIC according to the literature (Kurth et al., 2010;Naidich et al., 2004). We found a significant interaction effect between group (AIC, BDC, and NC) and task (BDT and DDT) in performance accuracy (F (2,21) = 5.19, p=0.015) and discrimination sensitivity (d') (F (2,21) = 4.77, p=0.023). Planned simple comparisons were conducted between groups for each task. For the BDT, patients with AIC lesions had significantly lower performance accuracy (58%, t (13)  A summary of the statistical results of the lesion study is provided in Table 8. Our results demonstrated significant impairment in discrimination ability when attending to bodily signals, but not to external visual input, in patients with AIC lesions.

Discussion
Using fMRI, we showed that the AIC is involved in interoceptive attention towards respiration, with the underlying connectivity between the AIC and the somatosensory cortex and visual areas modulated by interoceptive and exteroceptive attention, respectively. Notably, we confirmed the necessity of the AIC in supporting interoceptive attention by showing reduced behavioral performance on the interoceptive task in patients with focal AIC lesions. Thus, this study demonstrates that the AIC plays a critical role in interoceptive attention.

The necessity of the AIC in interoceptive attention
Previous functional neuroimaging studies have shown that the insula is activated by autonomic arousal and emotional reactions (Craig, 2002;Craig, 2003;Critchley et al., 2004) and emphasized the central role of the insula in interoceptive awareness. The achievement of interoceptive awareness depends on the integration of afferent bodily signals with higher-order contextual information attributable to the AIC (Craig, 2002;Craig, 2009;Critchley, 2005;Damasio et al., 2000;Mutschler et al., 2009). In this study, the increase in neural activation in the AIC and other related brain structures when focusing on breath rhythm indicates that the AIC supports attention toward bodily signals. Most importantly, participants' performance accuracy on the interoceptive task was significantly correlated with the activation of the AIC, further demonstrating the involvement of the AIC in interoceptive attention. Anatomically, the insula receives thalamo-insular projections of the interoceptive pathways (Craig, 2002). The AIC encodes subjective feelings (Craig, 2003;Craig, 2009;Flynn, 1999) and is critical for instantaneous representation of the state of the body (Allen et al., 2016;Allen and Friston, 2018;Cao et al., 2014;Gu et al., 2015). During the BDT, this is achieved by attending to bodily signals (i.e., breath rhythm) and matching them to external visual feedback (i.e., the breath curve). The present results provide additional support to previous finding that the activation of the  right AIC is related to accuracy in sensing the timing of one's bodily signals, for example, heartbeat . Consistent with the notion that the AIC contributes to accurate perception of bodily states (Bechara and Naqvi, 2004), the insula works as a hub to convey bodily information into internal feelings for maintaining homeostasis and to mediate the representations of visceral states that link to the representations of the external world (Farb et al., 2013a).
Our finding of a critical role of the AIC in interoceptive attention fits with a recent predictive coding account of the brain (Bastos et al., 2012;Friston and Kiebel, 2009;Rao and Ballard, 1999), which suggests that the brain actively tries to predict possible future states and to minimize the difference between actual and predicted states. In the context of interoceptive and embodied predictive coding (Allen and Friston, 2018;Gu et al., 2013), previous studies hypothesized that interoceptive predictions are computed within a network of brain regions with the AIC as the key node (Allen et al., 2016;Barrett and Simmons, 2015;Seth, 2013). Empirical evidence that directly supports this computational role of the insula is still rare. One such study using a tactile oddball paradigm and DCM of fMRI time series demonstrated that the AIC is the only region, among a network of body-related brain regions, that shows a reciprocal increase in connectivity with the somatosensory cortex (Allen et al., 2016). Our finding is consistent with this previous study and extends the role of the AIC in predictive coding to breathing-related interoception by using both fMRI and lesion approaches.
The role of the AIC in interoceptive attention identified by the fMRI studies was augmented by the data from patients with focal damage to the insula. Relative to non-insular lesion patients and healthy controls, AIC lesions led to a deficit in accuracy and sensitivity of interoceptive attention. These findings provide causal evidence demonstrating the critical role of the AIC in interoceptive attention. Traditionally, the insular cortex is considered as a limbic sensory region that participates in the intuitive processing of complex situations (Augustine, 1996;Butti and Hof, 2010;Menon and Uddin, 2010) by integrating the ascending visceromotor and somatosensory inputs with attention systems via intrinsic connectivity to identify and respond to salient stimuli (Menon and Uddin, 2010;Uddin, 2015). The AIC, in particular, is a node that mediates cognitive processes including bottomup control of attention (Corbetta et al., 2002;Corbetta et al., 2008;Wu et al., 2015) and conscious detection of signals arising from the autonomic nervous system (Craig, 2002;Critchley, 2004). Therefore, the behavioral deficit of interoceptive attention in patients with AIC lesions is due to the disruption in the integration of the somatic and visceral inputs with the abstract representation of the present internal state (i.e., the saliency of a certain type of signals). Consequently, it leads to failure in discriminating whether the displayed respiratory curve is different from internal states.
Most previous lesion studies indicated interoceptive deficits with AIC lesions (Critchley and Garfinkel, 2017;García-Cordero et al., 2016;Ibañez et al., 2010;Ronchi et al., 2015;Starr et al., 2009;Terasawa et al., 2015;Wang et al., 2014), supporting the conclusion that interoceptive accuracy relies on a widely distributed network with the insular cortex as a key node (Craig, 2002;Critchley and Harrison, 2013). However, the preservations of interoceptive processing (Khalsa et al., 2009) and self-awareness across a large battery of tests (Philippi et al., 2012) were documented in one patient with bilateral insular damages. These studies are mostly based on subjective report focusing on 'feeling/awareness' (Khalsa et al., 2009) that might be compensated by other brain structures such as the brainstem and subcortical structures, for example, nucleus tractus solitaries, the parabrachial nucleus, area postrema and hypothalamus , frontal and temporal regions, for example, amygdala, superior temporal gyrus, and temporal pole (García-Cordero et al., 2016;Shany-Ur et al., 2014). In the current study, the BDT challenged interoceptive attention that requires the integration of interoceptive awareness and accuracy. Our examination of interoceptive attention in patients with focal AIC lesions showed that lesions of the AIC were associated with a deficit in performance, indicating that the AIC is critical in supporting the precision of interoceptive processing.

Mechanisms of the AIC in relation to interoceptive attention
Interoceptive attention is the mechanism that coordinates the processing of bodily signals and higher-level representation of that information. The AIC reportedly encodes and represents bodily information (e.g., visceral states) and transmits this information to other neural systems for advanced computations in conscious perception and decision-making (Bechara and Naqvi, 2004;Flynn, 1999;   Gu and FitzGerald, 2014). The AIC is a key node of the large-scale network that detects information from multiple sources including objective visceral signals, generates subjective awareness (Craig, 2009;Gu et al., 2012;Gu et al., 2015;Kleckner et al., 2017;Seeley et al., 2007), and responds to the switch between networks that supports internal oriented processing and cognitive control (Menon, 2011;Menon and Uddin, 2010;Sridharan et al., 2008). Supporting this argument, we showed that the AIC is intrinsically connected to the somatosensory area of the PoCG and that this connection is positively modulated by interoceptive attention (relative to exteroceptive attention).
Other higher-level areas, for example, the ACC/SMA, FEF, and IFG of the so-called cognitive control network (CCN) (Fan, 2014;Wu et al., 2018), are also involved in the interoceptive process. This is supported by the results of the enhanced functional connectivity between the AIC and these regions. Both somatosensory afferents and a network that includes the AIC and the ACC are possible pathways of interoceptive attention (Khalsa et al., 2009). The AIC may play a central role in integrating sensory signals from the PoCG and visual cortex and sends top-down signals that guide sensation and perception through a dynamic interaction with sensory or bottom-up information. Somatosensory information concerning the internal state of the body is conveyed through the PoCG, as well as the visual signals in V2/3 containing the majority of external information. The topdown modulation of the AIC in interoceptive attention is accomplished by augmenting the efferent signals to the somatosensory cortices. This result is consistent with the argument that a first-order mapping of internal feeling is supported by insular and somatosensory cortices (Damasio, 2003) and that somatosensory information critically contributes to interoceptive attention (Khalsa et al., 2009).
In the BDT, interoceptive attention reflects a combination of the attention to the internal bodily signal (i.e., the breath) and the external visual stimulus (i.e., the curve). To coordinate perceptual processing, the AIC may distribute and balance the processes of external and internal information. The winning model and parameter inference from DCM provide evidence that interoceptive attention is achieved mainly by modulating the connectivity between the AIC and the somatosensory areas (PoCG), while exteroceptive attention is primarily modulated via the connectivity between the AIC and V2/3. We propose that the dynamic adjustment of the connectivity of the AIC to sensory cortices is the foundation of interoceptive attention for bodily signals, which is critical for homeostatic regulation, and of exteroceptive attention for external objects or inputs.

Interoceptive task in the respiratory domain
Although the neural correlates of interoceptive awareness have been studied by other tasks, such as the heartbeat detection task (Bechara and Naqvi, 2004;Critchley et al., 2004;Khalsa et al., 2009;Ring et al., 2015), the error rate arising from the difficulty in heartbeat counting or non-sensory process confounds are inherent to these cardioception designs (Kleckner et al., 2015;Ring et al., 2015). In contrast to cardioception, breath can be clearly perceived and autonomously controlled. This feature enabled us to design a task measuring interoceptive attention, which requires that the target of interoception be clearly and vividly perceivable by our consciousness. The positive correlations between objective interoceptive accuracy during the BDT and subjectively scored difficulty of interoceptive task relative to exteroceptive task (i.e., interoceptive awareness) further demonstrates that the BDT is valid in assessing interoceptive attention. We developed the respiratory interoception task as a non-intrusive measurement with low cognitive load that is more practical for patients with focal brain damage than the demanding cardiac interoception tasks.
The BDT may not represent a pure probe of interoception because respiratory processes can also be tracked using exteroceptive and proprioceptive information. Thus, the participants possibly relied on a mix of interoceptive, exteroceptive, and proprioceptive information to perform the task. In our design, we included the DDT for a measure of exteroception so that the cognitive subtraction of DDT from BDT leaves the interoceptive and proprioceptive processing components of interoception (Gu et al., 2013). In the BDT, the delayed manipulation in our study was fixed to 400 ms, approximately 1/10 of an average cycle of normal healthy people (i.e., 3-4 s/cycle). This delayed duration can be manipulated according to each individual's respiratory cycle in an effort to control subjective task difficulty across participants.

Interoceptive attention
Depending on the source of information, attention can be categorized into (1) interoceptive attention, which is directed toward bodily signals such as somatic and visceral signals (e.g., in a heartbeat detection or counting task); (2) exteroceptive attention, which is directed toward primary sensory inputs from outside (e.g., visual and auditory stimuli); and (3) executive control of attention, which coordinates thoughts and actions (e.g., in Color Stroop, flanker, and working memory tasks; see review Fan, 2014). Although extensive studies have focused on the attentional modulation of sensory and perceptual inputs and on the executive control of attention, interoceptive attention is difficult to study because the vast majority of intrinsic visceral activity, except breath effort, cannot be clearly perceived under normal conditions. Using the BDT to examine attentional deployment toward breath effort enabled us to reveal the neural mechanism of interoceptive attention. In general, the perceptible, controllable, measurable, and autonomous features of breathing guaranteed more accurate and reliable measurement of individual differences in interoceptive ability.
As a type of interoceptive attention that could be clearly perceived and autonomously controlled, breath plays a potentially important role in generating and regulating emotion. For example, mindfulness meditation, which is now well known for its role in emotion regulation and mental health (Khoury et al., 2015), can be viewed as a practice involving interoceptive attention. One of its primary methods is to bring one's attention (the processing) and then awareness (the outcome) to the current experiences of the movement of the abdomen when breathing in and out or the breath as it enters and exits the nostrils. Considering the revealed neural mechanisms of interoceptive attention in this study, we predict that the AIC plays an important role in meditation. Findings that meditation experience is associated with increased gray matter thickness in the AIC (Lazar et al., 2005) and increased gyrification (increase in folding) of the AIC (Luders et al., 2012) support this prediction. Meditation training may enhance interoceptive attention to focus on bodily signals so that accurate feelings can be generated based on the bodily responses, meanwhile the mind can be released from an intensive involvement of exteroceptive and executive control of attention (the internal attention for the coordination of thought processing) that consumes the majority of mental resources.

Conclusion
This study provided important evidence of the involvement of the AIC in interoceptive attention by the fMRI studies and further demonstrated that the AIC is critical for the process by the lesion study. The converging evidence also suggests that interoceptive attention is achieved through top-down modulation from the AIC to the somatosensory and sensory cortices. In addition, the implementation of the interoceptive task extends the research on interoceptive processing into the respiratory domain with the validity and reliability demonstrated. It may have significant applications in studying issues related to interoceptive attention in patients with neuropsychiatric disorders, such as anxiety (Avery et al., 2014) and autism Quattrocki and Friston, 2014), and in patients with substance use disorders (Sö nmez et al., 2017).

Task design Task implementations
A respiratory transducer (TSD201, MRI compatible, BIOPAC Systems Inc), which was fastened around the participants' upper chest, was utilized to record breathing effort by measuring thoracic changes in circumference during respiration. The signal for the change in circumference was sampled at 1000 Hz using the BIOPAC MP150/RSP100C system, passing through a DC amplifier with lowpass filtering at 1 Hz and high-pass filtering at 0.05 Hz, and gain set to 10 V. Analog signal was then digitized by an A/D converter (USB-1208HS-4AO, Measurement Computing, Inc) and sent to a USB port of the test computer (Figure 1a). The task program in E-Prime (Psychology Software Tools, Pittsburgh, PA, USA) served as an interface through which the digitized signal from the USB port was received and presented to the participants on the computer screen as a continuous blue breath curve extending from left to right as time elapsed (Figure 1b), which was representative of their breathing effort. The breath curve was presented either with or without a delay (Figure 1c).
For the engagement of interoceptive attention during BDT, the participants were required to judge whether the presented breath curve was delayed relative to the breath rhythm they perceived from their body. In half of the trials, the displayed breath curve was synchronized with the participant's own respiration. In the other half, the displayed breath curve was delivered after a 400 ms delay period relative to the participant's own respiration (i.e., the plotting of the point on the extending curve was actually the point saved 400 ms before the current time point). Note that the parameter of 400 ms delay was determined based on a proportion (~1/10) of an average respiratory cycle of normal healthy people (i.e., 3-4 s/cycle). For the engagement of exteroceptive attention, the DDT was performed. The participants were instructed to detect whether a red dot flashed on the respiratory curve at any time when the breath curve was displayed. In half of the trials, a red dot flashed (30 ms for the fMRI experiment, and 50 ms for the lesion study) at a randomized time point on the breath curve. Figure 1b illustrates the two tasks. The stimuli for these two tasks were the same, consisting of four trial types reflecting the combination of presence or absence of a delay and presence or absence of a dot (Figure 1c). The factor of attentional deployment involved directing attention interoceptively to 'respiration' or exteroceptively to 'dot' in the BDT and DDT, respectively. During the two tasks, the participants were instructed to breathe as usual without holding or forcing their breath.
The participants were asked to perform each task in a blocked fashion in the interoceptive and exteroceptive runs. The fMRI experiment consisted of two runs, with one run for the BDT and the other run for the DDT. Each run, which included 60 trials, began and ended with a 30 s blank display, and each trial lasted 18 s, with an average inter-trial interval of 2 s, for a total of 21 min per run. Each trial began with a 3 s "Relax" display, followed by a 12 s respiratory curve presented with or without a 400 ms delay, and ended with a 3 s response window during which the participants made a forced-choice button-press response, prompted by the presentation of two alternative choices to indicate their response (Figure 1b). After the fMRI scan, the participants were asked to indicate the subjective difficulty they felt for each task on a 1-10 scale, with higher value indicating higher difficulty. For the lesion study, the same tasks were employed, one run for each task, with 40 trials in each run.

Behavioral data analysis
Interoceptive attention is associated with the objective accuracy in detecting bodily signals, the subjective belief in one's ability to detect bodily signals in general (i.e., sensibility), and the correspondence between objective accuracy and subject report (i.e., metacognitive awareness about one's performance when detecting bodily signals) (Garfinkel and Critchley, 2013;Garfinkel et al., 2015). Objective interoceptive/exteroceptive accuracy was calculated as the overall correct response rate during the BDT/DDT. In addition, we used signal detection theory to index detection sensitivity and response bias. Signal detection theory characterizes how perceivers separate signal from noise, assuming that the perceiver has a distribution of internal responses for both signal and noise (Snodgrass and Corwin, 1988;Stanislaw and Todorov, 1999). A fundamental advantage of signal detection theory is the distinction between sensitivity (ability to discriminate alternatives) and bias (propensity to categorize 'signal' or 'noise'). For the BDT, the sensitivity index (d') was calculated as d'=Z hit rate -Z false alarm rate , where the hit rate is the proportion of trials with delayed breath curve and responded as 'yes', and the false alarm rate is the proportion of trials with non-delayed breath curve and responded as 'yes'. A higher value of d' indicates a better interoceptive accuracy, whereas a value of 0 represents that the performance is at the chance level. The response bias index b, which represents the position of the subjective decision criterion, was defined as b = exp (d'ÂC), where C = -(Z hit rate +Z false alarm rate )/2. Index b corresponds to the distance of participants' estimated criterion to ideal observer criterion, and a value of 1 indicates no bias. For the DDT, indices of d' and b were calculated using the same formula, with the dot present as 'signal' and dot absent as 'noise'. Relative interoceptive accuracy was defined as the difference in performance accuracy between the BDT and DDT to control for non-specific performance effects .
An individual's subjective account of how they experience internal sensation and perception represents an alternative aspect of interoceptive processing, namely sensibility (Garfinkel et al., 2015). In the fMRI experiment, the subjective sensibility of interoceptive processing was measured using the self-report questionnaire of Body Perception Questionnaire (BPQ) (Porges, 1993). The subjective perception of one's performance during the BDT represents the awareness aspect of interoceptive attention (Garfinkel et al., 2015), which was measured via the subjectively scored difficulty of the BDT relative to the DDT. Note that the confidence ratings of the performance would be more directly related to awareness of interoception, but the measures of subjective difficulty and confidence ratings should be closely related. The correlation between the relative interoceptive accuracy and these indices of subjective sensibility and awareness was calculated to examine the relationship between the perceived (subjective) and actual measured (objective) performance of interoceptive attention. In addition, we conducted BFs of these correlation coefficients using JASP (The JASP Team, 2018). A BF larger than three suggests a significant correlation, whereas a BF smaller than 1/3 indicates a null correlation.

fMRI experiments Participants
The fMRI experiments included two samples of participants: the first sample included 44 adults (23 females and 21 males, mean age ± standard deviation: 21.43 ± 2.51 years, age range: 19-29 years), and the second sample included additional 28 adults (14 females and 14 males, mean age ± standard deviation: 21.93 ± 2.11 years, age range: 18-26 years). All participants underwent the same experimental procedures, except that pulse and respiratory signals were recorded for the second sample using the pulse sensor (Siemens Peripheral Pulse Unit, PPU_098) of the scanner and BIOPAC, respectively. All participants were right-handed (except for one participant), reported normal or corrected-to-normal vision, and had no known neurological or visual disorders. All participants completed questionnaires indexing subjective interoceptive sensibility (BPQ), symptoms of anxiety (Hamilton anxiety scale, HAMA;Hamilton et al., 1976), depression (Beck Depression Inventory, BDI; Knight, 1984), and positive and negative affective experience (PANAS) (Watson, 1988). They provided written informed consent in accordance with the procedures and protocols approved by The Human Subjects Review Committee of Peking University.

fMRI data acquisition and preprocessing
During functional scanning, the participants performed the BDT and DDT in separate runs that required them to attend to either their respiration or a visual flash dot, respectively. All neuroimaging data were acquired on a MAGNETOM Prisma 3T MR scanner (Siemens, Erlangen, Germany) with a 64-channel phase-array head-neck coil. During the tasks, blood oxygen level-dependent (BOLD) signals were acquired with a prototype simultaneous multi-slices echo-planar imaging (EPI) sequence (echo time, 30 ms; repetition time, 2000 ms; field of view, 224 mm Â224 mm; matrix, 112 Â 112; inplane resolution, 2 mm Â2 mm; flip angle, 90 degree; slice thickness, 2.1 mm; gap, 10%; number of slices, 64; slice orientation, transversal; bandwidth, 2126 Hz/Pixel; slice acceleration factor, 2). For the second cohort, the thickness was changed to 2 mm with a gap of 15%, and the number of slices was changed to 62. Field map images were acquired using a vendor-provided Siemens gradient echo sequence (gre field mapping: echo time 1, 4.92 ms; echo time 2, 7.38 ms; repetition time, 635 ms; flip angle, 60 degree; bandwidth, 565 Hz/Pixel) with the same geometry and orientation as the EPI image. A high-resolution 3D T 1 structural image (3D magnetization-prepared rapid acquisition gradient echo; 0.5 mm Â0.5 mm Â 1 mm resolution) was also acquired. Image preprocessing was performed using Statistical Parametric Mapping package (SPM12, RRID: SCR_007037; Welcome Department of Imaging Neuroscience, London, United Kingdom). EPI volumes were realigned to the first volume, corrected for geometric distortions using the field map, coregistered to the T 1 image, normalized to a standard template (Montreal Neurological Institute, MNI), resampled to 2 Â 2 Â 2 mm 3 voxel size, and spatially smoothed with an isotropic 8 mm full-width at half-maximum Gaussian kernel.
fMRI: analysis of the first sample Image statistical parametric mapping Imaging data from the two samples were analyzed separately and independently, with the exploratory whole brain analysis conducted with the first sample and the confirmatory ROI analysis conducted with the second sample. For the whole brain analysis of the first sample, statistical inference was based on a random-effects approach (Penny and Holmes, 2007), which comprised two steps: first-level analyses estimating contrasts of interest for each subject followed by second-level analyses for statistical inference at the group level. For each participant, first-level statistical parametric maps of BOLD signals were modeled using general linear modeling (GLM) with regressors defined for each run with the four trial types: 2 breath curve delay (non-delayed, delayed)Â2 dot present (no dot, dot). Each trial was modeled as an epoch-related function by specifying an onset time and a duration of 12 s. The corresponding four regressors were generated by convolving the onset of each trial with the standard canonical hemodynamic response functions (HRF) with a duration of 12 s, that is, by convolving the trial block with HRF, equivalent to a box-car function. Six parameters generated during motion correction were entered as covariates of no interest. The time series for each voxel were high-pass filtered (1/128 Hz cutoff) to remove low-frequency noise and signal drift.
Contrast maps for interoceptive vs. exteroceptive attention (BDT -DDT), the presence of breath curve delay (delayed -non-delayed), and the interaction between them ([delayed -non-delayed] BDT -[delayed -non-delayed] DDT ) for each participant were entered into a second-level group analysis conducted with a random-effects model that accounts for inter-subject variability and permits population-based inferences. The statistical maps were corrected for multiple comparisons using Gaussian random field (GRF) theory (T > 3.29, cluster-wise p<0.05, GRF corrected) with a minimum cluster size of 420 resampled voxels. Note that changes in neural activity revealed by the main effect of interoceptive vs. exteroceptive attention (the contrast of BDT vs. DDT) could also reflect task-specific effects, such as differences in task difficulty or respiratory characteristics (i.e., amplitude and frequency) between the two tasks, in addition to effects of change in attentional focus. Although the main effect of interoceptive vs. exteroceptive attention (the contrast of BDT vs. DDT) is subject to confounding by the task-specific effects, the interaction effect can disentangle those effects (i.e., cancel out the breathing effort difference between the two tasks). This interaction reflected the brain response when directing attention to the feedback mismatch during interoceptive processing while controlling for the non-specific effect (i.e., the physical difference in feedback stimulus between delayed and non-delayed curves during exteroceptive processing). Therefore, a positive interaction effect represents brain response to the interoceptive processing above and beyond the physical feedback difference.

Correlation between interoceptive accuracy and the interaction effect of the AIC
To test for a linear correlation between AIC activation and behavioral performance on the BDT, we entered each participant's interaction contrast maps into the second-level random-effects group regression analysis, together with their individual accuracy in the BDT as the variable of interest and accuracy in the DDT as the covariate. Threshold of significance was GRF-corrected at p<0.05 (T > 2.42) with a cluster extent of 106 contiguous voxels (resampled), corrected using small-volume ROI correction. The mask image was generated from an anatomical template of the bilateral insular cortex based on the Automated Anatomical Labeling template (Tzourio-Mazoyer et al., 2002).

PPI analysis
PPI analysis provides a measure of change in functional connectivity between different brain regions under a specific psychological context (Friston et al., 1997). We conducted PPI analyses using a moderator derived from the product of the activity of a seed region (i.e., the AIC) and the psychological context (i.e., interoceptive in contrast to exteroceptive attention, BDT vs. DDT). The ROI selection was independent of the interoceptive attention process that was used as the psychological context: The left and right AICs were first identified from the main effect of the breath curve delay (the contrast of delayed versus non-delayed) in the GLM. We then conducted two whole-brain PPI tests for the right and left AIC, reflecting changes in functional connectivity between the seed region time series (physiological regressor) and other brain regions as a function of interoceptive relative to exteroceptive attention (BDT vs. DDT, psychological regressor). The AIC time series of each participant were extracted from a 6 mm-radius sphere centered at the peak of the AIC (right AIC: x = 30, y = 26, z = À4; left AIC: x = À30, y = 24, z = À4). The PPI term was calculated as an element-by-element product of the deconvolved physiological regressor and psychological regressor, which was then reconvolved with the canonical HRF. The generated PPI model included the PPI term, the physiological regressor, the psychological regressor, and nuisance regressors of six motion parameters. The threshold of significance for the second-level group data analysis of the images from the PPI regressor was determined the same as in the GLM. Regions identified as significant clusters have two possible interpretations: (1) the connectivity between the AIC and these regions was altered by the psychological context, or (2) the response of these regions to the psychological context was modulated by AIC activity. To simplify the explanation, we used the first interpretation throughout this article.

DCM analysis
DCM (Friston et al., 2003) is used to disambiguate different potential network structures by inferring hidden neuronal states from measurements of brain activity. DCM distinguishes between endogenous coupling and context-specific coupling, which could account for the effects of experimentally controlled network perturbations. Considering the inherent limited causal interpretability of the PPI analysis for the direction of interaction, we only conducted DCM to explain the potential mechanisms of interplay between AIC and other brain areas involved in interoceptive attention. The ROI of the right AIC in the DCM was the same as that in the PPI analysis. The other regions included in the DCM were selected based on significant positive and negative PPI results and with the coordinates of the ROIs identified by the group level T-contrast of all conditions versus baseline. Data from one participants were excluded from the DCM analysis because activity in one of the ROIs could not be identified.
A three-area DCM was specified for all participants with bidirectional endogenous connection between the right AIC and the other two ROIs, and with the main effect of 'all stimuli' as the driving input entering the other two ROIs. Five base models were generated by specifying possible modulations of interoceptive and exteroceptive attention (BDT and DDT, respectively) on the four endogenous connections between ROIs. These base models were then systematically elaborated to produce 52 variant models, which included all possible combinations of the modulation of interoceptive and exteroceptive attention (BDT and DDT, respectively) on endogenous connections between the right AIC and the two other ROIs.
Model comparison was implemented using random-effects BMS in DCM12 to determine the most likely model of the 52 models given the observed data from all participants (Stephan et al., 2009). The RFX analysis computes exceedance and posterior probabilities at the group level, and the exceedance probability of a given model denotes the probability that this model is more likely than all other models considered (Stephan et al., 2009). To summarize the strength of effective connectivity and its modulation quantitatively, we used random-effects BMA to obtain average connectivity estimates (weighted by their posterior model) across all models and all participants . We conducted one-sample t tests on the subject-specific BMA parameter estimates to assess their consistency across subjects with Bonferroni correction for multiple comparisons.

fMRI: ROI analyses of the second sample
Whereas whole brain analyses of the first sample aimed at identifying brain areas involved in interoceptive processing, ROI analyses of the second sample aimed to confirm that the effects found from the first sample were not confounded with the effects induced by other physiological signals. Change in BOLD signals can be due to direct neural activity (induced by experimental manipulation) or an indirect effect (such as vascular response, which would be considered as a confounding effect). For example, the cerebral vascular response is sensitive to the circulation of carbon dioxide (CO 2 ) and oxygen (O 2 ), and causes a change in global cerebral blood flow (CBF) and global BOLD signal. It is evident both in human and animals that the global CBF and global BOLD responses influence local stimulus-induced hemodynamic response to neural activation (Cohen et al., 2002;Friston et al., 1990;Ramsay et al., 1993). In general, a larger local stimulus-induced BOLD response occurs when global BOLD is lowered, whereas a smaller local stimulus-induced BOLD response occurs when global BOLD is elevated. In our study, the experimental manipulation of interoception was likely to cause a change in respiratory characteristics (i.e., circulation of CO 2 and O 2 ). The difference in physiologic states between the BDT and the DDT might cause a change in global BOLD signals. Thus, the effect resulting from local interoception-related BOLD responses would be confounded by the global hemodynamic influence.
To partial out the potential confounding, we processed physiological data, including cardiac pulsation and respiratory volume collected in the second sample, by using the PhLEM toolbox (http:// sites.google.com/site/phlemtoolbox/). Physiological noise correction consisted of (1) regressing out time-locked cardiac and respiratory effects, and their interaction effect using a modification of the conventional RETROICOR approach (Brooks et al., 2008;Glover et al., 2000), and (2) regressing out low-frequency respiratory and heart rate effects using the RVHRCOR approach (Verstynen and Deshpande, 2011). In RETROICOR, a cardiac phase calculated from a pulse oximeter was assigned to each acquired image in a time series (Hu et al., 1995), and a respiratory phase was assigned to a corresponding image using the histogram equalized transfer function that considers both the respiratory timing and depth of breathing (Glover et al., 2000). The conventional RETROICOR approach (Glover et al., 2000) defines low-order Fourier terms (i.e., sine and cosine values of the principal frequency and the 2 nd harmonic) to model the independent effects of the cardiac and respiratory fluctuation, which is considered insufficient to remove variations caused by physiological artifacts (Harley and Bielajew, 1992;Tijssen et al., 2014). Therefore, we used additional terms of higherorder Fourier expansions (i.e., to the 5th harmonics) in RETROICOR, and formed multiplicative sine/ cosine terms that consider the interaction between cardiac and respiratory effects. In specific, the interaction terms were calculated by Sinð' c AE ' r Þ and Cosð' c AE ' r Þ, where ' c; r is the cardiac or respiratory phase, consisting of a mixture of third-order cardiac and second-order respiratory harmonics. In RVHRCOR, two nuisance regressors were generated by convolving respiratory variations (RVs) and heart rate (HR) with 'respiration response function (RRF)' and the 'cardiac response function (CRF)' respectively. In specific, RV was computed as the root-mean-square amplitude of the respiration waveform across a 6 s sliding window, and HR was computed as the inverse of the average beat-to-beat interval in a 6 s sliding window (Chang et al., 2009). Therefore, the physiological correction contained a total of 46 regressors, of which 20 were from independent time-locked cardiac and respiratory effects, 24 were from interaction terms, and two were from low-frequency RV and HR effects. Statistical parametric maps were generated using the same GLM as in the whole brain analyses, with motion parameters and these physiological regressors entered as covariates of no interest.
To avoid double dipping, we defined the ROIs based on the first sample. In specific, the ROIs (i. e., left and right AICs) were the clusters of the second-level group analysis results of the interaction effect ([delayed -non-delayed] BDT -[delayed -non-delayed] DDT ). Parameter estimates were extracted from each ROI in the second sample under the four experimental conditions of each participant, and then entered into a two-way repeated-measures ANOVA. We also examined correlations between the interaction effect of each ROI and behavior measures (i.e., relative interoceptive accuracy) across participants.
To further examine the degree to which physiological correction impacted the whole brain activation, we conducted whole brain paired t-tests by comparing the contrast maps without and with physiological correction at an extremely permissive threshold (voxelwise p<0.05 uncorrected).

Lesion study Brain lesion patient and control groups
Six male patients (33-53 years old, mean 42.17 ± SD 7.31 years) with focal unilateral insular cortex lesions participated in the lesion study (see Table 9 for patient characteristics). Two patients had a right-side lesion, and four patients had a left-side lesion. In addition, six patients with focal lesions in regions other than the insular cortex (i.e., temporal pole, n = 3, lateral frontal cortex, n = 2, and superior temporal gyrus, n = 1) were recruited as BDCs, and 12 neurologically intact participants were recruited as NCs. All lesions were resulted from the surgical removal of low-grade gliomas. All patients were recruited from the Patient's Registry of Tiantan Hospital, Beijing, China. NC participants were recruited in the local community. All NC participants were right-handed, had normal color vision, and reported no previous or current neurological or psychiatric disorders. BDC patients matched with patients with insular cortex lesions in chronicity (t(10) = À0.36, p=0.38, BF = 0.48), and neither group significantly differed from the NC group in age (AIC vs. NC: t(9) = À1.01, p=0.18, BF = 0.61; BDC vs. NC: t(10) = À1.80, p=0.06, BF = 1.25) nor education (AIC vs. NC: t(10) = 0.77, p=0.25, BF = 0.52; BDC vs. NC: t(6) = À1.04, p=0.18, BF = 0.72). All six insular lesion patients were considered cognitively intact, as determined by Mini-Mental State Examination (MMSE), a measurement of cognitive impairment (Folstein et al., 1975), and the raw scores of MMSE were not significantly different from either BDCs (t(10) = À1.30, p=0.13, BF = 0.78) or NCs (t(10) = À1.76, p=0.06, BF = 1.17). Compared with the NCs, the patients with insular lesions did not show significant alteration in baseline mood indexed by the BDI score, compared to NCs (t(6) = 1.70, p=0.06, BF = 1.84) or BDCs (t(9) = 0.65, p=0.28, BF = 0.53). Demographic information of the groups can be found in Table 9. By chance, all the patients with AIC lesions were male. All participants were informed of the study requirements and provided written consent prior to participation. The patient study was approved by the Institutional Review Board of the Beijing Tiantan Hospital, Capital Medical University.

Lesion reconstruction
Two neurosurgeons, blinded to the experimental design and behavioral results, identified and mapped the lesions of each patient onto a template derived from a digital MRI volume of a normal control (ch2bet.nii) embedded in the MRIcro program (RRID: SCR_008264; http://www.cabiatl.com/ mricro/mricro/index.html). In each case, lesions evident on MRI were transcribed onto corresponding sections of the template to create a volume of interest image. This volume of interest image was then used to measure the location (in MNI coordinates) and volume (in mL) of individual lesions and to create within-group overlaps of lesions using the MRIcro program.

Behavioral data analysis of the lesion study
We used non-parametric analysis (Feys, 2016) to test the two-way interaction between group (AIC, BDC, and NC) and task (BDT and DDT) using R (R Development Core Team, 2013) because the small sample data sets did not meet the assumption of parametric tests. In specific, we used the npIntFactRep function (from the npIntFactRep package) that yielded an aligned rank test for interaction in the two-way mixed design with the group (AIC, BDC, and NC) as the between-subject factor and with the task (BDT and DDT) as the within-subject factor. If the interaction was significant, we used the non-parametric bootstrapping method to test the simple between group effects for each task separately. The bootstrapping procedure was conducted with 10,000 iterations (Hasson et al., 2003;Mooney and Duval, 1993). If the probability of obtaining the observed t-value was less than 5% (one-tailed), we considered the difference between the two groups to be significant. We used one-tailed tests because we hypothesized that lesions of a specific brain region (i.e., the AIC) would induce deficits in behavioral response. In addition, we calculated BFs with Cauchy prior to determine the relative strength of evidence for the null and alternative hypotheses (Dienes, 2014;Dienes and Mclatchie, 2018). The value of BF means that the data are BF times more likely under the alternative than under the null hypothesis. The standard value for assessing substantial evidence for the null is BF <1/3 and for the theory against null is BF >3, whereas values between 1/3 and 3 are counted as data insensitivity. The BFs were calculated using JASP (The JASP Team, 2018).