{'fraction_non_alphanumeric': 0.03753030007522777, 'fraction_numerical': 0.019642807389039035, 'mean_word_length': 5.640518038852914, 'pattern_counts': {'":': 0, '<': 7, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Spiritual experiences are related to engagement of a ventral frontotemporal functional brain network: Implications for prevention and treatment of behavioral and substance addictions

Background and aims Spirituality is an important component of 12-step programs for behavioral and substance addictions and has been linked to recovery processes. Understanding the neural correlates of spiritual experiences may help to promote efforts to enhance recovery processes in behavioral addictions. We recently used general linear model (GLM) analyses of functional magnetic resonance imaging data to examine neural correlates of spiritual experiences, with findings implicating cortical and subcortical brain regions. Although informative, the GLM-based approach does not provide insight into brain circuits that may underlie spiritual experiences. Methods Spatial independent component analysis (sICA) was used to identify functional brain networks specifically linked to spiritual (vs. stressful or neutral-relaxing) conditions using a previously validated guided imagery task in 27 young adults. Results Using sICA, engagement of a ventral frontotemporal network was identified that was engaged at the onset and conclusion of the spiritual condition in a manner distinct from engagement during the stress or neutral-relaxing conditions. Degree of engagement correlated with subjective reports of spirituality in the scanner (r = .71, p < .001) and an out-of-the-magnet measure of spirituality (r = .48, p < .018). Discussion and conclusion The current findings suggest a distributed functional neural network associated with spiritual experiences and provide a foundation for investigating brain mechanisms underlying the role of spirituality in recovery from behavioral addictions.


INTRODUCTION
Spiritual experiences, characterized by a felt union with a transcendent reality larger than oneself, represent a common element across various cultures and periods of human history (Eliade, 1959;James, 1902;Newberg & d'Aquili, 2008;Underhill, 1911). In an analysis of such experiences, James (1902) noted that they all share qualities of ineffability, transiency, passivity, and noesis. Examples range from a sense of oneness in nature to a transcendental state during communal worship, from a zone of intense physical activity to a feeling of inspiration and buoyancy during meditation or prayer. They occur during religious moments and in nonreligious contexts like during sporting competitions, musical performances, and experiencing of natural environments. While formal contemplative practices may induce such states, these experiences in which the ordinary sense of self is transcended may also occur without intention or forethought. They may range in intensity from gentle states of mindfulness to feelings of flow to deep mystical experiences (Davidson et al., 2003;Fredrickson, 2009;Maslow, 1964;Yaden, Haidt, Hood, Vago, & Newberg, 2017).
Spirituality has been linked to recovery processes in addictions. Twelve-step programs (e.g., Alcoholics Anonymous, Narcotics Anonymous, and Gamblers Anonymous) incorporate aspects of spirituality to varying degrees (Ferentzy, Skinner, & Antze, 2010;Galanter, 2018;Galanter, Dermatis, Post, & Sampson, 2013). Recently, we have been investigating spirituality as it relates to recovery in gambling disorders. We have found that spirituality correlates inversely with problem-gambling severity in individuals with gambling disorder at various stages of recovery (Gavriel-Fried, Moretta, & Potenza, 2019a). Furthermore, only spirituality and another positive psychology measures, recovery capital, remained significant in a model that also contained measures of anxiety, depression, and stressful life events in statistically predicting recovery from gambling disorder (Gavriel-Fried, Moretta, & Potenza, 2019b). Spirituality also appears linked to recovery from gambling disorder to a greater extent in younger as compared to older adults (Gavriel-Fried, Moretta, & Potenza, 2019c). As such, investigating spirituality has important implications for understanding recovery processes in gambling disorder, especially for younger and middle-aged adults.
Neuroimaging studies may provide valuable insight into the neural mechanisms of spiritual state experiences and their salutary effects. Previous research indicates that contemplative practices like focused awareness and open monitoring meditations are linked to structure and function of various brain regions, including the prefrontal cortex (PFC), anterior cingulate cortex (ACC), posterior cingulate cortex (PCC), insula, striatum, and amygdala, and meditative practices may promote capacities related to emotion regulation, attention control, and self-awareness (Tang, Holzel, & Posner, 2015;Vago & Zeidan, 2016). Yet, by design, this body of research does not specifically investigate the neural correlates of self-transcendent spiritual states.
Most studies of meditative and prayer practices utilize region-of-interest (ROI) or general linear model (GLM) analyses, which focus on brain regions. Alternate analytical approaches may identify functionally integrated activations or networks related to brain function, and these have been argued to possibly lead to more replicable findings (Bressler, 1995;Hermans, Henckens, Joëls, & Fernández, 2014;Xu, Calhoun, Worhunsky, et al., 2015;Xu et al., 2016). One such approach is spatial independent component analysis (sICA ;Calhoun, Adali, Pearlson, & Pekar, 2001;Zhang et al., 2015). sICA is a data-driven approach that can identify brain regions that show spatial patterns of systematic activation or deactivation and may identify overlapping activations and deactivations that cancel in GLM-based analyses Xu, Calhoun, Worhunsky, et al., 2015;Xu et al., , 2016. sICA has been proposed to generate findings consistent with brain properties such as balanced excitation/ inhibition Xu, Calhoun, Worhunsky, et al., 2015;Xu et al., , 2016. sICA has been applied to functional magnetic resonance imaging (fMRI) data to identify brain networks underlying specific processes like cognitive control, response inhibition, attention, and reward/loss processing Worhunsky, Potenza, & Rogers, 2017;Xu, Calhoun, Pearlson, & Potenza, 2014;. For these reasons, we believe sICA is well suited to identify brain networks associated with spiritual experiences. Given the importance of spirituality to recovery in behavioral addictions like gambling disorder (Gavriel-Fried et al., 2019a, 2019b, 2019c and the relevance of understanding neural correlates of treatment outcomes to refine interventions for behavioral addictions as has been proposed theoretically (Brand, Young, Laier, Wölfling, & Potenza, 2016Dong & Potenza, 2014) and demonstrated empirically in specific treatments , investigation of the neural correlates of spiritual experiences holds relevance for behavioral addictions.
Previously, we identified neural correlates of personalized spiritual experiences using a guided imagery task during fMRI . We demonstrated task validity in generating robust spiritual experiences of religious and non-religious qualities, and showed that recruitment of specific brain regions (e.g., the parietal cortex) differed in spiritual versus stressful and neutral-relaxing conditions . In this study, we apply sICA to the same data set to examine functional networks that may underlie spiritual as opposed to stressful or neutralrelaxing states.
Given perceptual, cognitive, and emotional aspects of spirituality, we hypothesized that networks including both cortical and subcortical areas would be implicated. Based on the GLM-based findings , we hypothesized that networks involving parietal function would be engaged at both the onset and ending of the conditions to a greater degree in the spiritual as compared to the stress and neutral-relaxing conditions. We further hypothesized that degree of engagement of the identified networks would correlate with subjective spiritual responses during the spiritual condition and out-of-the magnet scores on a measure of spirituality.

Participants
The study consisted of 27 community-recruited, native English-speaking participants, as described previously , as well as in "Supplementary Material." Participants completed the Spirituality Scale (Delaney, 2005), a global, holistic measure of human spirituality. Study procedures were approved by the Yale Human Investigations Committee and were in accordance with the Declaration of Helsinki. Participants provided written informed consent and received financial compensation for completing study procedures.
fMRI task, functional image acquisition, and processing The fMRI task, functional image acquisition, and image processing were performed as described previously  and detailed in the "Supplementary Material."

Independent component analysis and network selection
sICA was performed on the fMRI time series using the Group ICA of fMRI Toolbox (GIFT v1.3h; http://mialab. mrn.org/software/gift/). An estimate of 35 maximally independent components was determined using a minimum description length criterion (Li, Adali, & Calhoun, 2007). Data from all participants were concatenated into a single group and reduced through a principle component analysis. Neural network algorithms (Bell & Sejnowski, 1995) were used to extract the 35 components from the group aggregate, and extraction was repeated 50 times using ICASSO to assess stability and consistency of identified components (Himberg, Hyvarinen, & Esposito, 2004). Component time courses and spatial source maps were reconstructed and scaled to percent signal change for each session for each participant (Calhoun et al., 2001).
The spatial source maps across the six runs of each participant were averaged together, and the average spatial maps for each component were tested for significant regional loading using one-sample t-tests in SPM12 at a voxellevel family-wise error (FWE)-corrected p FWE < .01 and a cluster extent threshold of 100 voxels. Positive and negative clusters in functional networks represent inversely related component time courses, indicating simultaneous positive and negative engagement of different regions in the same network. That is, when a task-related β-weight is greater than zero, indicating positive network engagement, positive clusters exhibit positive signals (i.e., similar to the hemodynamic response function response) and negative clusters exhibit negative signals. By comparison, when a taskrelated β-weight is less than zero, indicating negative network engagement, positive clusters exhibit negative signals and negative clusters exhibit positive signals.
To identify components that were differentially associated with guided imagery conditions, multiple regression analyses were used to assess the task-relatedness of the component time courses. Guided imagery runs were subdivided into 20 s blocks of interest: the rest period immediately preceding audio (R1), the two blocks at the beginning of audio (B1 and B2), the final blocks of audio (E1 and E2), and the rest period immediately after audio (R2; Figure 1). We chose to focus on the beginning and end of the conditions based on our prior studies of emotional and motivational states (including gambling urges), given dynamic changes in brain function during these states (Potenza et al., 2003;Wexler et al., 2001). This approach also permitted examination of similar blocks across subjects, given slight variation in temporal lengths of audiotapes. Expected hemodynamic response functions were modeled for the six blocks and the single block of audio between B2 and E1 for each guided imagery run and compared to the time course of that run for each network for each participant. This process produced β-weights that represent a measure of "engagement" or "recruitment" of each network associated with each task event. The β-weights for each block from the two runs of each condition were averaged for each participant.
To identify components engaged during spirituality imagery, one-sample t-tests were performed on the β-weights of each block (B1, B2, E1, and E2) of the spirituality condition, identifying eight functional networks displaying significant spirituality-related engagement at a false discovery rate (FDR) corrected p FDR < .05. Two-way, four (block; B1-E2) by three (condition; spiritual, stressful, and neutral-relaxing) repeated-measures analyses of variance (ANOVAs) were then performed on β-weights to identify spirituality-related networks displaying a main effect of condition and/or a block-by-condition interaction, using Greenhouse-Geisser correction for sphericity violations. One functional network displayed a main effect of condition and a block-by-condition interaction (Bonferroni-corrected α .05/8 = .006).

Additional statistical procedures
Ratings of vividness following guided imagery conditions were tested using a repeated-measures ANOVA. Paired t-tests were performed on average ratings of anxiety and spirituality to compare ratings before and after guided imagery. Correlational analyses were performed to explore potential relationships between network engagement and subjective responses during the scanning sessions as well as a questionnaire assessment of spirituality on a Spirituality Scale (Delaney, 2005). Change in network engagement was computed as the difference in average β-weight during the final two blocks of each guided imagery (β E ) relative to the average β-weight from the first two blocks (β B ) of imagery. Differences in β-weights were correlated with changes in spirituality/anxiety ratings following guided imagery relative to ratings preceding imagery. The spirituality questionnaire was not completed by three participants; thus, correlations between network engagement and Spirituality Scale total scores were performed on a subset of 24 participants.

Ethics
All procedures performed in human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards.

Participants' characteristics and guided imagery ratings
Participants included 27 young adults (12 females), between the ages of 18 and 27 years, with an average of 22.4 (SD = 2.2) years, and 15.1 (SD = 1.7) years of education. Two participants provided breath samples with carbonmonoxide levels consistent with regular tobacco use, and one participant provided a urine sample positive for cannabis use. The average score on the Spirituality Scale was 105.5 (SD = 20.1).
Network engagement associated with spiritual experiences sICA identified a ventral frontotemporal network associated with spiritual imagery relative to stressful and neutral-relaxing guided imagery (Figure 2). Details regarding the regions integrated into the functional network (e.g., peak cluster locations, peak t values of spatial loading, cluster sizes, and anatomical labels) are presented in Table 1. The network was composed of temporally coherent positive signals in the middle and inferior frontal cortices, superior, middle and inferior temporal cortices, insula and frontal opercula, striatum, thalamus, brainstem, and cerebellum. Concurrent negative signals in the middle and posterior cingulate and parietal cortex were also integrated into the network.
The frontotemporal network displayed a within-subjects main effect of imagery condition (F 1.4, 37.4 = 7.69, p = .004) and a block-by-script interaction (F 6, 156 = 3.59, p = .002). Pairwise tests revealed greater engagement during spiritual imagery relative to stressful (p = .016) and neutral-relaxing (p = .004) imagery conditions. The frontotemporal network was positively engaged during the initial block (B1) of all imagery conditions, which demonstrated positive engagement during the ending blocks of spiritual imagery (B2-E2), and was negatively engaged or not engaged during the ending blocks of stressful and neutral-relaxing imagery. Changes in network engagement at the end relative to the beginning of the spiritual imagery condition (i.e., βE-βB) were correlated with changes in ratings of spirituality (r = .71, p < .001) following spiritual imagery and with total scores on the Spirituality Scale (r = .48, p = .018; Figure 3).

DISCUSSION
This is the first study to examine distributed functional brain networks involved in spiritual experiences. Spiritual experiences engaged a ventral frontotemporal network that includes middle and inferior frontal cortices, superior, middle and inferior temporal cortices, frontal opercula, and anterior insula, a finding that builds upon previous neuroimaging studies of self-transcendence, religiosity, and Figure 1. (a) Schematic of timeline and epochs used in fMRI data analysis. The R1 and R2 periods constitute the 20-s periods immediately prior to and immediately following the listening period when the audio file was being played and participants actively engaged in guided imagery processes. The B1 and B2 periods are the first and second 20-s periods during the beginning of the audio file listening. The E1 and E2 periods are the next-to-last and last 20-s periods at the end of the audio file listening. (b) Subjective responses on vividness of story, spiritual connection, and anxiety across the three conditions. Participants were questioned prior to listening to their individualized script ("Baseline") and following the script ("Post-imagery"). ***p < .001 spirituality. An associated network included subcortical regions such as the striatum, thalamus, and brainstem. Negative signals in the middle and posterior cingulate and parietal cortices were also integrated. The degree of engagement of this network at the beginning and end distinguished the spiritual from the stressful and neutral-relaxing conditions. The degree of network engagement during the spiritual condition was associated with subjective changes in the strengths of feelings of spiritual connectedness during the spiritual condition and to an out-of-magnet measurement of spirituality more broadly. As such, our hypotheses were supported. How this network may be understood in relationship to prior findings is described below and implications for behavioral addictions and their prevention and treatment are considered below.

Ventral frontotemporal network
The observed ventral frontotemporal network appears to correspond to a well-documented ventral attention network (VAN) composed of the temporoparietal junction and ventral frontal cortex and which responds to unexpected but  Note. Cluster details of the frontotemporal network at voxel-level p FWE < .01, k > 100. +/−: positive/negative signal integration; BA: Brodmann area; k: cluster size in contiguous voxels; x, y, z: cluster peak location in MNI coordinates; t: cluster peak t-score; R/L: right/left. behaviorally salient stimuli (Corbetta, Patel, & Shulman, 2008;Corbetta & Shulman, 2002;Macaluso, 2010;Vossel, Geng, & Fink, 2014). While multiple studies report a right-hemisphere-dominant VAN, some attentional tasks including those that involve particularly rare stimuli have implicated a bilateral network (Geng & Mangun, 2011;Serences et al., 2005;Vossel, Weidner, Thiel, & Fink, 2009). In a task that involves responding to a visual stimulus after a predictive spatial cue is provided, a separate dorsal attention network (DAN) is implicated when the stimulus location is precued, while the VAN becomes active when stimuli appear outside the previously cued area of focus (Kincade, Abrams, Astafiev, Shulman, & Corbetta, 2005;Vossel et al., 2009). These two primary attention networks may interact in a dynamic fashion (Vago & Zeidan, 2016). The DAN in a top-down fashion may inhibit pathways between the VAN and incoming sensory and perceptual information, which may effectively filter out unintended information. When an unexpected but personally relevant stimulus comes into one's field of perception, this inhibitory filter may be released, giving way to VAN-mediated bottom-up, stimulus-driven, and involuntary attentional processes (Corbetta & Shulman, 2002). That the observed pattern of engagement of the VAN, and not the DAN, is associated with the spiritual condition and not the other two conditions suggests that spiritual experiences, which James (1890) once described as "intense, voluminous, or sudden" may involve attentional processes that occur outside of Speculatively, this association may reflect the often unanticipated but personally meaningful nature of spiritual experiences. By extension and also speculatively, this finding suggests that neural networks underlying attentional processes may be important in the role of spirituality in recovery processes in behavioral addictions like gambling and gaming disorders. As both gambling and gaming disorders have been linked to attention-deficit hyperactivity disorder (Karaka, Canan, Saleh, & Potenza, 2017;Potenza et al., 2019), including with respect to the disorder persistence (Bruneau et al., 2016;Fatseas et al., 2016), the findings raise the possibility that impairments in attentional processing in behavioral addictions may extent to spiritual processes and may be targeted through interventions.
In addition to involvement in attentional processes, a ventral frontotemporal network has been implicated in processing representations of others and relational bondedness. The superior temporal sulcus has been linked to social perception (Lahnakoski et al., 2012), and co-activation of frontal and temporal regions has been implicated in understanding emotions and intentions of others, including those of G-d (Kapogiannis, Barbey, Su, Krueger, & Grafman, 2009;Lissek et al., 2008;Wang, Lee, Sigman, & Dapretto, 2007). Activations of both frontal and temporal lobes including middle and superior temporal gyri as well as the superior frontal gyrus and anterior insula have been related to bonding with intimate others, such as mothers and romantic partners (Decety & Jackson, 2004;Laurita, Hazan, & Spreng, 2017;Noriuchi, Kikuchi, & Senoo, 2008;Vanderwal, Hunyadi, Grupe, Connors, & Schultz, 2008). Taken together, the frontotemporal network engagement in the context of the literature raises the possibility, albeit speculatively, that a feeling of relational intimacy may accompany spiritual states. This interpretation also has implications for behavioral addictions as relational attachments have been implicated in addiction processes, including with respect to the transgenerational transmission of addictions (Alvarez-Monjaras, Mayes, Potenza, & Rutherford, 2019). As such, interventions examining spirituality in recovery from behavioral addictions should incorporate measures of social relatedness and attachment.
Previous research of large-scale neural networks have also linked frontotemporal activations to episodic memory retrieval (Barredo, Öztekin, & Badre, 2013;Simons & Spiers, 2003). It is therefore plausible that memories related to spiritual experiences are stronger and more salient, an interpretation consistent with the profound personal significance previously attributed to intense spiritual experiences (Maslow, 1962(Maslow, , 1964Waldron, 1998). This interpretation would suggest that the enhancement of positive emotional experiences in the treatment of behavioral addictions warrants consideration. In addition, a left-hemispheric frontotemporal network has been observed to be involved in producing and comprehending spoken language (Hagoort & Indefrey, 2014;Rodd, Longe, Randall, & Tyler, 2010). While it is possible that language might be processed differently during spiritual states, the present data identify a bilateral frontotemporal network activation, arguably lending less support to this interpretation.

Negative signals in the PCC and inferior parietal lobule (IPL)
The negative signal in the PCC suggests a possible reduction in self-related mental processes (Brewer & Garrison, 2014;Brewer et al., 2011;Panda et al., 2016). As a component of the default mode network and a hub with dense anatomical connections to multiple regions in the brain, the PCC has been implicated in multiple processes involving a selfreferential component. For example, the PCC is involved in activities ranging from self-evaluations to autobiographical memories to dealing with moral dilemmas involving oneself (Leech, Braga, & Sharp, 2012;Morey et al., 2012;Sporns, 2013), and PCC deactivation has been correlated with tasks that involve focused, present-oriented attention beyond one's sense of self (Brewer, Garrison, & Whitfield-Gabrieli, 2013;Garrison et al., 2013). In behavioral and substance addictions, PCC activation has been linked to reward-related and craving processes (Kosten et al., 2006;Yao et al., 2017). For these reasons, the PCC has been proposed as a neural target for addictions, especially for meditation and mindfulness-based interventions targeting stress reduction (Brewer & Garrison, 2014;Kral et al., 2019). As we have found that measures of spirituality relate inversely to stress-related activations in largely subcortical regions including the striatum and thalamus , the findings taken together suggest important roles for both cortical and subcortical processes in the beneficial effects of spirituality as well as potential neural targets in spirituality-related interventions for addictions.
Consistent with our hypothesis, the IPL also demonstrated a concurrent negative signal. The right lateralized IPL has been implicated in body representation in three-dimensional space (Berlucchi & Aglioti, 1997;Maguire et al., 1998;Vogeley & Fink, 2003). This result is consistent with an expanded sense of self in space that has been reported in previous phenomenological accounts of spiritual states (Newberg & Newberg, 2006), as well as the scripts generated for the current experiment. Reduced activity in the IPL has also been linked to a reduction in perceived distinction between self and other as well as a reduced sense of agency (Uddin, Molnar-Szakacs, Zaidel, & Iacoboni, 2006). Thus, these regional deactivations are consistent with the notion that perceptual changes occur. More specifically, identification with one's physical self-sense may become more relaxed while, simultaneously, perceived boundaries between self and other may become more diffuse. Taken together, these changes may give rise to a less bounded and more expanded sense of self. Supporting this interpretation, reductions in activity in the inferior parietal cortex (and PCC) have been associated with self-transcendence, religious and spiritual beliefs, and religious and spiritual commitments (Brewer et al., 2011;Crescentini, Aglioti, Fabbro, & Urgesi, 2014;Urgesi, Aglioti, Skrap, & Fabbro, 2010).

Subcortical involvement
The identified frontotemporal network also contained concurrent positive signals in subcortical regions that include the putamen and caudate nucleus, illuminating possible neural substrates of the concomitant emotional quality associated with spiritual states. These regions have been implicated in the processing of positive and rewarding emotions like love and bliss, which often accompany spiritual experiences (Bartels & Zeki, 2000;Beauregard & Paquette, 2006;Damasio et al., 2000;Langeslag, van der Veen, & Röder, 2014). Because the neutral-relaxing condition also consisted of positively valenced experiences (e.g., listening to pleasant music or relaxing with family members), the relative involvement of these regions suggests that the positive emotion associated with spiritual experience is distinct from and may be more pronounced than simply recalling a positive memory. Saver and Rabin (1997) previously proposed that subcortical limbic regions are responsible for salient aspects of spiritual experiences including a sense of unity and feelings of ecstasy. That dorsal striatal and not ventral striatal regions were implicated in this study suggests that people may seek these experiences not because of an anticipated reward but rather from an intrinsically motivated state, possibly reinforced through previous positive experiences (Balleine, Delgado, & Hikosaka, 2007;Delgado, 2007). Given the role of the dorsal striatum in other processes (e.g., habits) and overlap with regions of the striatum implicated in reward processing, alternate interpretations exist and warrant direct examination in future studies (Brewer & Potenza, 2008;Everitt & Robbins, 2005;Yin & Knowlton, 2006). Thus, while spiritual states may represent a form of positive reinforcement that may increase a person's likelihood of continuing to seek more (Otto, 1946;Underhill, 1911), given the wide range of processes in which the striatum participates, this and other possibilities should be studied directly in future experiments.

Spiritual experiences and volitional contemplative practices
Reduced PCC and IPL activity represents a common feature across a variety of meditative practices and spiritual experiences (Brewer et al., 2011;Fox et al., 2016;Urgesi et al., 2010), reflective of reduced conceptual, and self-related processing. Dorsal striatal activity, however, is not typically implicated in volitional contemplative practices (Farb et al., 2007;Fox et al., 2016;Garrison, Zeffiro, Scheinost, Constable, & Brewer, 2015), potentially reflecting an absence of concomitant positive emotion, although some exceptions do exist with certain practices (Hagerty et al., 2013;Tang et al., 2009). Furthermore, in contrast to our findings on spiritual states, volitional practices like focused attention and open-monitoring meditations activate neural regions involved in enhancing attentional control, particularly the ACC Hölzel et al., 2007;Tang, Tang, & Posner, 2016). The ACC monitors multiple channels of information and enables voluntary executive control of attention (Van Veen & Carter, 2002). Other areas related to executive attention like the dorsolateral and ventrolateral prefrontal cortices have also been implicated in meditative practices Hölzel et al., 2013). However, in a spiritual state, by contrast, there appear to be fewer correlates of volitional attention.

Implications for the prevention and treatment of addictions
Mindfulness-based practices often include spiritual elements and have empirical support in the treatment of addictions (Brewer, Bowen, Smith, Marlatt, & Potenza, 2010;Brewer et al., 2009;Witkiewitz et al., 2014). Similarly, mindfulness has been considered in the treatment of gambling disorder (Chen, Jindani, Perry, & Turner, 2014;de Lisle, Dowling, & Allen, 2012;Shonin, Van Gordon, & Griffiths, 2014), Internet gaming disorder (Zhang, Yao, Potenza, Xia, Lan, Liu, et al., 2016a;, and compulsive sexual behaviors (Blycker & Potenza, 2018). As such, the current findings suggests mechanisms by which aspects of such interventions may operate, and studies directly examining how interventions that include mindfulness-related elements should consider examining directly changes in spiritual processes and frontotemporal-network function. Such work could be extended to understanding how 12-step programs promote recovery (Ferentzy et al., 2010;Galanter, 2018;Galanter et al., 2013) and how mindfulness-based practice may promote health (Lin et al., 2019). As spirituality may mitigate against stress-related illnesses (Koenig, 2012) and stress has been linked to drug and behavioral addictions Sinha, 2008), including at neural levels (Potenza et al., 2012), direct examination of how spirituality may prevent substance and behavioral addictions and operate through the reduction of stress should be conducted. Preliminary data in this area suggest that spirituality may buffer against neural responses to stress, particularly at subcortical levels . However, the current findings suggest that spirituality may also promote positive health psychological benefits that operate through VAN-related neural processes. The extent to which such networks may also influence other positive health psychology domains (such as those operationalized as component of recovery capital; Gavriel-Fried et al., 2019a, 2019b, 2019c and promote resiliency against and recovery from behavioral addictions warrants direct examination in future studies.

Limitations, future directions, and conclusions
Limitations should be noted. Although participants represented a variety of religious backgrounds, the study drew from a moderately small sample of young adults, which may limit the ability to detect more subtle between-condition differences and increase the possibility of false-positive findings. However, the sample size is standard for withinsubject designs, and the current findings survived clusterlevel corrections, although these findings warrant replication before drawing definitive conclusions. Further investigations that employ larger samples representing a range of cultures and age groups would also allow for broader generalization and for analyses of specific categories of spiritual experiences. Because the neutral-relaxing condition also consisted of positively valenced experiences (e.g., listening to pleasant music or relaxing with family members), the relative involvement of these regions suggests that the positive emotion associated with spiritual experience is distinct from and may be more pronounced than simply recalling a positive memory (e.g., communal worship experiences vs. connections with nature). Due to the relatively sparse literature on both spiritual experience and ventral frontotemporal networks, reverse inferences that were discussed should be considered cautiously.
Despite limitations, the current findings suggest a pattern of neural networks associated with spiritual experiences and have potentially important implications for understanding recovery from behavioral addictions. Correlated with both perceived subjective changes in spirituality and overall spirituality, the identified functional network is spatially consistent with networks implicated in a range of processes, chiefly: (a) involuntary reorientation of attention, (b) a more interconnected perception of oneself, and (c) enhanced positive affect consistent with attachment or bonding. Given that spirituality has been linked to better mental health across contexts and conditions, future studies should examine the extent to which the findings may relate to improved mental health, especially within clinical settings relating to behavioral addictions.
Funding sources: This work was supported in part by the National Institutes of Health (R01DA039136), the Connecticut State Department of Mental Health and Addiction Services, the Connecticut Mental Health Center, the Connecticut Council on Problem Gambling, a Center of Excellence in Gambling Research Award from the National Center for Responsible Gaming, and the Peter Boris Centre for Addictions Research.
Authors' contribution: LM and MNP conceived original study idea. RS designed guided imagery paradigm and LM, MNP, RS, CHMC, and IMB contributed to study design. IMB and RS coordinated data collection. JX preprocessed data and performed ICA. CHMC and PDW interpreted results and wrote the manuscript. All authors revised the manuscript critically and approved the final version. CHMC, PDW, LM, and MNP contributed equally to the generation of the manuscript.

Conflict of interest:
None of the authors have any relevant financial disclosures. Dr. MNP has received financial support or compensation for the following: he has consulted for RiverMend Health, Game Day Data, the Addiction Policy Forum, and Opiant Pharmaceuticals; has received research support from Mohegan Sun Casino and the National Center for Responsible Gaming; and has consulted for gambling and legal entities on issues related to addictive disorders.