Threat-Related Corticocortical Connectivity Elicited by Rapid Auditory Looms

Threat-Related Corticocortical Connectivity Elicited by Rapid Auditory Looms | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results Threat-Related Corticocortical Connectivity Elicited by Rapid Auditory Looms Karolina Ignatiadis , Roberto Barumerli , Gustavo Deco , Brigitta Tóth , View ORCID Profile Robert Baumgartner doi: https://doi.org/10.1101/2024.06.21.600106 Karolina Ignatiadis 1 Acoustics Research Institute, Austrian Academy of Sciences , Wohllebengasse 12-14, Vienna, 1040, Austria Find this author on Google Scholar Find this author on PubMed Search for this author on this site For correspondence: karolina.ignatiadis{at}oeaw.ac.at robert.baumgartner{at}oeaw.ac.at Roberto Barumerli 1 Acoustics Research Institute, Austrian Academy of Sciences , Wohllebengasse 12-14, Vienna, 1040, Austria Find this author on Google Scholar Find this author on PubMed Search for this author on this site Gustavo Deco 2 Computational Neuroscience Research Group, Pompeu Fabra University , carrer de la Mercè, 12, Barcelona, 08002, Spain Find this author on Google Scholar Find this author on PubMed Search for this author on this site Brigitta Tóth 3 Institute of Cognitive Neuroscience and Psychology, HUN-REN Research Centre for Natural Sciences , Magyar tudósok körútja 2, Budapest, 1117, Hungary Find this author on Google Scholar Find this author on PubMed Search for this author on this site Robert Baumgartner 1 Acoustics Research Institute, Austrian Academy of Sciences , Wohllebengasse 12-14, Vienna, 1040, Austria Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Robert Baumgartner For correspondence: karolina.ignatiadis{at}oeaw.ac.at robert.baumgartner{at}oeaw.ac.at Abstract Full Text Info/History Metrics Data/Code Preview PDF Abstract While sounds of approaching objects are generally more salient than those of receding ones, the traditional association of this auditory looming bias with threat perception is subject to debate. Differences between looming and receding sounds may also be learned through non-threatening multisensory information, or influenced by confounding stimulus characteristics. To investigate, we analyzed corticocortical connectivity patterns from electroencephalography, examining the preferential processing of looming sounds under different attentional states. To simulate rapid distance changes we used complementary distance cues, previously studied in the looming bias literature. Notably, despite the absence of conscious threat perception, we observed crucial involvement of frontal cortical regions typically associated with threat and fear responses. Our findings suggest an underlying bias towards the ventral ‘what’ stream over the dorsal ‘where’ stream in auditory information processing, even when the participants’ task was solely focused on the discrimination of movement direction. These results support the idea, that the perceptual bias towards looming sounds reflects an auditory threat detection mechanism, while offering insights into the neural function involved in processing ecologically relevant environmental cues. 1 Introduction If a car is approaching from a distance, its timely detection and avoidance are essential to our survival. It is presumably to improve warning capacity, that approaching stimuli are more salient than receding ones. This perceptual asymmetry, often referred to as the auditory looming bias, has been found present across species [ 1 – 5 ] and ages [ 6 – 9 ], making it a rather universal trait. As an effective warning mechanism, the looming bias should have the capacity to readily capture attention and be rather universal across cue types. Corroborating this hypothesis, signatures of the bias have indeed been found across attentional states and auditory distance cue types; they were present already at the level of Heschl’s gyrus (HG), housing the primary auditory cortex [ 10 ]. Yet the bias’ relationship to threat detection remained a hypothesis, and discrepancies in behavioral performance, timing, and attentional amplification, suggest that there are differences in cortical processing depending on these factors. The notion of stimulus-specificity has frequently been put forth with regards to the selective advantage of the looming bias, and its function as a warning mechanism for organisms facing potential collisions with sound sources. It tends to be observed more consistently in response to stimuli with a natural overtone structure, in contrast to Gaussian white noise stimuli, which sound arguably more artificial [ 4 , 9 , 11 ]. However, studies have also demonstrated looming biases in response to noise stimuli when accounting for the natural acoustic filtering properties of listeners [ 12 , 13 ]. This suggests that the absence of natural spatial cues, rather than the source’s identity, may be responsible for the failure to elicit the bias in certain cases. Additionally, some investigations into the looming bias have employed auditory distance changes as short as 10 ms [ 10 , 12 , 13 ], prompting questions about the necessary identity and ecological validity required to evoke this effect. From a neuroimaging perspective, increased amygdala activation in response to slowly rising sound intensities has been an important argument for the bias’ warning function [ 14 ]. Apart from Heschl’s gyrus and the amygdala, functional magnetic resonance imaging (fMRI) has highlighted the involvement of the temporal plane, superior temporal sulcus (STS), prefrontal cortex (PFC), and inferior parietal lobe (IPL) in the preferential processing of looming sounds [ 14 , 15 ]. In general, auditory stimuli have been hypothesized to follow two parallel cortical processing streams: one following a ventral and the other a dorsal path [ 16 ]. Originally stemming from visual research, the dorsal pathway is associated with spatial perception (”where”), while the ventral stream with object identification (”what”) [ 16 , 17 ]. The inferior parietal lobe, as part of the dorsal auditory pathway, is thought to play a crucial role in spatial hearing [ 18 , 19 ] and sound motion processing in particular [ 20 ], whereas superior temporal sulcus and prefrontal cortex belong to the ventral pathway. Based on these findings, the looming bias circuit emerges as an extended distributed cortical network. Besides the mere activation increase induced by looming sounds, one crucial aspect is the way in which the involved regions are at interplay. This question can be addressed through functional connectivity investigations [ 21 ]; namely computational methods exploring the information exchange among regions of interest (ROIs). Unlike structural connectivity, which describes anatomical connections linking sets of neural elements, functional connectivity is dynamic in nature. It represents changes in statistical inter-dependencies between or among brain regions, within a specific time interval and connected to an event of interest. The observed brain regions that are found to contribute the most with regards to connectivity, relative to the others or the combinations thereof, are defined as functional hubs. Those are also dynamic and may deviate from an anatomical definition, as they can be a part of different functional clusters [ 22 ]. Findings stemming from available connectivity analyses of the auditory looming bias circuit are inconclusive: a study based on intensity ramps argues that top-down directional causal influence from prefrontal cortex to Heschl’s gyrus enhances processing of looming versus receding sounds [ 23 ], while prior investigations on spectral stimuli argue for a bottom-up, temporofrontal connectivity [ 13 ]. Functional connectivity methods employed in previous research focus on a bidirectional analysis process, albeit relying on a small preselection of brain regions (Granger [ 24 ], Phase Transfer Entropy [ 25 ]). Although insightful regarding the interplay of the considered ROI pairs, further methods may offer an approach that is closer to a network structure. They are nevertheless limited by either the number of regions that can be considered (conditional Granger Causality [ 26 ]), the dimensionality (multivariate Granger [ 27 ]) or the amount of constraints enforced by parameters of a model (e.g., Dynamic Causal Modeling [ 28 ]). Contrary to that, recent frameworks offer both a holistic as well as data-driven approach [ 29 , 30 ]. They provide the possibility to investigate the whole brain on different levels, without the necessity of a predefined set of ROIs or network structure parameters. Of those, the INSIDEOUT approach [ 29 ] relies on the observation, that the environment drives hierarchically lower, sensory regions, stronger than hierarchically higher ones. A system in equilibrium has seamless transitions between different states; thermodynamically, it is reversible in time. Should the system get driven out of equilibrium, the transitions between states become non-reversible and an arrow of time emerges. Measuring the effects of the extrinsic environment on the intrinsic brain dynamics through the non-reversibility of the system, here the brain, can therefore help uncover variations in brain states under different conditions. The framework of normalized directed transfer entropy (NDTE) [ 30 ], contrarily, works on a mesoscopic level: By considering the interconnectivity of all defined brain regions, it draws assumptions about the most essential contributors, or functional hubs, of the underlying networks. In the current study, we investigate the cortical connectivity network underlying the auditory looming bias under the individual factors of cue type and attention, in search for overlapping patterns along spatial and/or identity-related cortical processing streams. Through the high temporal resolution of electroencephalography (EEG) in combination with recently proposed, data-driven approaches for connectivity analyses, we investigate the brain at different levels of granularity [ 29 , 30 ]: First as a whole, and subsequently in search for the functional hubs that act as essential contributors in the looming network. High spatial resolution is achieved by complementing source localisation of high-density EEG with individual brain anatomies and electrode locations [ 31 ]. The present analyses are based on prior collected data, studying the auditory looming bias at the level of HG under the aforementioned factors of attention and cue type [ 10 ]. In that paradigm ( Fig. 1 ), participants listened to broadband harmonic tones that rapidly changed in their simulated distance from the listener and thereby elicited a looming or receding percept. Distance cues comprised either overall sound intensity or spectral shape changes. Listeners were first passively exposed to the stimuli while watching a silent subtitled movie and later had to discriminate the sonic motion direction. We find that there is to be higher sensitivity for intensity stimuli, while different main hubs, traditionally connected to threat and fear perception, emerge based on the factors considered. Download figure Open in new tab Fig. 1: Experimental design. A ) Looming and receding percepts created through simulated transition between two sounds of different intensities (top, blue) or spectral shapes (bottom, red). Thick arrows represent 50% transition probability for motion trials (dark = looming; light = receding). B ) Sound intensity over time (top panel) and magnitude spectrum (bottom panel) of all implemented stimuli. Figure adapted from [ 10 ]. 2 Results For the present connectivity investigations we extracted the source-localized EEG time series of all cortical regions, as defined by the Desikan-Killiany parcellation [ 32 ]. We considered the time interval between 0 and 300 ms relative to the event of distance change. This choice was made based on the finding, that this time window has shown significant biases evoked in HG in previous investigations [ 10 ]. 2.1 Intensity looms induce stronger non-reversibility in cortical processing INSIDEOUT reflects how the environment (extrinsic, outside) affects the dynamics and equilibrium of the underlying brain state (intrinsic, inside) [ 29 ], by measuring the non-reversibility of a considered system. We implemented this framework by accounting for the set of all ROIs of the considered parcellation, hence the cortex as a whole ( Sec. 4.4.1 ). Higher non-reversibility is thus understood as a quantification of the amount of change in causal interactions of the brain under each considered condition. Figure 2 shows the distributions of non-reversibility measures obtained for every considered condition. An ANOVA with the factors of attention, cue type and motion direction revealed the latter to be a significant factor ( F 1,27 = 8.34, η 2 = 0.24, p = 0.008), as well as its interaction with the cue type ( F 1,27 = 5.11, η 2 = 0.16, p = 0.032). To further investigate this interaction, we performed a separate ANOVA for each cue type and adjusted the p-values for multiple comparisons using Bonferroni correction. For intensity, looming sounds were found to elicit higher non-reversibility than receding sounds ( F 1,27 = 9.17, η 2 = 0.25, p = 0.011). No significant factors or interactions thereof appeared for the spectral condition. Hence, looming stimuli appeared to disrupt the intrinsic equilibrium more than receding ones, in particular when they are based on intensity changes. Download figure Open in new tab Fig. 2: Effects of cue type and movement direction on the temporal non-reversibility of cortical processing, as calculated by the INSIDEOUT framework. Points and bars represent the means and their standard errors within each violin plot. Asterisks indicate statistical significance (p < 0.05). 2.2 Connectivity hubs relevant to the auditory looming bias To better understand the dependencies, we then applied the NDTE framework, as it offers a more granular view on the interacting brain regions. Following the procedure suggested in [ 30 ], we considered each ROI’s connection to every other ROI in the cortical parcellation. The connectivity between each pair of regions was calculated on the actual data, and its significance assessed through a distribution of surrogate data stemming from the same ROI-pair. Aggregation of the connectivity information across subjects allowed for the construction of connectivity hubs, namely regions, or sets thereof, that are, as a whole, more connected compared to any other considered set comprising the same number of ROIs. We performed a connectivity analysis on the bias data by considering the factors of attention and cue type. The two quantities of essence in this framework are termed inflow ( G in in [ 30 ]) and outflow ( G out in [ 30 ]); they respectively represent the connectivity incoming to or outflowing from a ROI. If a set of ROIs is considered as a network, inflow is the sum of all incoming connectivity across all its constituent ROIs. The respective holds for the outflow. We determined the major inflow and outflow hubs per condition by following the concept and search procedure of functional rich clubs ( Sec. 4.4.2 ; FRICs in [ 30 ]). Following the procedure for their definition based on the inflow, we respectively defined the major hubs based on the connectivity outflow ( Fig. 3 , Sec. 4.4.2 ). As demonstrated in figure 3A , the major inflow was attributed to one region, except for the active intensity condition. The ROIs receiving the most inflow spanned over temporal regions (STG), frontal regions (pars opercularis - IFGoper, frontal pole - Fpole), and both hemispheres. Across both active conditions, only the frontal pole emerged as a crucial inflow hub for looming bias. In the passive conditions, the relevant inflow hubs comprised the right precentral gyrus (PreCG) for intensity and the left pars triangularis (IFGtriang) for spectral stimuli. Regarding the outflow hubs, one region emerged per condition and all regions were located in the left hemisphere. Apart from the active intensity condition, where the insular gyrus (IG) was identified as the main hub, temporal regions were identified for the remaining cases: superior temporal gyrus (STG) for active spectral, transverse temporal gyrus (Heschl’s Gyrus, HG) for passive intensity, and the banks of the superior temporal sulcus (BanksSTS) for passive spectral. Download figure Open in new tab Fig. 3: Major inflow and outflow hubs of looming bias identified per considered condition. Fpole (magenta) is activated by both active intensity and spectral. Fpole - frontal pole, IFGtriang - pars triangularis, STG - superior temporal gyrus, PreCG - precentral gyrus, IFGoper - pars opercularis, IG - insular gyrus, BanksSTS - banks of the superior temporal sulcus, HG - Heschl’s gyrus. We further extracted the pattern of hub connections that emerged, separately for inflow and outflow connectivity in each considered condition ( Fig. 4 ). On a large scale, the inflow hubs, localized to the frontal cortices (Fpole, IFG, PreCG), dominantly received information from more distant regions of the sensory temporal regions such as the superior temporal gyrus, Heschl’s gyrus or inferior temporal cortices. In contrast, the outflow hubs, localized to the temporal regions (IG, STG, HG, STS), tended to send information to more local areas within the temporal cortex. All those connections occurred mainly within hemisphere, reflecting rather weak inter-hemispheric connections to and from the hubs. Download figure Open in new tab Fig. 4: Detailed connectivity for major inflow and outflow hubs identified per condition. Displayed connections correspond to the top 30% of connectivity strength per condition. Color bars were adjusted accordingly and normalized connectivity values were scaled up by a factor of 100. 3 Discussion Looming sounds exhibit remarkable salience, consistently eliciting a perceptual bias compared to equivalent receding sounds. This bias is commonly hypothesized to signify a mechanism for threat detection and hazard protection. In this study, we employed two data-driven, state-of-the-art functional connectivity approaches, to examine cortical responses to simulations of rapid changes in auditory distance. We examined the responses during both passive and active listening, and while utilizing either intensity or spectral cues. Analyzing macroscopic brain states, temporal non-reversibility tests revealed a more pronounced impact of looming sounds on the overall functional connectivity when intensity was employed to simulate motion in distance. Taking a more granular approach, we identified functional connectivity hubs for each condition, shedding light on the intricate neural networks underlying these perceptual processes. Throughout our conditions, frontal regions emerged as the main inflow hubs (frontal pole, IFG including the pars opercularis, pars triangularis also known as Broca’s area-BA45), and temporal regions as the major outflow hubs (primary and secondary auditory cortices such as the superior temporal gyrus, transverse temporal gyrus including the PAC, banks STS). Prior studies have argued regarding the directionality of the connectivity in the looming bias, with results divided between bottom-up processing [ 13 ] or a top-down intervention [ 23 ]. The here-identified hubs support the bottom-up notion, as temporal regions act more as sources (outflow hubs) and frontal regions as receivers (inflow hubs) in response to looming sounds. Moreover, those hubs mostly appear to belong to the ventral auditory pathway. As seen in the case of passive intensity, the motor cortex is additionally present, while no hubs appear in our findings that could be clearly associated with the dorsal auditory pathway. This result is intriguing, as it suggests that, in the considered context, priority is given to recognizing and identifying the auditory sound source (”what”) rather than its location (”where”); even though listeners were solving a spatial auditory task during the active listening part of their study participation. The to be expected high involvement of the dorsal pathway in spatial perception might have been canceled out by the bias calculation (difference between looming and receding). The importance of source identity on the elicitation of auditory looming bias has previously also been demonstrated by comparisons between different types of source stimuli: tones elicited stronger biases than noise [ 4 , 5 , 9 , 33 , 34 ]. The here identified connectivity hubs seem to reflect the crucial role of source identity in the manifestation of the auditory looming bias. 3.1 Emerging looming bias hubs relate to fear responses The analysis of sound source identity is governed by cortical hubs that are congruent to various literature findings in the context of threat and fear. Among the frontal regions, the anterior part of the frontal cortex (Fpole) emerged as the major inflow hub in both active conditions, which is in agreement with its recognized function for executive cognitive processing and action selection [ 35 , 36 ]. Other auditory studies also demonstrated its involvement in the context of threatening sentences and emotionally salient pictures [ 37 ]. Apart from the frontal pole, parts of the inferior frontal gyrus appear as essential in facilitating the bias: the right pars orbitalis for the active intensity and left pars triangularis for passive spectral. More generally, studies of visually presented threat-related words have reported activation of the left inferior frontal gyrus [ 37 , 38 ]. Inferior frontal gyrus has additionally been activated in the context of fear conditioning [ 39 ] and the downregulation of psychophysiological reactions to threat [ 40 ]. The information received by the major inflow hubs primarily appeared as coming from primary and secondary auditory cortical regions on the temporal lobe. Concordantly, the main outflow hubs we identified were localized in the temporal areas. Being an auditory task, areas such as the PAC may be especially involved; yet previous research has also implicated it in threatening sound paradigms [ 37 ]. Superior temporal ROIs, emerging in our considered active (attentive) conditions, have additionally been connected to attention linked to threatening voices [ 41 ]. In contrast to mainly temporal regions as outflow hubs, the insula appeared as a key outflow node in the intensity stimuli. As a ROI, it has been associated with fear and anxiety conditioning [ 37 , 39 , 42 ], while animal studies have demonstrated its implication in fear or extinction memory [ 43 ]. Intensity stimuli have behaviorally and neurally emerged as more salient than their spectral counterparts, in both adult and newborn listeners [ 10 ]. Implication of the insula only in the perception of intensity stimuli might act as a contributor of that manifestation. Finally, frontotemporal activations have generally been linked to the basolateral amygdala (BLA), an essential hub of the limbic system, in the context of automatic fear detection [ 44 ]. The amygdala itself has, in turn, been further implicated in the looming perception [ 14 ] in a warning role. As we conducted our study with the use of EEG, subcortical activations, consequently also the amygdala, are either inaccessible or unreliable; the direct verification of the BLA-frontotemporal link in the context of the bias can thus not be made through our findings. Yet the emerging frontal and temporal connectivity hubs may be a manifestation of the BLA-frontotemporal exchange, congruently to previous findings: Invasive studies on animals have specifically implicated the medial prefrontal cortex and BLA in the discrimination between harmful and safe stimuli, and highlighted that the corticocortical dialogue between sensory and prefrontal areas is essential for fear-discrimination processes [ 45 ]. Taken together, the functional relevance of the major hubs we identified along the ventral auditory pathway suggests that, regardless of cue type, looming sounds elicit the perceptual bias by rapidly recognizing the sound as a potential threat. 3.2 Methodology and limitations In the current study, we utilized direct (NDTE) and indirect (INSIDEOUT) connectivity metrics in order to obtain an image of the bias-related processes on the cortical surface. Depending on the method at hand, investigations can be done at different levels of granularity. INSIDEOUT captures the breaking of causal connections through non-reversibility and the arrow of time in order to measure brain connectivity. Compared to other approaches, it has the big advantage that no underlying constraints (e.g., ROIs or networks) or models (e.g., directionality or node assumptions) are necessary for its implementation. It can additionally give a coarse representation of the different brain states based on the whole cortex in a significantly less computationally complex and time-consuming manner than conventional approaches would demand. In terms of non-reversibility, the looming bias was found mainly for the intensity stimuli. Broadly considered in looming studies, intensity stimuli have generally appeared more salient than spectral ones; the latter seem to be more complex in their understanding and cognitively processed in a much more subtle manner [ 46 – 48 ]. As INSIDEOUT is reflective of subjective conscious awareness [ 29 ], our result corroborates the difference in perception depending on cue type. The greater intervention of intensity stimuli, in terms of disruption in causal interactions, highlights their salience as already emerged through prior behavioral as well as neural studies [ 6 , 7 , 9 , 10 , 12 , 14 , 15 , 23 , 49 , 50 ]. The effects we found from the INSIDEOUT framework, although present, are small in size. This is likely due to our highly specific paradigm (auditory looming bias), rendering the brain states only subtly, but not fundamentally, different. Despite this highly specific approach, though, INSIDEOUT still revealed significant effects in line with previous findings. Contrary to the coarse granularity offered by INSIDEOUT, the fine-grained method of NDTE yielded insights into which regions are the main hubs in manifesting the looming bias, and does so in a data-driven way. By considering all ROIs of a given parcellation, the cross-connectivity is calculated. By, then, ranking regions based on their outflow (sources) or inflow (receivers) and iteratively comparing networks ( Sec. 4.4.2 ), conclusions about ROIs, or networks thereof, with the most essential contribution per considered condition emerge. It should be noted that the timescale of all effects is defined by the calculated minimum of the autocorrelation function. As shown in previous research, this is a solid approach to our investigations [ 29 , 30 ]. In a more ideal way, though, and although computationally significantly more costly, this parameter could be set individually for each considered time series. In our investigation we adhered to the rather coarse parcellation of the Desikan-Killiany atlas [ 32 ]. Our selection relies on both aiming to compare outcomes to prior literature [ 13 , 23 ] as well as reduce complexity, especially in the case of NDTE calculations. Finer parcellations, such as the one from Destrieux [ 51 ], could offer different insights depending on the question at hand; yet they come with higher amount of regions and therefore complexity. Finer parcellations may additionally be more prone to wrongful activity attribution if the precision of the source localization is insufficient. Although we used individualised anatomical information to aid the performed EEG source localization [ 31 ], spatial imprecisions are inevitable. An example thereof is the depth-weighting done by algorithms for sources that are intricately placed on the cortex. Additionally, should activity arise from subcortical surfaces at greater distances from the sensors, EEG may wrongfully attribute the recorded activity. Our results are in good agreement with relevant literature, yet different imaging methods, selected parcellations or implemented algorithms may lead to slightly altered outcomes. 4 Methods 4.1 Participants Thirty-five healthy young adults were invited for study participation. Exclusion criteria during recruitment comprised self-reported indications of psychological and neurological disorders or acute or chronic heavy respiratory diseases that may prevent the participant from sitting still during the EEG recording. All invited participants signed informed consent prior to testing, were neither deceived nor harmed in any way and were informed that they could abort the experiment at any time without any justification or consequences. The study was conducted in accordance with the standards of the Declaration of Helsinki (2000). No additional ethics committee approval was required given the non-medical non-invasive nature of our study, as per the Austrian Universities Act of 2002. In total, experiments lasted around five hours per participant and participants received monetary compensation in return for their time. Prior to conducting the experiment, participants’ hearing thresholds between 1 and 12.5 kHz were measured via pure tone audiometry (Sennheiser HDA200; AGRA Expsuite application [ 52 ]), with a deviation of more than 20 dB from the age mean [ 53 ] leading to subject exclusion. Six participants were excluded (29 remaining subjects, 15 females: 25.0 ± 2.60 years old (mean ± standard deviation); 14 males: 25.1 ± 2.77 years old). An error rate in recognition of static sounds (catch-trials) exceeding 20% resulted in one additional exclusion (female, 45.2% errors). In total, 28 participants were included in this study. They were informed of the procedure and their rights (no deception nor harm, freedom to interrupt the experiment without justification or repercussions) and signed informed consent prior to testing. The study was conducted in accordance with the standards of the Declaration of Helsinki. No additional ethics committee approval was required given the non-medical non-invasive nature of our study, as per the Austrian Universities Act of 2002. Experiments lasted around five hours per subject and participants were remunerated after testing. 4.2 Stimuli The auditory stimuli were complex harmonic tones [ 54 ] ( F 0 = 100 Hz, bandwidth 1 − 16 kHz), filtered with listener-specific head-related transfer functions (HRTFs) to sound as coming from either the right or left direction on the interaural axis when presented over earphones. Stimulus duration was 1.2 s with 10 ms onset and offset ramps of raised-cosine shape. Inter-stimulus intervals lasted 500 ms. Trials were randomized throughout the experiment and balanced over blocks, with 50% looming and 50% receding sounds. Those were created by either modifying the intensity or the spectral shape of a sound and crossfading between the final simulated sound source positions (from far to near for the looming, and near to far for the receding condition). The intensity manipulation resulted in a sound appearing to recede while its intensity decreased with time. We presented sounds crossfading between +2.5 dB (near position) and −2.5 dB (far position) to induce looming and receding sensations (e.g., 11). For changes in spectral shape, we manipulated the individually recorded HRTFs following the procedure introduced in [ 12 ]. The different cue types were applied block-wise. Apart from movement and spatial cue type, we block-wise manipulated whether the sound source was presented from the left or the right side of the listener. The experiment consisted of two parts: an initial passive listening part, during which subjects were watching a silent subtitled movie while being exposed to 600 trials and a subsequent active part, in which subjects performed a spatial discrimination task on the presented sounds. Stimuli and experimental procedures were programmed in MATLAB (R2018b, Mathworks, Natick, Massachusetts) with the use of the Auditory Modeling Tool-box [ 55 ] and Psychtoolbox [ 56 ]. 4.3 Recordings and processing EEG recordings were done with a 128-channel system (actiCAP with actiCHamp; Brain Products GmbH, Gilching, Germany) at a sampling rate of 1 kHz. Noisy channels were being noted during the recordings. All saved EEG data were visually inspected to detect potential additional noisy channels, which were then spherically interpolated. Inspected data were bandpass-filtered between 0.5 − 100 Hz (Kaiser window, β = 7.2, n = 462) and epoched ([−200, 1500] ms) relative to stimulus onset. We applied hard thresholds at −200 and 800 µV to detect and inspect extremely noisy trials. An additional check for the identification of additional bad channels was implemented, via an automatic channel rejection step; detected channels would then be visually inspected and interpolated. No additional noisy channels were detected for any of the subjects at this step. Independent component analysis (ICA) was followed by a manual artifact inspection and rejection of oculomotor artifacts (removal of up to 3 components per subject). The cleaned data were thereafter re-referenced to their average. Within each subject, trials were equalized to match the condition with the minimum amount within the subject after trial rejection. This was achieved within each subject by pseudo-selection aiming to maintain an equal distribution across the recordings. This resulted in an average of 569 clean trials (SD = 27.7) per subject. All preprocessing steps were undertaken on the EEGLAB free software (57) in MATLAB (R2018b, Mathworks, Natick, Massachusetts). Twenty-five (25) of 28 participants had their individual anatomical structures and electrode positions recorded. Anatomical magnetic resonance images (MRIs) were segmented via Freesurfer [ 58 ] and used to create a study protocol on Brainstorm [ 59 ]. For the remaining 3 subjects, the default anatomical models of Brainstorm were used (ICBM152 brain template); individual MRIs could not be recorded due to incompatibilities with the scanner (suspicion of metallic parts in the body). Anatomical models were created via OpenMEEG [ 60 ] with following parameters: boundary element model (BEM) surfaces had 1922 vertices per layer for scalp, outer skull and inner skull, and a skull thickness of 4 mm. The relative conductivity was set to 0.0125 for the outer skull and to 1 for the remaining layers. Manual co-registration between head model and individual electrode locations was done for each subject individually. Recorded activity was inferred to the cortical surface via dynamic statistical parametric mapping (dSPM) [ 61 ]. The noise covariance was calculated from a 200 ms pre-stimulus interval. Dipole orientations were considered constrained to the surface and source signals were reconstructed at 15000 vertices describing the pial surface. Following previous literature [ 13 , 23 ], cortical mapping was done according to the Desikan-Killiany parcellation [ 32 ]. We extracted all ROI time series from the 68 areas of the atlas, as defined in Brainstorm. Based on the evoked time courses at the level of the transverse temporal gyrus, taken from [ 10 ], a time window of 300 ms post-change was defined as the time window of interest. 4.4 Connectivity calculations Our NDTE connectivity analyses, being based on Granger causality, assume stationary signals as input. In order to fulfill this stationarity requirement, we tested our time courses for this property. Following the recommendations of Brainstorm [ 59 ], each time-series was subjected to both the Kwiatkowski-Phillips-Schmidt-Shin test (KPSS) for trend-stationarity and the unit root Augmented Dickey Fuller test (ADF), as implemented in MATLAB 2018b (kpsstest, adftest; Mathworks, Natick, Massachusetts). As broadband EEG signals are highly non-stationary, stationarity of all signals was restored through double differencing of the individual time-series [ 24 ]. 4.4.1 INSIDEOUT The INSIDEOUT framework [ 29 ] is based on the time-shifted correlation matrices between each considered time series and its time-reversed version, thereby echoing the asymmetry in temporal processing. The arrow of time captures the interaction with the environment: a system that remains unperturbed by external factors maintains its intrinsic equilibrium and is therefore characterised by high reversibility. Higher dissimilarity of the forward and reverse time series corresponds to higher non-reversibility, and thereby higher impact of the external environment on the intrinsic dynamics. Reversed time series were obtained by inverting the original ones in time, for each condition, subject, ROI and trial. Correlations between time series were calculated through the MATLAB function corr, for both the forward as well as the reverse time-shifted correlations. If FS forward ( T ) and FS reversal ( T ), expressed as mutual information based on the time-shifted correlations, are the matrices representing the causal dependencies of the system, here across ROIs, the non-reversibility (non-equilibrium) per condition is calculated as and is hence equal to the mean of the absolute squared difference between the forward and reversed matrices (cf. [ 29 ] for detailed calculations). Time-shift T is defined as the decay to the first minimum of the autocorrelation function across conditions and subjects [ 29 , 30 ]. Statistical differences among the conditions were assessed based on ANOVA with the factors of attention, cue type and motion direction. 4.4.2 NDTE The data used in the assessment of NDTE were based on the extracted looming bias on a single trial basis for each ROI. All subsequent calculations were done following the pipeline described by Deco et al. [ 30 ]. In it, the statistical causal interaction between any two ROIs is assessed based on the measure of mutual information. Considering X and Y to denote the activity of the source and target ROIs, respectively, the mutual information is calculated as where I ( Y i +1 ; X i |Y i ) corresponds to the degree of statistical dependence between the source’s past X i = [ X i , X i− 1 , …, X i− ( T − 1) ] and the target’s immediate future Y i +1 [ 30 ]. H ( Y i +1 |Y i ) and H ( Y i +1 |X i , Y i ) express the respective conditional entropies, that are, in the implemented framework, estimated based on covariance matrices [ 62 ]. The time interval defined by T stems from the autocorrelation of the time series. Following Deco et al. [ 30 ], the corresponding order (”maximum lag”) was calculated based on the decay to the first minimum of the autocorrelation function across conditions and subjects; in our case T = 6. In order to be able to compare and combine NDTE values across ROI pairs, calculated connectivity values were normalised as namely by the total mutual information the past of both source X and target Y hold about the future of target Y . This is the quantity considered throughout our calculations and yields, for each trial, an NDTE matrix with the bidirectional flow among all 68 ROIs of the parcellation. As the whole cortical surface and all bidirectional connections therein are considered, the high amount of comparisons is susceptible to spurious correlation outcomes. For that reason, circular-shift surrogate data were generated for each considered ROI pair. P-values for each connection were assessed based on the distribution of connectivity data resulting from 100 independent circular time-shifted surrogate iterations. Statistical significance of connections between ROIs was calculated through p-value aggregation done by Stouffer’s method [ 63 ] in two steps; initially at a subject level with a within-condition aggregation across trials of a subject, and subsequently at a group level, with aggregation across all subjects of a condition. For each considered condition, the multiple comparison correction was performed by the false discovery rate method (FDR) [ 64 ]. The corrected values were then used as a binary “significance mask”, to select the significant connections per condition. The resulting data comprised one NDTE matrix of dimensions ROI x ROI per condition, containing the averaged connectivity values for the ROI pairs that survived the significance evaluations. Inflow ROIs are positioned along the first, while outflow ROIs along the second dimension of the NDTE matrix, termed C All . For each ROI i of C All , the total inflow from all remaining ROIs j of the cortical parcellation is defined as the sum of connectivity across all columns of the matrix: G in ( i ) = Σ j C Alli,j . The respective holds for the total outflow per ROI j : G out ( j ) = Σ i C Alli,j . The major hubs are identified through an iterative process. After sorting the regions based on their inflow (for inflow hubs) or outflow (for outflow hubs), an algorithm searches for the largest subset of ROIs k that have a value G hub significantly larger than any other set, comprising the same amount of regions. The significance value of each G hub ( k ) is assessed via 1000 Monte Carlo simulations, where for each permutation, one member of the current subset k is substituted with any of the remaining ones from the parcellation, and the G hub ( k ) is calculated anew. The in- and outflow values are calculated as where Σ k C Allk is the total flow within the considered subset, Σ k G in ( k ) represents the total inflow to the considered subset from all ROIs of the parcellation and Σ k G out ( k ) the total outflow of the subset to the rest of the ROIs. For the inflow hubs the mutlipliers are [ a = 1; b = 1] and for the outflow hubs [ a = −1; b = −1]. By progressively adding one ROI (”node” in [ 30 ]) to the considered subset ( k = [ i 1 , …, i l ], where l is the whole set of ROIs), the major hubs emerge as the set for which the in- or outflow is still within significance limits (i.e., smaller than 0.05). Data availability Data are available under https://osf.io/4gdy2/ . Code availability Analysis scripts are available and will be updated under https://osf.io/4gdy2/ . Codes for the NDTE framework [ 30 ] are available under https://github.com/decolab/nhb-ndte . Author contributions K.I. and R.Baum. conceived the study. K.I. performed the EEG processing and connectivity analyses. R.Baru. performed the statistical analyses and offered insights during the study design. G.D. contributed to concept refinement and offered consultation throughout the implementations. K.I., R.Baru, R.Baum and B.T. designed the data presentation and wrote the manuscript. All authors revised the manuscript and approved the current version. Conflict of interest The authors declare no competing interests. Acknowledgments We would like to thank Dr. Guenther Koliander for his insightful discussions on the statistical basis and background of the studied approaches, and Elvira Del Agua Banyeres for her support in the computational implementations. This research was funded by the Austrian Science Fund (FWF, I 4294-B and ZK66), the National Research Development and Innovation Office grant (NKFIH; ANN 131305), and the Hungarian Academy of Sciences [Magyar Tudományos Akadémia (MTA)] through the János Bolyai grant (BO/00237/19/2). Footnotes https://osf.io/4gdy2/ References [1]. ↵ Pereira , A. G. & Moita , M. A. Is there anybody out there? Neural circuits of threat detection in vertebrates . Current Opinion in Neurobiology 41 , 179 – 187 , DOI: 10.1016/j.conb.2016.09.011 ( 2016 ). OpenUrl CrossRef PubMed [2]. Deneux , T. , Kempf , A. , Daret , A. , Ponsot , E. & Bathellier , B. Temporal asymmetries in auditory coding and perception reflect multi-layered nonlinearities . Nature Communications 7 , DOI: 10.1038/ncomms12682 ( 2016 ). OpenUrl CrossRef PubMed [3]. Maier , J. X. , Neuhoff , J. G. , Logothetis , N. K. & Ghazanfar , A. A. Multisensory Integration of Looming Signals by Rhesus Monkeys . Neuron 43 , 177 – 181 , DOI: 10.1016/j.neuron.2004.06.027 ( 2004 ). OpenUrl CrossRef PubMed Web of Science [4]. ↵ Ghazanfar , A. A. , Neuhoff , J. G. & Logothetis , N. K. Auditory looming perception in rhesus monkeys . Proceedings of the National Academy of Sciences of the United States of America 99 , 15755 – 15757 , DOI: 10.1073/pnas.242469699 ( 2002 ). OpenUrl Abstract / FREE Full Text [5]. ↵ Maier , J. X. & Ghazanfar , A. A. Looming Biases in Monkey Auditory Cortex . The Journal of Neuroscience 27 , 4093 – 4100 , DOI: 10.1523/JNEUROSCI.0330-07.2007 ( 2007 ). OpenUrl Abstract / FREE Full Text [6]. ↵ Orioli , G. , Bremner , A. J. & Farroni , T. Multisensory perception of looming and receding objects in human newborns . Current Biology 28 , R1294 – R1295 , DOI: 10.1016/j.cub.2018.10.004 ( 2018 ). OpenUrl CrossRef [7]. ↵ Morrongiello , B. A. , Hewitt , K. L. & Gotowiec , A. Infants’ discrimination of relative distance in the auditory modality: Approaching versus receding sound sources . Infant Behavior and Development 14 , 187 – 208 , DOI: 10.1016/0163-6383(91)90005-D ( 1991 ). OpenUrl CrossRef Web of Science [8]. Freiberg , K. , Tually , K. & Crassini , B. Use of an auditory looming task to test infants’ sensitivity to sound pressure level as an auditory distance cue . British Journal of Developmental Psychology 19 , 1 – 10 , DOI: 10.1348/026151001165903 ( 2001 ). OpenUrl CrossRef [9]. ↵ Neuhoff , J. G. Perceptual bias for rising tones . Nature 395 , 123 – 124 , DOI: 10.1038/25862 ( 1998 ). OpenUrl CrossRef PubMed Web of Science [10]. ↵ Ignatiadis , K. et al. Cortical signatures of auditory looming bias show cue-specific adaptation between newborns and young adults . Communications Psychology 2 , 1 – 15 , DOI: 10.1038/s44271-024-00105-5 . OpenUrl CrossRef [11]. ↵ Neuhoff , J. G. An adaptive bias in the perception of looming auditory motion . Ecological Psychology 13 , 87 – 110 , DOI: 10.1207/S15326969ECO1302 2 ( 2001 ). OpenUrl CrossRef Web of Science [12]. ↵ Baumgartner , R. et al. Asymmetries in behavioral and neural responses to spectral cues demonstrate the generality of auditory looming bias . Proceedings of the National Academy of Sciences 114 , 9743 – 9748 , DOI: 10.1073/pnas.1703247114 ( 2017 ). OpenUrl Abstract / FREE Full Text [13]. ↵ Ignatiadis , K. , Baier , D. , Tóth , B. & Baumgartner , R. Neural Mechanisms Underlying the Auditory Looming Bias . Auditory Perception & Cognition 4 , 60 – 73 , DOI: 10.1080/25742442.2021.1977582 ( 2021 ). OpenUrl CrossRef [14]. ↵ Bach , D. R. et al. Rising sound intensity: An intrinsic warning cue activating the amygdala . Cerebral Cortex 18 , 145 – 150 , DOI: 10.1093/cercor/bhm040 ( 2008 ). OpenUrl CrossRef PubMed Web of Science [15]. ↵ Seifritz , E. et al. Neural processing of auditory looming in the human brain . Current Biology 12 , 2147 – 2151 , DOI: 10.1016/S0960-9822(02)01356-8 ( 2002 ). OpenUrl CrossRef PubMed Web of Science [16]. ↵ Olman , C. Chapter 14.2 Central Processing, Sound identity vs. location ( 2022 ). [17]. ↵ Saur , D. et al. Ventral and dorsal pathways for language . Proceedings of the National Academy of Sciences 105 , 18035 – 18040 , DOI: 10.1073/pnas.0805234105 ( 2008 ). Publisher: Proceedings of the National Academy of Sciences. OpenUrl Abstract / FREE Full Text [18]. ↵ Blauert , J. Braasch , J. Majdak , P. , Baumgartner , R. & Jenny , C. Formation of three-dimensional auditory space . In Blauert , J. & Braasch , J. (eds.) The Technology of Binaural Understanding, Modern Acoustics and Signal Processing , 115 – 149 , DOI: 10.1007/978-3-030-00386-9 5 ( Springer International Publishing , Cham , 2020 ). OpenUrl CrossRef [19]. ↵ van der Heijden , K. , Rauschecker , J. P. , de Gelder , B. & Formisano , E. Cortical mechanisms of spatial hearing . Nature Reviews Neuroscience 20 , 609 – 623 , DOI: 10.1038/s41583-019-0206-5 ( 2019 ). Number: 10 Publisher: Nature Publishing Group. OpenUrl CrossRef PubMed [20]. ↵ Warren , J. D. , Zielinski , B. A. , Green , G. G. R. , Rauschecker , J. P. & Griffiths , T. D. Perception of sound-source motion by the human brain . Neuron 34 , 139 – 148 , DOI: 10.1016/s0896-6273(02)00637-2 ( 2002 ). OpenUrl CrossRef PubMed Web of Science [21]. ↵ Stam , C. J. & van Straaten , E. C. W. The organization of physiological brain networks . Clinical Neurophysiology: Official Journal of the International Federation of Clinical Neurophysiology 123 , 1067 – 1087 , DOI: 10.1016/j.clinph.2012.01.011 ( 2012 ). OpenUrl CrossRef PubMed [22]. ↵ Sporns , O. , Honey , C. J. & Kötter , R. Identification and Classification of Hubs in Brain Networks . PLOS ONE 2 , e1049 , DOI: 10.1371/journal.pone.0001049 ( 2007 ). Publisher: Public Library of Science. OpenUrl CrossRef PubMed [23]. ↵ Bidelman , G. M. & Myers , M. H. Frontal cortex selectively overrides auditory processing to bias perception for looming sonic motion . Brain Research 1726 , 146507 , DOI: 10.1016/j.brainres.2019.146507 ( 2020 ). OpenUrl CrossRef [24]. ↵ Seth , A. K. , Barrett , A. B. & Barnett , L. Granger Causality Analysis in Neuroscience and Neuroimaging . The Journal of Neuroscience 35 , 3293 – 3297 , DOI: 10.1523/JNEUROSCI.4399-14.2015 ( 2015 ). OpenUrl FREE Full Text [25]. ↵ Lobier , M. , Siebenhühner , F. , Palva , S. & Palva , J. M. Phase transfer entropy: A novel phase-based measure for directed connectivity in networks coupled by oscillatory interactions . NeuroImage 85 , 853 – 872 , DOI: 10.1016/j.neuroimage.2013.08.056 ( 2014 ). OpenUrl CrossRef PubMed [26]. ↵ Chen , Y. , Bressler , S. L. & Ding , M. Frequency decomposition of conditional Granger causality and application to multivariate neural field potential data . Journal of Neuroscience Methods 150 , 228 – 237 , DOI: 10.1016/j.jneumeth.2005.06.011 ( 2006 ). OpenUrl CrossRef PubMed Web of Science [27]. ↵ Schmidt , C. et al. A Multivariate Granger Causality Concept towards Full Brain Functional Connectivity . PLoS ONE 11 , e0153105 , DOI: 10.1371/journal.pone.0153105 ( 2016 ). OpenUrl CrossRef [28]. ↵ Friston , K. J. , Harrison , L. & Penny , W. Dynamic causal modelling . NeuroImage 19 , 1273 – 1302 , DOI: 10.1016/S1053-8119(03)00202-7 ( 2003 ). OpenUrl CrossRef PubMed Web of Science [29]. ↵ Deco , G. , Sanz Perl , Y. , Bocaccio , H. , Tagliazucchi , E. & Kringelbach , M. L. The INSIDEOUT framework provides precise signatures of the balance of intrinsic and extrinsic dynamics in brain states . Communications Biology 5 , 1 – 13 , DOI: 10.1038/s42003-022-03505-7 ( 2022 ). Number: 1 Publisher: Nature Publishing Group. OpenUrl CrossRef [30]. ↵ Deco , G. , Vidaurre , D. & Kringelbach , M. L. Revisiting the global workspace orchestrating the hierarchical organization of the human brain . Nature Human Behaviour 5 , 497 – 511 , DOI: 10.1038/s41562-020-01003-6 ( 2021 ). Number: 4 Publisher: Nature Publishing Group. OpenUrl CrossRef [31]. ↵ Ignatiadis , K. , Barumerli , R. , Tóth , B. & Baumgartner , R. Effects of individualized brain anatomies and EEG electrode positions on inferred activity of the primary auditory cortex . Frontiers in Neuroinformatics 16 ( 2022 ). [32]. ↵ Desikan , R. S. et al. An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest . NeuroImage 31 , 968 – 980 , DOI: 10.1016/j.neuroimage.2006.01.021 ( 2006 ). OpenUrl CrossRef PubMed Web of Science [33]. ↵ Romei , V. , Murray , M. M. , Cappe , C. & Thut , G. Preperceptual and Stimulus-Selective Enhancement of Low-Level Human Visual Cortex Excitability by Sounds . Current Biology 19 , 1799 – 1805 , DOI: 10.1016/j.cub.2009.09.027 ( 2009 ). OpenUrl CrossRef PubMed Web of Science [34]. ↵ Bach , D. R. , Furl , N. , Barnes , G. & Dolan , R. J. Sustained Magnetic Responses in Temporal Cortex Reflect Instantaneous Significance of Approaching and Receding Sounds . PLOS ONE 10 , e0134060 , DOI: 10.1371/journal.pone.0134060 ( 2015 ). OpenUrl CrossRef PubMed [35]. ↵ Kovach , C. K. et al. Anterior Prefrontal Cortex Contributes to Action Selection through Tracking of Recent Reward Trends . The Journal of Neuroscience 32 , 8434 – 8442 , DOI: 10.1523/JNEUROSCI.5468-11.2012 ( 2012 ). OpenUrl Abstract / FREE Full Text [36]. ↵ Koechlin , E. Frontal pole function: what is specifically human? Trends in Cognitive Sciences 15 , 241 , DOI: 10.1016/j.tics.2011.04.005 ( 2011 ). Publisher: Elsevier. OpenUrl CrossRef PubMed Web of Science [37]. ↵ Farrow , T. et al. Neural correlates of the behavioral-autonomic interaction response to potentially threatening stimuli . Frontiers in Human Neuroscience 6 ( 2013 ). [38]. ↵ Blackwood , N. J. et al. Imaging attentional and attributional bias: an fMRI approach to the paranoid delusion . Psychological Medicine 30 , 873 – 883 , DOI: 10.1017/S0033291799002421 ( 2000 ). Publisher: Cambridge University Press. OpenUrl CrossRef PubMed Web of Science [39]. ↵ Yin , S. , Liu , Y. , Petro , N. M. , Keil , A. & Ding , M. Amygdala Adaptation and Temporal Dynamics of the Salience Network in Conditioned Fear: A Single-Trial fMRI Study . eNeuro 5 , DOI: 10.1523/ENEURO.0445-17.2018 ( 2018 ). Publisher: Society for Neuroscience Section: New Research. OpenUrl Abstract / FREE Full Text [40]. ↵ Herrmann , M. J. et al. Modulation of sustained fear by transcranial direct current stimulation (tDCS) of the right inferior frontal cortex (rIFC) . Biological Psychology 139 , 173 – 177 , DOI: 10.1016/j.biopsycho.2018.10.013 ( 2018 ). OpenUrl CrossRef [41]. ↵ Mothes-Lasch , M. , Becker , M. P. I. , Miltner , W. H. R. & Straube , T. Neural basis of processing threatening voices in a crowded auditory world . Social Cognitive and Affective Neuroscience 11 , 821 – 828 , DOI: 10.1093/scan/nsw022 ( 2016 ). OpenUrl CrossRef PubMed [42]. ↵ Rodríguez , M. , et al. Interoceptive Insular Cortex Mediates Both Innate Fear and Contextual Threat Conditioning to Predator Odor . Frontiers in Behavioral Neuroscience 13 , 283 , DOI: 10.3389/fnbeh.2019.00283 ( 2020 ). OpenUrl CrossRef [43]. ↵ Wang , Q. et al. Insular cortical circuits as an executive gateway to decipher threat or extinction memory via distinct subcortical pathways . Nature Communications 13 , 5540 , DOI: 10.1038/s41467-022-33241-9 ( 2022 ). Number: 1 Publisher: Nature Publishing Group. OpenUrl CrossRef [44]. ↵ Hortensius , R. et al. The Basolateral Amygdalae and Frontotemporal Network Functions for Threat Perception . eNeuro 4 , ENEURO.0314–16.2016 , DOI: 10.1523/ENEURO.0314-16.2016 ( 2017 ). OpenUrl Abstract / FREE Full Text [45]. ↵ Concina , G. , Cambiaghi , M. , Renna , A. & Sacchetti , B. Coherent Activity between the Prelimbic and Auditory Cortex in the Slow-Gamma Band Underlies Fear Discrimination . Journal of Neuroscience 38 , 8313 – 8328 , DOI: 10.1523/JNEUROSCI.0540-18.2018 ( 2018 ). Publisher: Society for Neuroscience Section: Research Articles. OpenUrl Abstract / FREE Full Text [46]. ↵ Baumgartner , R. & Majdak , P. Decision making in auditory externalization perception: model predictions for static conditions . Acta Acustica 5 , 59 , DOI: 10.1051/aacus/2021053 ( 2021 ). Publisher: EDP Sciences. OpenUrl CrossRef [47]. Best , V. , Baumgartner , R. , Lavandier , M. , Majdak , P. & Kopčo , N. Sound externalization: A review of recent research . Trends in Hearing 24 , 2331216520948390 , DOI: 10.1177/2331216520948390 ( 2020 ). Publisher: SAGE Publications Inc. OpenUrl CrossRef [48]. ↵ Kolarik , A. J. , Moore , B. C. J. , Zahorik , P. , Cirstea , S. & Pardhan , S. Auditory distance perception in humans: a review of cues, development, neuronal bases, and effects of sensory loss . Attention, Perception, & Psychophysics 78 , 373 – 395 , DOI: 10.3758/s13414-015-1015-1 ( 2016 ). OpenUrl CrossRef PubMed [49]. ↵ Neuhoff , J. G. Looming sounds are perceived as faster than receding sounds . Cognitive Research 1 , 1 – 9 , DOI: 10.1186/s41235-016-0017-4 ( 2016 ). OpenUrl CrossRef [50]. ↵ Hall , D. A. & Moore , D. R. Auditory Neuroscience: The Salience of Looming Sounds . Current Biology 13 , R91 – R93 , DOI: 10.1016/S0960-9822(03)00034-4 ( 2003 ). OpenUrl CrossRef PubMed Web of Science [51]. ↵ Destrieux , C. , Fischl , B. , Dale , A. & Halgren , E. Automatic parcellation of human cortical gyri and sulci using standard anatomical nomenclature . NeuroImage 53 , 1 – 15 , DOI: 10.1016/j.neuroimage.2010.06.010 ( 2010 ). OpenUrl CrossRef PubMed Web of Science [52]. ↵ Mihocic , M. & Majdak , P. ExpSuite - Browse /Applications/Agra at Source-Forge.net . [53]. ↵ Rodríguez Valiente , A. , Trinidad , A. , García Berrocal , J. R. , Górriz , C. & Ramírez Camacho , R. Extended high-frequency (9-20 kHz) audiometry reference thresholds in 645 healthy subjects . International Journal of Audiology 53 , 531 – 545 , DOI: 10.3109/14992027.2014.893375 ( 2014 ). OpenUrl CrossRef PubMed [54]. ↵ Schroeder , M. Synthesis of low-peak-factor signals and binary sequences with low autocorrelation (Corresp .). IEEE Transactions on Information Theory 16 , 85 – 89 , DOI: 10.1109/TIT.1970.1054411 ( 1970 ). Conference Name: IEEE Transactions on Information Theory. OpenUrl CrossRef [55]. ↵ Majdak , P. , Hollomey , C. & Baumgartner , R. AMT 1.x: A toolbox for reproducible research in auditory modeling . Acta Acustica 6 , 19 , DOI: 10.1051/aacus/2022011 ( 2022 ). Publisher: EDP Sciences. OpenUrl CrossRef [56]. ↵ Brainard , D. H. The Psychophysics Toolbox . Spatial Vision 10 , 433 – 436 , DOI: 10.1163/156856897X00357 ( 1997 ). Publisher: Brill. OpenUrl CrossRef PubMed Web of Science [57]. Delorme , A. & Makeig , S. EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis . Journal of Neuroscience Methods 134 , 9 – 21 , DOI: 10.1016/j.jneumeth.2003.10.009 ( 2004 ). OpenUrl CrossRef PubMed Web of Science [58]. ↵ Fischl , B. FreeSurfer . NeuroImage 62 , 774 – 781 , DOI: 10.1016/j.neuroimage.2012.01.021 ( 2012 ). OpenUrl CrossRef PubMed Web of Science [59]. ↵ Tadel , F. , Baillet , S. , Mosher , J. C. , Pantazis , D. & Leahy , R. M. Brainstorm: A User-Friendly Application for MEG/EEG Analysis . Computational Intelligence and Neuroscience 2011 , 1 – 13 , DOI: 10.1155/2011/879716 ( 2011 ). OpenUrl CrossRef [60]. ↵ Gramfort , A. , Papadopoulo , T. , Olivi , E. & Clerc , M. OpenMEEG: opensource software for quasistatic bioelectromagnetics . BioMedical Engineering OnLine 9 , 45 , DOI: 10.1186/1475-925X-9-45 ( 2010 ). OpenUrl CrossRef PubMed [61]. ↵ Dale , A. M. et al. Dynamic Statistical Parametric Mapping: Combining fMRI and MEG for High-Resolution Imaging of Cortical Activity . Neuron 26 , 55 – 67 , DOI: 10.1016/S0896-6273(00)81138-1 ( 2000 ). OpenUrl CrossRef PubMed Web of Science [62]. ↵ Brovelli , A. , Chicharro , D. , Badier , J.-M. , Wang , H. & Jirsa , V. Characterization of Cortical Networks and Corticocortical Functional Connectivity Mediating Arbitrary Visuomotor Mapping . The Journal of Neuroscience: The Official Journal of the Society for Neuroscience 35 , 12643 – 12658 , DOI: 10.1523/JNEUROSCI.4892-14.2015 ( 2015 ). OpenUrl Abstract / FREE Full Text [63]. ↵ Stouffer , S. A. , Suchman , E. A. , Devinney , L. C. , Star , S. A. & Williams Jr. , R. M. The American soldier: Adjustment during army life. (Studies in social psychology in World War II), Vol. 1 . The American soldier: Adjustment during army life. (Studies in social psychology in World War II), Vol. 1 ( Princeton Univ. Press , Oxford, England , 1949 ). Pages: xii , 599. OpenUrl [64]. ↵ Benjamini , Y. & Hochberg , Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing . Journal of the Royal Statistical Society. Series B (Methodological) 57 , 289 – 300 ( 1995 ). Publisher: [Royal Statistical Society, Wiley]. OpenUrl CrossRef PubMed Web of Science View the discussion thread. Back to top Previous Next Posted June 24, 2024. Download PDF Data/Code Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Threat-Related Corticocortical Connectivity Elicited by Rapid Auditory Looms Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Threat-Related Corticocortical Connectivity Elicited by Rapid Auditory Looms Karolina Ignatiadis , Roberto Barumerli , Gustavo Deco , Brigitta Tóth , Robert Baumgartner bioRxiv 2024.06.21.600106; doi: https://doi.org/10.1101/2024.06.21.600106 Share This Article: Copy Citation Tools Threat-Related Corticocortical Connectivity Elicited by Rapid Auditory Looms Karolina Ignatiadis , Roberto Barumerli , Gustavo Deco , Brigitta Tóth , Robert Baumgartner bioRxiv 2024.06.21.600106; doi: https://doi.org/10.1101/2024.06.21.600106 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Neuroscience Subject Areas All Articles Animal Behavior and Cognition (7644) Biochemistry (17726) Bioengineering (13916) Bioinformatics (42033) Biophysics (21486) Cancer Biology (18635) Cell Biology (25549) Clinical Trials (138) Developmental Biology (13397) Ecology (19940) Epidemiology (2067) Evolutionary Biology (24361) Genetics (15620) Genomics (22541) Immunology (17763) Microbiology (40468) Molecular Biology (17207) Neuroscience (88739) Paleontology (667) Pathology (2842) Pharmacology and Toxicology (4834) Physiology (7659) Plant Biology (15175) Scientific Communication and Education (2047) Synthetic Biology (4304) Systems Biology (9834) Zoology (2272)