Introduction

The capacity to flexibly update our memories in light of new information is fundamental to maintaining an accurate model of the world around us. This flexibility requires adaptable memory networks that can be reconfigured upon acquiring new insights. Previous research provided direct evidence for insight-induced reconfigurations of memory representations and showed that insight into the connection of initially separate events propels the integration of these events into coherent episodes (Collin et al., 2015; Milivojevic et al., 2015). Such mnemonic integration allows novel inferences (Spalding et al., 2018; Zeithamova et al., 2012) that aid efficient navigation (Coutanche et al., 2013; Fernandez et al., 2023; He et al., 2022) and decision-making (Boorman et al., 2021; Kumaran et al., 2009; Shohamy & Daw, 2015). Importantly, in everyday life, the inference about the relationship between seemingly unrelated events is often not inferred via direct observation but through imagination. For instance, when reading a book, we gain insight into the plot and possible twists through our imagination, which then prompts us to update our memory representations. Even when new insights are derived from direct observation, the integration process requires imaginative capacities to bind the previously separate memories into a coherent narrative. At the neural level, the hippocampus has been shown to play a pivotal role in (imagination-based) mnemonic integration (Cohn-Sheehy et al., 2021; Collin et al., 2015; Griffiths & Fuentemilla, 2020; Grob et al., 2023a; Milivojevic et al., 2015). However, while the hippocampus appears to be crucial for mnemonic integration, it does not act in isolation but operates in collaboration with cortical areas to accomplish this complex process (Backus et al., 2016; Milivojevic et al., 2015; Pehrs et al., 2018; Schlichting & Preston, 2015; Spalding et al., 2018). Yet, our understanding of the specific areas implicated in the insight-driven reconfiguration of memory representations, beyond the hippocampus, remains limited. Moreover, existing data on the neural underpinnings of mnemonic integration are mainly correlational in nature and which areas are causally involved in the integration of initially unrelated memories into cohesive representations is completely unknown.

One promising candidate that may contribute to insight-driven memory reconfiguration is the angular gyrus. The angular gyrus has extensive structural and functional connections to many other brain regions (Petit et al., 2023), including the hippocampus (Coughlan et al., 2023; Thakral et al., 2020; Uddin et al., 2010; Wang et al., 2014). Accordingly, the angular gyrus has been implicated in a myriad of cognitive functions, including mental arithmetic, visuospatial processing, inhibitory control, and theory-of-mind (Cattaneo et al., 2009; Grabner et al., 2009; Lewis et al., 2019; Schurz et al., 2014). Moreover, there is accumulating evidence pointing to a key role of the angular gyrus in long-term memory (Bellana et al., 2017; Bonnici et al., 2018; Kwon et al., 2022; Wang et al., 2014) and imagination (Ramanan et al., 2018; Thakral et al., 2017, 2020). How these putative functions of the angular gyrus relate to one another, however, remained unclear. We reasoned that these functions might be directly linked, enabling the angular gyrus to drive the integration of (imagination-related) insights into long-term memory. In line with this idea, recent theories propose that the angular gyrus acts as dynamic buffer for spatiotemporal representations (Humphreys et al., 2021), which may allow the angular gyrus to transiently maintain the initially separate events and to integrate these into cohesive narratives. This buffering function of the angular gyrus may be particularly relevant for imagination-based linking. Thus, we hypothesized that the angular gyrus plays a crucial role in integrating imagination-related insights into long-term memory and hence in the dynamic reconfiguration of memory representations in light of new information.

To test this hypothesis and determine the causal role of this area in insight-related memory reconfigurations, we conducted a preregistered study combining a life-like video-based narrative-insight task (NIT; Milivojevic et al., 2015; Fig. 1), probing insight-related reconfigurations of memory, with representational similarity analysis of EEG data and (double-blind) ‘neuro-navigated’ TMS to an area of the left angular gyrus that was implicated in imaginative processing before (Thakral et al., 2017). We predicted that cTBS to the left angular gyrus would reduce the impact of (imagination-based) insight into the link of initially unrelated events on memory performance as well as the insight-driven neural reconfiguration of linked and non-linked event representations.

Modified narrative-insight task (NIT) and procedure. During the pre-phase, participants viewed video events (A, B, and X) from ten different storylines. Each event was preceded by a title (1s) and repeated 18 times. The inter- stimulus interval (ISI) was ∼ 1000 ms. The subsequent insight-phase consisted of two parts. In one part, participants gained insight through a written imagination instruction (I) interspersed with a control instruction (CI). In the other part, they gained insight through a linking video (L) interspersed with a control video (C). The order of gaining insight through imagination or video observation was counterbalanced across participants. Before each insight part, participants received, depending on the experimental group, either a sham or cTBS stimulation over the left angular gyrus (MNI: −48, −67, 30). After the insight-phase, participants had a 30-minute break and then completed a free recall for a maximum of 20 minutes in a different room. In the post-phase, all video events were presented in the same manner as the pre-phase.

Results

Inhibition of the angular gyrus reduces insight-related memory boost

The angular gyrus has been implicated in a myriad of tasks and functions, including long-term memory (Bonnici et al., 2018; Kwon et al., 2022; Wang et al., 2014) and imagination (Ramanan et al., 2018; Thakral et al., 2017). Here we hypothesized that these functions of the angular gyrus are directly linked to one another. Specifically, we postulated that the angular gyrus plays a crucial role in the integration of imagination-related insights into long-term memory representations and that it thus represents a key player in the dynamic reconfiguration of memory in light of new information. To test this hypothesis and the causal role of the angular gyrus in insight-related memory reconfigurations, we combined the life-like video-based NIT with representational similarity analysis of EEG data and (double-blind) neuro-navigated TMS over the left angular gyrus. During the NIT, participants first saw three video events (A, B, and X; pre-phase), which were then either linked into a narrative (A and B) or not (A and X) in a subsequent insight-phase. Critically, before the insight-phase, we applied either sham stimulation (31 participants, 15 females) or inhibitory continuous theta burst stimulation (cTBS; 34 participants, 16 females) to the left angular gyrus. Following the insight-phase and a 30-minute break to mitigate potential TMS aftereffects (Huang et al., 2005; Jannati et al., 2023), participants completed a free recall task, which provided a measure of insight-related changes in subsequent memory. Thereafter, participants saw the same video events (A, B, and X) again in a post-phase. EEG was measured during all stages of the NIT. Contrasting neural representation patterns from the pre- and post-phases allowed us to assess insight-related memory reconfiguration and its modulation by angular gyrus inhibition. Due to its specific relevance in imaginative processes (Ramanan et al., 2018; Thakral et al., 2017, 2020), we expected that the angular gyrus would be particularly relevant if insight relies strongly on imagination. Therefore, participants gained insight into half of the stories by imagining the link themselves, while they observed the link as a video in the other half of the stories. Participants’ ratings showed that they adhered well to these instructions during the linking phase. When linking events via imagination, they reported imagining the linking events very well (M = 3.38, SD = 0.47) and their imagination as depictive (M = 3.35, SD = 0.46). When linking via observation, they reported a high level of understanding of the linking events (M = 3.37, SD = 0.51) and found the linking events meaningful (M = 3.35, SD = 0.52) on a 1-4 Likert scale. Furthermore, participants demonstrated a high level of attention throughout the NIT, responding to target stimuli with near-ceiling performance (M = 99.25 %; SD = 1.40 %) without any group differences (t(63.00) = 0.42, p = 0.675, d = −0.10).

Importantly, participants were unaware of the allocation to the cTBS or sham condition, as indicated by the treatment guess at the end of the experiment (Fisher’s exact test; p = 0.597). Furthermore, TMS stimulation did not affect participants’ subjective mood, wakefulness or arousal (mood: group × time: F(1, 63) = 0.76, p = 0.386, ηG = 0.00; wakefulness: group × time: F(1, 63) = 0.01, p = 0.921, ηG = 0.00; arousal: group × time: F(1, 63) = 0.01, p = 0.921, ηG = 0.00).

As expected, all participants gained insight into which events were linked in the NIT, as they rated the belongingness of linked events higher than non-linked events from pre- to post-insight, as indicated by a linear mixed model (LMM: time × link: β = 2.49, 95% CI [2.15, 2.83], t(418.44) = 14.01, p < 0.001; Figure 2 – Figure supplement 1). Post-hoc tests showed increasing belongingness ratings for linked events and decreasing belongingness ratings for non-linked events from pre- to post-insight (LMM: link: β = 1.49, 95% CI [1.33, 1.64], t(418) = 24.49, p < 0.001; non-link: β = −.91, 95% CI [−1.07, −.76], t(418.00) = −15.03, p < 0.001). This insight was further reflected in the multi-arrangements task (MAT), in which participants were instructed to arrange representative images (A, B, and X) from each story based on their relatedness. In this task, all participants arranged linked events closer together than non-linked events (MAT; LMM: link: β = −1.33, 95% CI [−1.59, −1.07], t(177.00) = −9.81, p < 0.001; Figure 2 – Figure supplement 2). The strong insight gained by all participants was further reflected in their near-ceiling performance in the forced-choice recognition task, in which participants were instructed to identify the event (B or X) that was linked with A. Participants accurately indicated whether B or X was linked to A (sham: M = 94.65%, SD = 9.00%; cTBS: M = 97.34%, SD = 6.10%; Figure 2 – Figure supplement 3). Importantly, there were no group differences in any of these measures (LMM: NIT: group × time × link: β = −0.02, 95% CI [−0.49, 0.45], t(418.27) = −0.09, p = 0.929; LMM: MAT: group × link: β = 0.11, 95% CI [−0.26, 0.48], t(177.00) = 0.58, p = 0.561; LMM: Forced-choice recognition: group: β = 0.23, 95% CI [−0.26, 0.72], t(113.37) = 0.90, p = 0.368), indicating that all participants successfully gained insight into which events were linked and that the (left) angular gyrus did not play a critical role in the process of gaining insight itself.

To investigate the causal role of the left angular gyrus in insight-related episodic memory integration, the key question of this study, we first analyzed the detailedness of participants’ memory for both linked and non-linked events during free recall. Across groups, linked events were generally recalled in more detail than non-linked events (LMM: link: β = 1.20, 95% CI [0.86, 1.54], t(406.00) = 6.75, p < 0.001), suggesting a memory boost for integrated narratives. Most interestingly, inhibiting the left angular gyrus through cTBS reduced this insight-related memory boost for linked events significantly (LMM: group × link: β = −0.54, 95% CI [−1.02, −0.06], t(406.00) = −2.17, p = 0.030; Fig. 2A). Pairwise comparisons revealed a significantly lower number of recalled details for linked events in the cTBS compared to the sham group, while there was no significant difference for non-linked events (LMM: link: β = −0.40, 95% CI [−0.79, −0.02], t(85.50) = −2.79, p = 0.033; LMM: non-link: β = −0.09, 95% CI [−0.29, 0.42], t(406.00) = −0.62, p = 0.926). Additionally, we observed that all participants showed better memory for central compared to peripheral details of the plot when recalling linked events, which was not observed to the same extent for non-linked events (LMM: link × detail: β = 0.61, 95% CI [0.12, 1.09], t(406.00) = 2.43, p = 0.016).

Behavioral results. A, Significantly reduced recall of details for linked events in the cTBS group compared to the sham group, with no significant difference for the non-linked events. B, Significantly reduced recall of the linking events in the cTBS group compared to the sham group. C, Schematic overview of electric field modeling: Simulation was performed for the angular gyrus coordinate (MNI: x = −48, y = −67, z = 30) using a Magstim 70mm figure-of eight coil at 80 % of individual motor thresholds, reflecting the applied setup. The resulting electric field was averaged within a 10mm spherical ROI and centered on the target coordinate and extracted for subsequent analyses. D, Significantly reduced number of details recalled for linked events specifically in the high cTBS group (based on a median-split on simulated electric field strengths). E, Significantly reduced recall of the linking events specifically in the high cTBS group (based on a median-split on simulated electric field strengths). Boxplots show the median for each group. Boxplot whiskers extend to the minimum or maximum value within 1.5 times the interquartile range. Points within the boxplot indicate individual data points per each group. Density plots indicate data distribution per group. The belongingness ratings for the linked and non-linked events are shown in Figure 2 – Figure supplement 1, the data of the multiple arrangements task in Figure 2 – Figure supplement 2, and the data of the forced-choice recognition test in Figure 2 – Figure supplement 3. Statistical differences stem from pairwise post-hoc tests of marginal means. *p < 0.05, ***p < 0.001.

In a second step, we analyzed whether cTBS to the angular gyrus affected, in addition to memory detailedness for initially separate but now linked events, also the memory for the linking events themselves. Our results showed that cTBS (vs. sham) significantly reduced the frequency with which participants recalled the linking events (LMM: group: β = −0.66, 95% CI [−1.13, −0.18], t(98.13) = −2.71, p = 0.008; Fig. 2B). Interestingly, this TMS effect appeared to be particularly pronounced when events were linked via imagination (cTBS vs. sham: t(61.46) = −2.53, p = 0.014, d = −0.63) and was less prominent when they were linked via direct observation (cTBS vs. sham: t(58.59) = − 1.63, p = 0.107, d = −0.40), although it is important to note that the interaction was not significant (LMM: group × mode: β = 0.30, 95% CI [−0.18, 0.77], t(62) = 1.23, p = 0.225).

To assess the effect of cTBS stimulation on the angular gyrus (Pizem et al., 2022; Zhang et al., 2022), we performed electric field simulations at 80 % of the individual motor threshold, averaging the estimated field strength within a 10mm sphere centered around the angular gyrus coordinate (MNI: −48, −67, 30). In order to examine whether the behavioral effects were dependent on the simulated electric field strength (Fig. 2C), we next included electric field strength (strong vs. weak via median split) and repeated the previous linear mixed model predicting the number of details for linked events including a group factor reflecting stimulation strength (sham, low, high). This model yielded a significant group × link interaction (LMM: β = −0.78, 95% CI [−1.35, −0.21,], t(399.00) = −2.63, p = 0.009; Fig. 2D), suggesting a dependency of memory on stimulation strength. Pairwise comparisons for linked events confirmed that a stronger electric field induction in the angular gyrus significantly reduced the memory boost for linked events, while there was no such effect for weak cTBS stimulation (LMM: sham vs. low: β = 0.05, 95% CI [−0.45, 0.55], t(87.70) = 0.28, p = 1.000; sham vs. high: β = 0.74, 95% CI [0.25, 1.23], t(87.7) = 4.45, p < 0.001, low vs. high: β = 0.69, 95% CI [0.13, 1.26], t(87.7) = 3.60, p = 0.007). We further included the electric field strength (strong vs. weak via median split) and repeated the previous linear mixed model predicting the naming of the linking events including the group factor stimulation strength (sham, low, high). This analysis yielded a significant effect of group (LMM: β = −0.92, 95% CI [−1.50, −0.35], t(97.79) = −3.11, p = 0.003; Fig. 2E), suggesting that the memory for the linking events was dependent on the angular gyrus stimulation strength.

Angular gyrus inhibition disrupts neural pattern reconfiguration following imagination-based insight

Our behavioral data showed that cTBS to the angular gyrus reduced the insight-related memory boost. In a next step, we tested whether inhibitory stimulation of the angular gyrus may also alter the insight-related reconfiguration of neural memory representations, taking the mode of insight (i.e. imagination vs. observation) into account. To this end, we leveraged representational similarity analysis (RSA) of EEG data and compared changes in multivariate oscillatory theta power patterns for linked and non-linked events from pre- to post-insight (Fig. 3A). We focused exclusively on the theta band since theta has been shown to hold a key role in episodic memory integration (Backus et al., 2016; Nicolás et al., 2021). For this analysis, similarity maps (time × time) were computed by correlating story-specific theta frequency patterns within linked (A with B) and within non-linked (A with X) events in the pre- and post- phase, separately.

Representational pattern changes. A, Conceptual overview of the representational similarity analysis (RSA) on theta oscillations. First, time-frequency data was computed, and the theta power values (4-7 Hz) were extracted. Using these feature vectors, Pearson’s correlations were computed to compare the power patterns across time points of events (here: event A and B). These correlations resulted in a time × time similarity map. B, Significant cluster, denoted by white dotted line for illustrative purposes, for the change from post-pre and imagination-observation between the cTBS and sham groups using an independent sample cluster-based permutation t-test for linked events (A and B). In the middle panel, follow-up tests on stories linked via imagination revealed increased similarity for the sham group, while no significant effect was observed for the cTBS group. In the lower panel, follow-up tests on stories linked via observation showed decreased similarity for the sham group and increased similarity for the cTBS group. C, Significant cluster, denoted by white dotted line for illustrative purposes, for the change from post-pre and imagination-observation between the cTBS and sham groups using an independent sample cluster-based permutation t-test for non-linked events (A and X). In the middle panel, follow-up tests on stories linked via imagination revealed increased similarity for the cTBS group, while no significant effect was observed for the sham group. In the lower panel, follow-up tests on stories linked via observation showed decreased similarity for the cTBS group and no significant effect for the sham cTBS group. Boxplots show the median similarity for each group at each time point. Boxplot whiskers extend to the minimum or maximum value within 1.5 times the interquartile range. Points within the boxplot indicate individual data points in each group. Density plots indicate data distribution per group and time. *p < 0.05, **p < 0.01, ***p < 0.001.

We examined insight-induced effects on neural representations for linked events by comparing the change from pre- to post-insight (post - pre) and the difference between imagination and observation (imagination – observation) between cTBS and sham groups using an independent sample cluster-based permutation t-test. This analysis yielded a negative cluster (p = 0.032, ci-range = 0.00, SD = 0.00) in the parieto-temporal region (electrodes: T7, Tp7, P7; Fig. 3B). Follow-up tests on the extracted similarity cluster analyzed the representational pattern change and its modulation by TMS separately for the imagination and observation condition. For stories linked via imagination, we obtained an increase in representational similarity from pre- to post-insight in the sham group (t(30) = 3.48, p = 0.002, drepeated measures = 0.62), whereas there was no such increase and even a trend for a decrease in representational similarity for linked events from pre- to post-insight in the cTBS group (t(30) = −2.01, p = 0.053, drepeated measures = −0.36; group × time: F(1, 60) = 14.03, p < 0.001, ηG= 0.09; Fig. 3B middle panel). Interestingly, we observed that a lower change (post - pre) in representational similarity of events linked via imagination (vs. observation) was associated, across groups, with a reduced probability of recall of the linking events (r = 0.27, t(59) = 2.17, p = 0.034), suggesting a direct association between neural pattern reconfiguration and subsequent memory.

For stories that were linked via observation, we observed a seemingly opposite pattern (group × time: F(1, 60) = 19.21, p < 0.001, ηG= 0.12): decreased similarity in the sham group (t(30) = −3.94, p < 0.001, drepeated measures = −0.62) but increased representational similarity in the cTBS group (t(30) = 2.30, p = 0.029, drepeated measures = 0.62; Fig. 3B lower panel). However, these changes in representational similarity for the observation condition should be interpreted with caution, as these seemingly opposite changes appeared to be at least in part driven by group differences already in the pre-phase, before participants gained insight.

Interestingly, we observed a different pattern of insight-related representational pattern changes for non-linked events. Similarly to linked events, we compared the change from pre- to post-insight and the difference between imagination and observation between cTBS and sham using an independent sample cluster-based permutation t-test. This analysis yielded a positive cluster (p = 0.035, ci-range = 0.00, SD = 0.00) in a fronto-temporal region (electrode: FT7; Fig. 3C). Again, we pursued this effect with separate follow-up tests for the imagination and observation conditions. In the imagination condition, the sham group did not show any representational changes for non-linked events (t(30) = −1.35, p = 0.187; drepeated measures = −0.23), while we observed increased neural similarity for non-linked events from pre- to post- insight in the cTBS group (t(30) = 3.61, p = 0.001, drepeated measures = 0.67; Fig. 3C middle panel). Thus, participants who received cTBS to the angular gyrus showed a pattern of pre- to post- insight representational changes for non-linked events that resembled the pattern observed in the sham group for events linked via imagination, suggesting that angular gyrus inhibition before gaining imagination-based insight interfered with efficient pruning of the integrated narrative.

For stories linked via observation, we observed, again, a seemingly opposite pattern (group × time: F(1, 60) = 10.32, p = 0.002, ηG= 0.07): no representational change in the sham group (t(30) = 1.65, p = 0.110, drepeated measures = 0.34) but decreased neural similarity for non-linked events from pre- to post-insight in the cTBS group (t(30) = −2.40, p = 0.023, drepeated measures = − 0.42). Again, these representational changes should be interpreted with caution, as the differences appear to be at least in part driven by group differences in the pre-insight phase.

Inhibiting the angular gyrus diminishes fronto-temporal connectivity associated with imagination-based insight

To further elucidate the neural mechanisms involved in the changes in insight-related memory reconfiguration after cTBS to the left angular gyrus, we next examined changes in functional connectivity using the same contrast as in the previous RSA analyses. More specifically, we computed imaginary coherence for the mean theta frequency using a sliding window approach and tested the change in connectivity from pre- to post-insight for linking via imagination vs. observation between the cTBS and the sham groups using an independent sample cluster-based permutation t-test.

This analysis yielded a negative cluster (p = 0.044, ci-range =0 .00, SD = 0.00; time window: 1.25 – 1.75s; Fig. 4A) between centro-temporal and frontal regions (C4 – Fp1; C6 – Fp1; T8- Fp1; T8 – AF7). For stories linked via imagination, follow-up tests indicated decreased functional connectivity between these regions in the cTBS group (t(30) = −4.25, p < 0.001, drepeated measures = −0.70), while there was no change in the sham group (t(30) = 0.02, p = 0.987, drepeated measures = 0.00; group × time: F(1, 60) = 8.05, p = 0.006, ηG= 0.050; Fig. 4B). Interestingly, across both groups higher coherence between these areas in the post-relative to the pre-phase for stories linked via imagination relative to observation was associated with better recall of details (central and peripheral) for stories linked via imagination (r = 0.31, t(59) = 2.50, p = 0.015; Fig. 4C), suggesting that the reduced crosstalk between these regions was linked to impaired subsequent memory.

Connectivity change for linked events. A, Conceptual depiction of the negative cluster (centro-temporal and frontal). In the upper panel, the connections are presented separately for each electrode pair. The lower panel illustrates the cluster coherence based on the electrode locations relative to the brain. B, Follow-up tests for events linked via imagination indicated decreased coherence between centro-temporal and frontal electrodes for the cTBS group, while no change was observed for the sham group. C, Significant positive correlation, indicating that the less the coherence between centro-temporal and frontal electrodes changed from pre to post for imagination (vs. observation), the fewer details for linked events via imagination were recalled. Please note, that this correlation became even stronger when the outlier was removed (r = 0.38, t(58) = 3.11, p = 0.003). D, Follow-up tests for events linked via observation indicated decreased coherence centro-temporal and frontal electrodes for the sham group, while the cTBS group showed a significant increase in coherence. Boxplots show the median coherence for each group at each time point. Boxplot whiskers extend to the minimum or maximum value within 1.5 times the interquartile range. Points within the boxplot indicate individual data points in each group. Density plots indicate data distribution per group and time. *p < 0.05, **p < 0.01, ***p < 0.001.

For stories linked via observation, we found a seemingly opposite pattern (group × time: F(1, 60) = 12.73, p < .001, ηG= .080; Fig. 4D): decreased functional connectivity for the sham group (t(30) = −2.75, p = 0.010, drepeated measures = −0.47), while the cTBS group exhibited increased connectivity (t(30) = 2.28, p = 0.030, drepeated measures = 0.36). Since these differences appeared to be again already present prior to gaining insight, the functional connectivity changes obtained for the observation condition should be interpreted with caution. Regarding non-linked events (X), we did not find any significant cluster in this coherence analysis (all clusters p > 0.221), indicating that the reported connectivity changes were specific to linked events.

Control variables

Overall, levels of subjective chronic stress, anxiety, and depressive mood were relatively low and not different between groups. The groups did further not differ in terms of personality traits or imagination capacity (all p > 0.056; see Table 1).

Control variables.

Discussion

Updating our memory representations in light of new information is key to keeping an accurate model of the world. Given the previously described role of the angular gyrus in episodic memory and imagination (Benoit & Schacter, 2015; Bonnici et al., 2018; Thakral et al., 2017), we hypothesized that this area plays a pivotal role in (imagination-based) insight-driven memory reconfiguration. To probe the causal role of the angular gyrus in insight-induced mnemonic changes and related neural pattern reconfigurations, we combined a life-like narrative insight task with cTBS to the angular gyrus. Our results show that inhibition of the angular gyrus reduced the insight-related memory boost for linked events and decreased memory for the linking events themselves. At the neural level, angular gyrus inhibition reduced the coupling between centro-temporal and frontal regions and abolished insight-induced neural representational changes for events linked via imagination, which points to impaired memory integration.

Unsurprisingly, all participants successfully gained insight into which events were linked in the NIT, as reflected in their reliable recognition of the linked events in the forced-choice recognition task and their relatedness ratings in the MAT. Importantly this insight had a direct impact on memory, with better memory for linked relative to non-linked events. This memory boost for linked events aligns with previous findings (Cohn-Sheehy et al., 2021; Grob et al., 2023a, 2023b; Wang et al., 2015) and supports the notion that the brain stores episodic memories as coherent narratives (Tulving, 1983), for which any element can cue the entire episode (Horner et al., 2015; Nakazawa et al., 2002). Crucially, while cTBS to the angular gyrus had no effect on basic insight, inhibiting the angular gyrus prior to gaining insight into the relationship of previously separate events specifically reduced this memory boost for linked events, particularly when electric field stimulation of the angular gyrus was strong. The diminished memory boost for linked events following effective angular gyrus inhibition demonstrates a causal involvement of the angular gyrus in prioritizing linked narratives in memory and, by implication, the reconfiguration of memories in light of new information.

Inhibition of the angular gyrus prior to gaining insight impaired memory not only for the linked events but also for the linking events themselves. Again, this was exclusively observed for inhibitory stimulation with a strong electric field, ensuring effective inhibition of the angular gyrus. The resulting impairment in memory for linking events following angular gyrus inhibition highlights its causal role in integrating newly acquired information into pre-existing memory representations, suggesting that the linking events – whether observed or imagined – serve as binding information for creating a coherent narrative. This interpretation is in line with a recent study showing that online TMS over the angular gyrus during a reading-based task impaired the integration of contextual information (Branzi et al., 2021). Dovetailing with recent accounts proposing the angular gyrus as a spatiotemporal buffering region for integrating our continuous stream of experiences (Humphreys et al., 2021), inhibition of the angular gyrus may have disrupted the reactivation and maintenance of initially separate events A and B during the linking phase. This disruption may have impeded the binding process that is likely crucial for the observed memory boost. At the same time, the lacking reactivation and integration of previously separate events A and B during the linking phase may have diminished the significance of the linking event, leading to reduced memory for the linking event itself.

Consistent with previous neuroimaging studies demonstrating distinct representational changes associated with insight into the relationship of initially separate events (Collin et al., 2015; Milivojevic et al., 2015), we observed, in addition to insight-related memory changes, increased similarity in neural representations of events linked via imagination in the theta band from pre- to post-insight. This finding aligns with studies linking theta band activity to hippocampal involvement in memory integration (Backus et al., 2016; Herweg et al., 2020; Nicolás et al., 2021). Mechanistically, hippocampal theta oscillations may facilitate memory integration by promoting more accurate representations of stimulus-specific information while gaining insight (Pacheco Estefan et al., 2021). Critically, inhibiting the angular gyrus abolished these insight-induced changes in neural pattern similarity between linked events in the theta band, which was directly associated with reduced recall of the linking events. Beyond its causal involvement in the insight-driven neural reconfiguration of linked event representations, angular gyrus inhibition led to representational changes for non-linked events that closely resembled the changes observed in the sham group for linked events. This finding indicates that inhibition of the angular gyrus impairs the separation of linked and non-linked events, suggesting its causal role in effectively pruning out events that are not part of the integrated narrative.

Given the consistent findings linking the angular gyrus to scene imagination (Addis et al., 2007; Hassabis et al., 2007; Ramanan et al., 2018; Thakral et al., 2017), our study focused on the role of the angular gyrus in imagination-based mnemonic integration. It is important to note that the neural findings discussed thus far exclusively pertain to imagination-based insight, aligning with the primary focus of our study. When gaining insight via imagination, participants had to retrieve memory representations of initially separate events to construct the linking event in their mind (Schacter et al., 2008). We assume that the retrieval and reinstatement necessary to construct the imagined linking event led to a high degree of co-activation of the neural patterns associated with each event, leading to more similar representations (Wammes et al., 2022). Most interestingly, inhibiting the angular gyrus’ buffering function (Humphreys et al., 2021) disrupted the observed increase in similarity for events linked via imagination, lending further support to its causal role in imagination-based memory integration. Together, our findings demonstrate that inhibiting the angular gyrus disrupts neural pattern reconfiguration following imagination-based insight.

Although we specifically targeted the left angular gyrus and identified its causal role in insight-driven memory reconfiguration, the angular gyrus does not act in isolation. As a major connector hub, it integrates information from various brain regions, with connections including the prefrontal and parietal regions (Frey et al., 2008; Makris et al., 2005, 2007), striatum (Petit et al., 2023; Uddin et al., 2010), sensory-motor areas (Bonner et al., 2013), and the medial-temporal lobe, including the hippocampus (Uddin et al., 2010). Some of these regions play a key role in insight-induced memory reconfigurations, in particular the hippocampus and prefrontal cortex (Collin et al., 2015; Milivojevic et al., 2015). Consistent with this, inhibiting the angular gyrus resulted in decreased theta band connectivity between centro-temporal and frontal regions for events linked via imagination. Importantly, this reduction in connectivity was directly related to a subsequent memory decline. This finding, specific to linked events, aligns with prior research indicating that stronger (and inter-hemispheric) functional coupling in the theta band between frontal and posterior regions is linked to enhanced associative memory (Cruzat et al., 2021; Summerfield & Mangels, 2005; Wu et al., 2007). Furthermore, in line with our finding that is specific to imagination-based linking, the importance of theta band synchronization has been highlighted in research focused on mental imagery (Li et al., 2009). These results suggest that by disrupting the angular gyrus, long-range theta synchronization for events linked via imagination is also disrupted, directly impacting memory performance for these events. This underscores the role of the angular gyrus as a connector hub, integrating information from various regions, particularly during the process of imagination-based memory integration.

To compare the angular gyrus’ involvement in imagination-based linking with observation-based linking, where mental construction of linking events is not required, participants linked half of the stories by observing the linking events as videos. For events linked via observation, we found decreased pattern similarity in the theta band between linked events from the pre-to the post-phase. It may be tempting to speculate that viewing the linking event, reactivated only elements of the memory representations of A and B, resulting in a moderate co-activation of their memory representations, which typically leads to decreased similarity between events (Wammes et al., 2022) and more distinctive memory representations in line with previous studies (Grob et al., 2023a, 2023b; Heinbockel et al., 2022). Inhibition of the angular gyrus disrupted this decrease, potentially preventing full access to the original memory representation of A and B in all their detailed distinctiveness during the observation of the linking event, subsequently failing to induce a decrease in similarity (Ramanan et al., 2018). Following cTBS to the angular gyrus, we further observed decreased pattern similarity for non-linked events, resembling the pattern change observed in the sham group for linked events, which may highlight the role of the angular gyrus in representational separation during observation-based linking. Furthermore, for events linked via observation, we observed decreased theta coherence in the sham group. It may be speculated that decreased theta synchronization for events linked via observation indicated less working memory demands in line with prior research (Fell & Axmacher, 2011; Kawasaki et al., 2014; Sarnthein et al., 1998). Importantly, angular gyrus inhibition led to increased functional connectivity between centro-temporal and frontal regions, potentially indicating increased working memory demands. However, caution is warranted when interpreting these findings for observation-based insight because these appeared to be driven at least in part by group differences already in the pre-phase, i.e. before participants’ gained insight and before the TMS manipulation.

In conclusion, our data point to a causal involvement of the angular gyrus in (imagination-based) insight-driven memory reconfiguration. These results provide novel insights into the neural mechanisms of memory integration and bridge the traditionally separate functions attributed to the angular gyrus, namely memory and imagination. Beyond their relevance for understanding fundamental memory processes, these findings may have relevant implications for promoting the integration of fragmented memories in mental disorders, such as posttraumatic stress disorder.

Materials and Methods

Preregistration

This study was preregistered before the start of data collection at the German Clinical Trials Register (DRKS-ID: DRKS00025202; https://drks.de/search/de/trial/DRKS00025202).

Participants

Sixty-five healthy right-handed individuals (34 males, 31 females, age: M = 24.49 years, SD = 4.29 years) with normal or corrected-to-normal vision volunteered to participate in this study. Participants were screened using a standardized interview for exclusion criteria that comprised a history of neurological and psychiatric disease, medication use and substance abuse, cardiovascular, thyroid, or renal disease, body mass index below 19 or above 26 kg/m2, evidence of COVID-19 infection or exposure, and any contraindications to MRI examination or TMS. All participants gave written informed consent before participation and received a monetary compensation at the end of the experiment. The procedures were approved by the local ethics committee (Faculty of Psychology and Human Movement Science, Universität Hamburg, Hamburg, Germany, 2020_301 Grob Schwabe) and adhered to the Declaration of Helsinki. The sample size is in line with recent studies on episodic memory integration using the same task (Grob et al., 2023a, 2023b). Additionally, an a-priori power calculation using G*Power (Faul et al., 2007) indicates that a sample size of N = 54 is sufficient for detecting a medium-sized group × link effect (f = 0.25) with a power of 0.95.

We implemented a mixed-design including the within-subject factors link (linked vs. non-linked events), session (pre- vs. post-link), and mode (imagination vs. observation) as well as the between-subjects factor group (cTBS to the angular gyrus vs. sham). Participants were pseudo-randomly assigned to the cTBS group (n=34, 16 females) and the sham group (n = 31, 15 females) to achieve a comparable distribution of men and women in each group. Due to technical issues, three cTBS participants were excluded from EEG analyses.

Procedure

After obtaining participants’ written informed consent, we determined their individual motor thresholds for transcranial magnetic stimulation (TMS). Thereafter, they completed a training session of the modified narrative-insight task (NIT; Milivojevic et al., 2015), a life-like video-based task that tests the integration of initially separate events into coherent episodes (see below). During this training, participants were equipped with electroencephalography (EEG) caps and electrodes. Following the training session, participants completed the pre-phase of the NIT. After completing a German mood questionnaire (MDBF; Steyer et al., 1997), participants underwent either sham or cTBS targeting the left angular gyrus before commencing the insight-phase, with an additional stimulation session administered before the second half of the insight-phase. Crucially, this study was double-blind, ensuring that both the participant and the experimenter were unaware of the stimulation condition. Upon completion of the insight-phase, participants transitioned to another experimental room where they were given a 30-minute break, during which they completed the German mood questionnaire (MDBF) again, along with assessments of their imagination capacity (FFIS; Zabelina & Condon, 2019), trait-anxiety (STAI-T; Laux et al., 1981), and state-anxiety (STAI-S; Laux et al., 1981), depressive symptoms (BDI; Hautzinger et al., 2006), chronic stress (TICS; Schulz & Schlotz, 1999), and personality dimensions (BFI-2; Danner et al., 2016). This break was crucial for minimizing potential aftereffects of TMS that could have affected performance during the post-phase and memory tasks (Huang et al., 2005; Jannati et al., 2023). After this break, participants engaged in a self-paced free recall task, after which they returned to the EEG room to complete the post-phase of the NIT. The comparison of neural activity patterns between the pre- and post-phases allowed the analysis of insight-related changes in neural memory representations. EEG recordings were obtained during the pre-, insight-, and post-phase of the NIT. Finally, participants completed a multi-arrangements task (MAT; Kriegeskorte & Mur, 2012) and a forced-choice recognition task. In total, the experiment took about 4.5 hours per participant.

Narrative-Insight Task

To examine insight-related mnemonic integration processes, participants completed a modified version of the narrative-insight task (NIT; Milivojevic et al., 2015; Fig. 1), while their brain activity was measured using EEG. The task involved watching life-like videos from the computer game The Sims 3, representing different storylines. Each storyline consisted of events that could either be integrated (A and B) into narratives or not (A and X). Participants were unaware that each narrative had two versions. The two versions shared event A but had different events B. Event X from one version served as event B in the other version. Thus, all participants viewed the same events A, B, and X, with 37 participants linking events A and X and 28 participants linking events A and B. Counterbalancing which events were linked over participants controlled for non-specific stimulus effects and visual similarity. Throughout this manuscript, the linked events are referred to as events A and B, while the non-linked event is referred to as X. Given the role of the angular gyrus in imagination (Benoit & Schacter, 2015; Thakral et al., 2017) and recent research highlighting distinct neural underpinnings when gaining insight via imagination vs. observation, we introduced two insight modes: imagination-based linking for half of the stories and observation-based linking for the other half.

The NIT consisted of three phases (Fig. 1): the pre-phase, the insight-phase, and the post- phase. The task included ten stories featuring three videos (A, B, and X) in both the pre- and post-phase, and two videos (L, C) and two imagination instructions (I, CI) in the insight-phase. In the pre-phase, participants viewed events A, B, and X for two seconds each, with inter-trial intervals (ITIs) between 700-1300 ms (∼1000 ms). Each video was preceded by a brief title (one second) and presented 18 times in pseudorandom order. The presentation order ensured that each video appeared before the next round of presentations began, and consecutive trials did not feature the same video. Following the pre-phase, participants rated the extent to which they perceived the events as belonging together on a scale from not at all (1) to very much (4). The subsequent insight-phase comprised two parts. For five stories, participants viewed the linking video event (L) interspersed with a control video event (C), each presented for eight seconds and repeated nine times with ITIs between 700-1300 ms (on average ∼1000 ms). For the other five stories, participants observed an instruction (I) to imagine a particular linking scene, with the video titles indicating the intended person, alternating with a control instruction (CI). Each instruction was presented for eight seconds and repeated nine times with ITIs between 700-1300 ms (∼1000 ms). A total of 32 participants first linked events via imagination and later via observation, while 33 participants first linked events via observation and then via imagination. Importantly, participants were stimulated with cTBS or sham before both parts of the insight-phase to maintain the stimulation effect throughout the insight-phase. For detailed description of the TMS procedure see below.

Participants were instructed to imagine specific scenes only when specifically asked to do so; otherwise, they were told to simply relax and watch the videos or answer the rating questions. In the observation condition, the linking video (L) depicted the main characters from videos A and B interacting with each other, while the control video (C) featured an unrelated activity involving an unknown character (e.g., two women engaged in conversation). In the imagination condition, a written linking instruction (I) prompted participants to imagine the main characters from video A and B interacting with each other, while the control instruction (CI) asked them to imagine an unknown character engaged in an unrelated activity (e.g., two women talking to each other). Following the insight-phase, participants provided ratings regarding their comprehension of the link and adherence to instructions on a scale ranging from 1 (not at all) to 4 (very much). After the insight-phase, participants were taken to another room for a 30-minute break, during which they completed several questionnaires. After the break, they performed the free recall task (see below). Subsequently, in the final post-phase, events A, B, and X were again displayed for two seconds, repeated 18 times with ITIs between 700-1300 ms (on average ∼1000 ms). Each video was preceded by a one-second title. The post-phase order was pseudorandomized to minimize sequence effects. Its purpose was to examine neural representation changes for events A and B after participants learned that these were linked. Participants then rated the extent to which they perceived the events as belonging together on a scale of 1 (not at all) to 4 (very much). Participants received visual feedback in the form of highlighted selected responses when entering a rating question. In addition to the presentation of A, B, and X events in the pre- and post-phases, we included target events to ensure sustained attention throughout the experiment. These target events, accounting for 11% of pre- and post-phase trials, required participants to press a button in response to a two-second animated video of a girl on a pink scooter.

Free recall

To assess the extent to which insight into the relationship of initially unrelated events affects subsequent memory, participants performed a free recall test in which they were instructed to recall all presented events in as much detail as possible (Fig. 1). During free recall, participants were voice recorded for a maximum of 20 min. To assess the level of detailedness of the integrated episodes, audio recordings from free recall were scored according to how much detail of the different video events (A, B, and X) were recalled from day one and whether the linking events (L, and I) were named. A rating system was employed that allowed for distinct coding of details associated with each specific event (A, B, X, L, and I) and distinguished between central and peripheral details. Central details refer to elements that are crucial to the plot and directly impact the linking process. These details include significant aspects such as distinctive features of the protagonist in each event. Peripheral details encompass any observable details in the video events that are not central to the plot. For example, these could include features like the presence of a carpet in a room or the color of the curtains. Importantly, there was no difference in the number of details that could be named among the video events A, B, and X across different stories (event: F(1.35, 12.14) = 2.09, p = 0.173). However, it is worth noting that a greater number of peripheral details could be named compared to central details (detail: F(1, 9) = 83.24, p < 0.001), which was expected as there were more details visible in the video events that were unrelated to the plot and, consequently, had no direct influence on the linking process. We engaged four independent raters and instructed them to assign details only to events for which it was clear that they belonged exclusively to that event, thereby avoiding any confusion between different events. The raters further scored whether the participants named the linking events or not. The scoring process involved two raters evaluating the first half of the data, while another two raters assessed the second half. All raters were blinded to the experimental conditions. To assess inter-rater reliability, all raters rated the first five participants, and on average, these ratings were highly correlated with each other (mean correlation = 0.80, SD = 0.14). To enhance inter-subjectivity, these ratings were averaged. The details of the different event types (A, B, and X) were combined across stories to generate a comprehensive rating of event details for both imagined and observed links. The average rating of the linked events (A and B) was then calculated to represent the overall measure of linked events. The non-linked event (X) remained unchanged. The naming of the linking events were combined across stories, separately for imagination and observation. Following the free recall, participants proceeded to the post-phase of the NIT.

Multi-arrangements Task

In order to ensure that participants accurately retained the structure of the events they gained insight into, we assessed their representational structure through a multi-arrangements task (MAT; Kriegeskorte & Mur, 2012). In this task, participants were instructed to arrange representative images (A, B, and X) from each story based on their relatedness. Using a computer mouse, participants dragged and dropped the images within a circular two-dimensional arena displayed on the computer screen. This task served to assess whether participants could successfully bring the linked images (A and B) closer together than the non-linked images (A and X). Each trial was self-paced and could be concluded by the participant by selecting “Done”. In the first trial, participants had to arrange all images by similarity and were instructed to do so carefully. Subsequent trials consisted of subsets of the first trial selected based on an adaptive procedure designed to minimize uncertainty and better approximate the high-dimensional perceptual representational space. This procedure is based on an algorithm optimized to provide optimal evidence for the dissimilarity estimates (Kriegeskorte & Mur, 2012). Distances in this MA task were computed by initially computing the squared on-screen distance (Euclidian distance) between all items in the first trial to produce a roughly estimated representative dissimilarity matrix (RDM) and by iteratively updating this RDM by the weighted average of scaled trial estimates. The completion of the MA task required approximately 15 minutes. Distances for linked (A and B) and non-linked events (A and X) were averaged across stories for both imagined and observed links.

Forced-choice recognition

To further ensure participants’ accurate identification of linked and non-linked events following the NIT, a forced-choice recognition task was administered to assess participants’ comprehension. They were presented with an image of event A at the top of the computer screen and had to indicate whether the image of B or X in the bottom half of the screen was linked to A. Participants were presented with these forced-choice options for each of the stories they had seen before. After indicating for a story which event was linked to event A, they had to rate how confident they were in their answer. Confidence was rated on a scale from not at all (1) to very sure (4). This process was repeated for all ten stories. Participants were presented with the forced-choice question and the confidence rating for 5.5 seconds each, with one-second inter-stimulus intervals. Participants received visual feedback when submitting their ratings, as the selected response was highlighted. The forced-choice recognition test lasted approximately five minutes. Data from the forced-choice recognition task were pooled across stories and the percentage of correct responses was calculated, separately for imagined and observed links.

Transcranial Magnetic Stimulation

Transcranial magnetic stimulation (TMS) was applied over the left angular gyrus before participants gained insight into the relationship of initially unrelated events. We used a PowerMAG Research 100 stimulator (MAG & More GmbH, München, Germany) for stimulation, that is specifically designed for delivering repetitive transcranial magnetic stimulation (rTMS) in both clinical and research applications. Two different 70 mm figure-of-eight-shaped coils were used depending on the TMS condition: The PMD70-pCool coil with a 2T maximum field strength was used for cTBS, while the PMD70-pCool-SHAM coil, with minimal magnetic field strength, was employed for sham, providing a similar sensory experience, with stimulation pulses being scattered over the scalp.

Motor Threshold Determination

The motor threshold (MT) was assessed at the beginning of the experiment while participants were at rest, wearing an EEG cap without electrodes attached. This measurement was utilized to determine the appropriate strength of TMS required to pass through the cap. Disposable Ag/AgCL surface electromyography (EMG) electrodes were placed on the right abductor pollicis brevis (ABP) muscle, with the reference electrode on the bony landmark of the index finger and the ground electrode on the right elbow. To locate the motor hotspot, we identified the center of the head and moved 5 cm to the left and 3.5 cm forward at a 45° angle, marking it as the center point of a 9-point grid search area with each point spaced 1 cm apart from adjacent points. Starting at 40% of the maximum stimulator output (MSO), we gradually increased the intensity in 5% increments while positioning the TMS coil at a 45° angle and moving it around the search area, delivering single pulses until we identified the motor hotspot. Once the motor hotspot was located, the MT was determined at that site. It was defined as the minimum percentage of maximum stimulator output (MSO) over the left motor cortex needed to elicit motor evoked potentials (MEPs) with a peak-to-peak amplitude of 50 μV in eight out of 16 consecutive pulses.

Neuro-navigation

Before the experimental session, we obtained individual T1-weighted structural MR images using a 3T Siemens PRISMA scanner from each participant. These images were used for neuro-navigation with the PowerMag View! System (MAG & More GmbH, München, Germany). The system utilizes two infrared cameras (Polaris Spectra) to track the positions of the participant’s head and TMS coil in space. Based on the individual T1 MR images, we created 3D reconstructions of the participants’ heads, allowing us to precisely locate the left angular gyrus coordinate (MNI: −48, −67, 30) derived from previous work (Thakral et al., 2017) for TMS stimulation. The coordinates were entered as TAL coordinates. Once the TAL coordinate was entered, the coil was positioned in accordance with the template provided by the neuro-navigation system, aiming for a brain-to-target distance of less than 3 cm. This procedure ensured precise coil placement tailored to the unique anatomy of each participant, while maintaining the shortest and therefore optimal distance to the cortex.

Continuous Theta Burst Stimulation

Depending on the experimental condition, we administered continuous Theta Burst Stimulation (cTBS) using either the cTBS or the sham figure-of-eight coil at 80% of the motor threshold (MT) intensity. The experiment was conducted in a double-blind manner, where neither the participant nor the experimenter were aware of the stimulation condition (cTBS vs. sham). Previous evidence has demonstrated the inhibitory effect of cTBS on the targeted brain region under stimulation (Huang et al., 2005; Jannati et al., 2023). Following the standard TBS protocol, participants received a series of bursts comprising three magnetic pulses (pulse triplets) at a frequency of 50 Hz, with the triplets repeated at a rate of five Hz (i.e., five pulse triplets per second). Each participant received a total of 600 magnetic pulses delivered over a 40-second duration. The coil was positioned tangentially on the head and mechanically fixed in a coil holder to maintain its position. Throughout the stimulation, it was ensured via neuro-navigation that the brain-to-target distance remained below 3 cm from the left angular gyrus coordinate (MNI: −48, −67, 30).

Electric field modeling

Electric field simulations were performed in SimNIBS v4.0.1 to perform TMS simulations for the cTBS group. To assess the potential stimulation strength based on individual motor thresholds and T1 images, we segmented and meshed these MRI scans into tetrahedral head models using the SimNIBS charm pipeline. All head models were visually inspected to exclude segmentation errors. In a next step, we performed the TMS simulation at 80 % of individual motor thresholds (M = 54.82 %; SD = 12.93 %). We modeled the Magstim 70mm figure-of-eight coil placed over the left angular gyrus target coordinates (MNI: −48, −67, 30), accounting for the presence of the EEG cap during stimulation. Next, to estimate the average field strength in the region of interest (ROI), we extracted the gray matter regions and created a 10mm spherical ROI centered around the target coordinate and averaged the estimated field strength for the sphere. This approach enabled us to evaluate the potential stimulation strength and its impact on the target brain area (Pizem et al., 2022; Zhang et al., 2022).

EEG data acquisition

EEG data was recorded using a 64-electrode BioSemi ActiveTwo system (BioSemi B.V., Amsterdam, the Netherlands) following the international 10-20 system. The sampling rate was set to 1024 Hz, and a band-pass filter of 0.03-100 Hz was applied online. Additional electrodes were placed at the mastoids, above and below the orbital ridge of the right eye, and at the outer canthi of both eyes. Electrode DC offsets were kept within ± 20 μV range using the common mode sense (CMS) and driven right leg (DRL) electrodes, serving as active reference and ground, respectively.

Behavioral data analysis

To illuminate the impact of gaining insight into the relationship between initially unrelated events on subsequent memory, we subjected the number of details remembered during free recall to a linear-mixed model (LMM) implemented with the lme4 package (Bates et al., 2015) including group (cTBS/sham), mode (imagination/observation), link (link/non-link), and detail (central/peripheral) and their interactions as fixed effects, with a random intercept per participant. As a follow-up analysis, we calculated a median split on the cTBS group, based on simulations of the electric field strength and re-analyzed this data with a new group variable (sham/low stim/high stim), mode (imagination/observation), link (link/non-link), and detail (central/peripheral) and their interactions as fixed effects, with a random intercept per participant. To further examine mnemonic integration based on memory of the linking phase, naming of the linking events were entered into an LMM including group (cTBS/sham), and mode (imagination/observation), and their interactions as fixed effects. Additionally, a random intercept per participant was included to account for individual variability. As a follow-up analysis, we re-analyzed this data with a new group variable (sham/low stim/high stim), mode (imagination/observation), link (link/non-link), and detail (central/peripheral) and their interactions as fixed effects, with a random intercept per participant to estimate the effect of stimulation strength on this outcome. To verify that all participants acquired insight into the relationship between events, we analyzed the ratings for the event duplets of interest (linked events AB and non-linked events AX) from the pre- and post-phase of the NIT. This analysis was conducted using a LMM including group (cTBS/sham), mode (imagination/observation), time (pre/post), and link (link/non-link) and their interactions as fixed effects and a random intercept per participant. To confirm the retention of the representational structure of the NIT in memory, Euclidean distance estimates were extracted from the multi-arrangements task for linked (AB) and non-linked events (AX), averaged across stories, and then entered into an LMM including group (cTBS/sham), mode (imagination/observation), and link (link/non-link), and their interactions as fixed effects and a random intercept per participant. Additionally, to further ensure participants’ accurate identification of linked and non-linked events, we assessed performance in the forced-choice recognition task by calculating the proportion of correct answers. These performance measures (in %) were then entered into an LMM including group (cTBS/sham), and mode (imagination/observation), and their interactions as fixed effects and a random intercept per participant.

All analyses were performed in R version 4.0.4 and for all analyses standardized betas are reported. Prior to the analysis, the data were examined for outliers, defined as mean +/- 3 SD. For the modified NIT analysis, three outliers (two from the sham and one from the cTBS group) were identified and excluded. For the analysis of the forced-choice recognition task, two outliers (one from each group) were identified and excluded. For the analysis of the multi- arrangements task, four outliers (one from the sham and three from the cTBS group) were identified and excluded. For the free recall analysis, five outliers (two from the sham and three from the cTBS group) were identified and excluded.

EEG preprocessing

The offline analysis of EEG data from the narrative-insight task (NIT) was conducted using the FieldTrip toolbox (Oostenveld et al., 2011) and custom scripts implemented in Matlab (TheMathWorks). Pre- and post-phase trials were segmented from −2 to 3 seconds relative to stimulus onset and then re-referenced to the mean average of all scalp electrodes. The data were demeaned based on the average signal of the entire trial and de-trended. To eliminate power-line noise, a discrete Fourier-Transform filter (DFT) at 50 Hz was applied. Any electrodes that did not record or exhibited constant noise were removed (max. one per participant) and interpolated using weighted neighboring electrodes. Noisy trials were removed after visual inspection, on average 2.32 (+/- SD 1.34) of the 540 pre-phase trials and 2.63 (+/- SD 1.74) of the 540 post-phase trials. Following artifact rejection, the epochs were down-sampled to 256 Hz. Next, we performed an extended infomax independent component analysis (ICA) using the ’runica’ method with a stop criterion of weight change < 10-7 to identify and reject components associated with eye blinks and other sources of noise. In a two-step procedure, we first correlated the signals from the horizontal and vertical EOG electrodes with each independent component. Components exhibiting a correlation higher than 0.9 were immediately removed from further analysis. In a second step, the remaining components were identified through visual inspection of their time courses and corresponding brain topographies. On average, 3.32 (+/- SD 1.38) components were removed before back projecting the signals into electrode space.

Representational similarity analysis

To investigate how the brain processes insight-induced changes in the relationships between unrelated events, we conducted a representational similarity analysis (RSA) at the EEG electrode level (Heinbockel et al., 2022; Pacheco Estefan et al., 2021). RSA allows us to estimate neural activity patterns associated with specific events by measuring their correlations, thus providing insights into the underlying neural processes (Kriegeskorte et al., 2008). To measure the insight-induced representational changes, we focused on assessing the similarity of linked and non-linked events before and after gaining insight, separately for events that were linked via imagination and events that were linked via observation. We performed the RSA in the theta frequency range as prior evidence highlighted the key role of theta activity in episodic memory integration (Backus et al., 2016; Nicolás et al., 2021).

To perform this analysis, we first calculated time-frequency representations utilizing spectral decomposition using sliding Hanning windows on the preprocessed EEG data. The frequency range was set from 2 Hz to 45 Hz, with 1 Hz increments and a five-cycle window. The analysis was conducted within a time interval of −2 to 3 seconds relative to stimulus onset. For each participant, single trial power estimates were then averaged across stories and baseline corrected using absolute baseline correction with a time window of −1.8 to −1 seconds relative to stimulus onset. The time-frequency data was then appended into separate data files for the pre- and post-phase, as well as for imagined and observed stories. In a second step, we utilized the time-frequency data obtained in the theta range (4-7 Hz) to conduct RSA. These theta power values were then combined to create representational feature vectors, which consisted of the power values for four frequencies (4-7 Hz) × 41 time points (0-2 seconds) × 64 electrodes. We then calculated Pearson’s correlations to compare the power patterns across theta frequency between the time points of linked events (A with B), as well as between the time points of non-linked events (A with X) for the pre- and the post-phase separately, separately for stories linked via imagination and via observation. To ensure unbiased results, we took precautions not to correlate the same combination of stories twice, which prevented potential inflation of the data. To facilitate statistical comparisons, we applied a Fisher z- transform to the Pearson’s rho values at each time point. This yielded a global measure of similarity on each electrode site. We, thus, obtained time × time similarity maps for the linked events (A and B) and the non-linked events (A and X) in the pre- and post-phases, separately for insight gained through imagination and observation. In total, this analysis produced eight Representational Dissimilarity Matrices (RDMs) for each participant.

We performed statistical analyses on the RDMs using cluster-based permutation t-tests in the Fieldtrip toolbox (10.000 permutations; Oostenveld et al., 2011). This approach allows for testing statistical differences while controlling for multiple comparisons without spatial constraints (Maris & Oostenveld, 2007). The samples were clustered at a level of αcluster = 0.001 to allow for more refined clusters. Clusters with a corrected Monte Carlo p-value < 0.05 were considered statistically significant. The RDMs for the change from pre- to post-phase in linked events (post - pre) that were linked via imagination (vs. observation; imagination - observation) were contrasted between the cTBS (vs. the sham) groups via an independent sample cluster-based permutation t-test. Similarly, the RDMs for the change from the pre- to post-phase (post - pre) in non-linked events that were linked via imagination (vs. observation) were contrasted between the cTBS (vs. sham) groups via an independent sample cluster-based permutation t-test.

Coherence analysis

To gain a deeper understanding of the connectivity changes between linked (A and B) events before and after gaining insight, we conducted a sliding window coherence analysis in electrode-space on the same contrasts that we found in our RSA analyses. Therefore, this analysis specifically focused on the comparison between imagined and observed links, utilizing the concept of imaginary coherence (Nolte et al., 2004). Imaginary coherence quantifies the synchronization between two electrodes, accounting for phase-lag at a specific frequency and minimizing the influence of volume conduction effects. We first computed a frequency analysis focusing on the mean theta frequency (5.5 Hz; dpss-taper = 1.5 Hz). Then, we computed imaginary coherence for all possible electrode combinations. To capture the temporal dynamics, we employed a sliding window approach that spanned the duration of the video display (0-2 seconds) in 500ms windows, sliding forward in steps of 50ms. By applying this sliding window analysis, we obtained a coherence spectrum matrix of 64 (electrodes) × 64 (electrodes) for each of the 31 time windows in each participant. We calculated the average coherence matrices for both A and B, resulting in a single coherence matrix that represents the connectivity patterns of the linked events. Subsequently, we focused on the same interaction that was yielded by the previous RSA by examining the changes in coherence spectra from pre- to post-phase (post - pre) for the imagined (vs. observed) linked events. To determine the statistical significance of the observed differences between groups (sham vs. cTBS), we utilized an independent sample cluster-based permutation t-test across all 31 time windows, correcting for the multiple comparisons of channels and time windows.

In order to investigate the changes in connectivity from pre- to post-phase for the non-linked event (X), we conducted a sliding window coherence analysis following the same procedure as described above. Please note that the non-linked event was treated separately and not averaged with any other event to maintain its distinct characteristics. Similar to the linked events analysis, we obtained 31 time windows representing the change from pre- to post-phase (post - pre) for the imagined (vs. observed) non-linked event, which was the same interaction as in the previous RSA. To evaluate the statistical significance of these changes between the groups (sham vs. cTBS), we employed an independent sample cluster-based permutation t-test.

Correlational analysis

To relate the findings from the RSA and the coherence analyses to the behavioral results, we extracted the significant clusters. We then proceeded to estimate the correlation between the extracted neural cluster activity and behavioral outcomes. Specifically, we correlated neural similarity and activity with behavioral outcomes separately for the cTBS and sham groups. Subsequently, we compared these correlations to determine if they differed significantly from each other (Eid et al., 2017).

Data availability

The data generated in this study as well as all original code is available at https://doi.org/10.25592/uhhfdm.12928. For any additional information needed to reanalyze the reported data, please contact the lead researcher directly.

Author contributions

A.G.: Conceptualization, Data curation, Formal analysis, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing; H.H.: Formal analysis, Methodology, Visualization, Writing – review and editing; B.M.: Conceptualization, Writing – review and editing; C.D.: Conceptualization, Writing – review and editing; L.S.: Funding acquisition, Conceptualization, Resources, Supervision, Writing – original draft; Writing – review and editing