Autonomic nervous system modulation during self‑induced non‑ordinary states of consciousness

Different non-ordinary state of consciousness practices utilize voluntary control of breathing to induce spe- cific increases or decreases in HRV by manipulating breathing rate and variability. For example, practices like Samatha and Vipassana meditation or slow breathing techniques aim to increase HRV, while practices like holotropic breathwork aim to decrease it 28–32 . However, not all non-ordinary state of consciousness practices involve the voluntary control of breathing. In some forms of meditation and hypnosis, increased vagal control has been observed without specific control over breathing 22,33,34 . Nevertheless, the literature shows some divergent effects depending on the specific tech- nique used. For instance, certain types of meditation induce parasympathetic withdrawal (i.e., decreased HRV), while others produce the opposite effect (i.e., increased HRV)35 . Previous studies have also shown that hypnosis induces a systematic increase in parasympathetic control 22,36,37 . These different physiological findings related to autonomic nervous system modulation in the context of non-ordinary state of consciousness may suggest relevant clinical opportunities. In particular, drug-free self-induced and volitional non-ordinary states of consciousness are being increas- ingly tested and integrated into various clinical settings. These practices are employed in the treatment of mood disorders, pain reduction, and lead to improvement of quality of life 38,39 . Among these non-ordinary states of consciousness practices, self-induced trance states, which have historically received less attention, are believed to possess significant clinical potential. Trance states are often characterized along two dimensions: psychophysi- ological and cultural 40 . Various types of trance, such as hypnotic, shamanic, and possession trance, have been documented, and while they exhibit differences, they also share common phenomenological and psychological features 41,42 . Self-induced cognitive trance (SICT) has recently been introduced as a volitional non-ordinary state of consciousness. It is adapted from traditional Mongolian shamanic trance, abstracted from any ritual or spiritual expression1,42 . SICT and shamanic trance are characterized by a lucid but narrowed awareness of external surroundings, a hyper-focused immersive experience of flow, expanded inner imagery, modified somatosensory processing, altered sense of self, subjective changes in time and space perception, modified body awareness, and cognitive and emotional states, along with disengagement from the sensory environment 1,43–47 . Unlike shamanic trance, SICT does not require the presence of a shaman to access the trance state. It is learned through a standard- ized training, initially using specific sound loops and later individualized movements or vocalizations that allow for deliberate and controlled induction of non-ordinary state of consciousness 41,42 . However, the neurophysi- ological and autonomic underpinnings of SICT remains unknown. Autonomic regulation, along with several other physiological processes, plays a significant role in changes of consciousness states. Therefore, we hypothesized that SICT is associated with modulations of autonomic regulation. To test this hypothesis, we monitored HRV and breathing during SICT. The aim of this study is to investigate cardio-respiratory changes during SICT, with a specific focus on identifying the features among cardio-respiratory metrics that characterize the SICT state in comparison to a baseline resting state and an imagination control task condition. Methods Participants and training We recruited a total of 27 adult participants who were native French speakers and had expertise in SICT. The participants had an average practice time of 32.6 months (± 50; min = 9; max = 216). It was a requirement for par- ticipants to be capable of remaining motionless during their trance. Data analysis was conducted on 26 subjects (with 23 females, mean age of 45 yrs, ± 13.33; min = 24; max = 72). One subject was excluded from the analysis due to excessive artifacts in the recorded data, such as movements and poor recording quality. All participants received standardized training in SICT, which involved a system of techniques based on a sound loop developed by Corine Sombrun, an ethno-musician, and the TranceScience Research Institute (https:// trancescience.org). This training program, inherited from shamanic Mongolian traditional practice, utilized the sound loop as a means for participants to enter into SICT. The objective of the training was to enable participants to voluntarily induce trance without the continuous use of sound. Following a training period consisting of two 3 Vol.:(0123456789) Scientific Reports | (2023) 13:15811 | https://doi.org/10.1038/s41598-023-42393-7 www.nature.com/scientificreports/ weekends of practice with this standardized program, most participants were generally able to induce trance voluntarily. Participants were then given the opportunity to continue practicing autonomously at home. Prior to their participation in the study, participants were fully informed about the study’s objectives and provided written consent. No incentives were offered and informed consent was obtained from all participants. The study received approval from the Ethics Committee of the Faculty of Medicine at the University of Liege and all experiments were performed in accordance with relevant guidelines and regulations. Procedure Prior to the experiment, demographic data including age and gender were collected from each participant. The experimental session consisted of five conditions: ordinary conscious resting state (’Rest’), ordinary conscious state with auditory stimulations (’Auditory’), imagination task of a previous intense trance without entering a trance state (’Imag’), SICT, and SICT with the same auditory stimulations as the ordinary conscious state auditory condition (’Auditory-SICT’). The first three conditions were counterbalanced among participants, while the last two conditions were always conducted after the first three to avoid potential after-effects of SICT. The order of the last two conditions was consistent for all participants, starting with SICT without external stimulation, as it is commonly practiced in their daily life. During each condition, participants were instructed to keep their eyes closed. In the Rest and Auditory conditions, they were asked to let their thoughts freely come and go. In the Imag condition, participants were instructed to imagine a previous intense SICT experience without actually entering a trance state. In the SICT and Auditory-SICT conditions, participants were instructed to induce and maintain SICT. The specific techniques used to induce SICT varied based on the participants’ preferences and habits, involving body movements and/or vocalizations and lasting between 2 and 10 min. Once participants reached the trance state, they were instructed to remain motionless until the end of the recording. However, if the trance started fading away, participants were allowed to reinduce the trance and the recording was extended accordingly. Each condition had a duration of approximately 12 min. Following the experimental procedure, participants completed a self-report to indicate whether they reached a trance state and rate the intensity of their experience on a Likert scale ranging from 0 (no trance) to 10 (the most intense trance ever experienced). It is important to note that this work excluded the two ’Auditory’ condi- tions, as they are the focus of another related project. Therefore, the present study focused on three conditions: baseline resting state (Rest), imagination (Imag), and SICT. Electrophysiological recordings During the experiment, we collected data on various body parameters, including electrocardiogram (ECG) and respiration using the EGI polygraph input box. The ECG data was obtained by placing two electrodes above and below the heart. For respiration monitoring, two belts were used, one placed around the chest and another around the belly. These measurements allowed us to capture physiological changes during the different experimental conditions and gain insights into the participants’ physiological responses during the study. Data preprocessing and feature extraction Cardiac analysis For the cardiac analysis, the ECG time series data for each condition (Rest, Imag, and SICT) underwent pre- processing steps. The raw recordings were first cleaned by applying a high-pass Butterworth filter with a cut-off frequency of 0.5 Hz (order = 5). Powerline filtering was then performed to remove any interference from the electrical power source. After the cleaning and preprocessing steps, 11 min of data (the longest duration available across all partici- pants after removing the induction time) were selected for further analysis. R-peak detection was performed to locate the R peaks in the ECG signal, and the RR intervals (the time intervals between successive R peaks) were computed. Manual inspection was carried out to identify and correct any abnormal detections, such as aberrant RR intervals. HRV features were extracted using Neurokit2, a Python package that provides advanced biosignal processing capabilities 48 . The following time domain features of HRV were computed: heart rate, mean of RR intervals (RR mean), standard deviation of RR intervals (SDNN), root mean square of the standard deviation of RR intervals (RMSSD), and the probability of RR intervals greater than 50 ms (pNN50). Frequency domain analysis was performed by applying the Welch method to compute the power spectral density of the RR intervals. The low frequency band (LF, 0.04–0.15 Hz) and the high frequency band (HF, 0.15–0.4 Hz) were extracted. Normal- ized powers (nLF, nHF) were obtained by dividing the power in a given frequency band by the total power. The logarithmic transformation of HF (lnHF) was also computed, along with the LF/HF ratio. Nonlinear metrics of HRV were assessed using the approximate entropy measure (ApEn), which quantifies the complexity of the HRV time series. These preprocessing steps and feature extraction procedures allowed for the characterization of various car- diac parameters and HRV indices during the different experimental conditions (Rest, Imag, and SICT). Table 1 provides additional details about the extracted features and their interpretation. HRV correction In the HRV literature, several authors have highlighted the influence of heart rate and respiration rate on HRV metrics 30,49–52 . To account for the correlation between heart rate and HRV, we performed a correction on the HRV frequency domain metrics. Specifically, we multiplied each frequency metric by the mean RR interval for all 4 Vol:.(1234567890) Scientific Reports | (2023) 13:15811 | https://doi.org/10.1038/s41598-023-42393-7 www.nature.com/scientificreports/ conditions, as described by Sacha (2014). Following this correction, no significant dependence (p > 0.05) between HRV and heart rate, as well as between HRV and respiration rate, was found using correlation Pearson test. Phasic and tonic HRV According to the vagal reservoir theory (Laborde et al., 2018b), HRV should be assessed not only at rest but also in response to a task or during a particular state (such as effort or stress). The ability to modulate HRV during a task compared to its resting state is referred to as phasic HRV. To compute phasic HRV, we subtracted the values of SICT minus Rest and Imag minus Rest conditions as follows: Subsequently, Spearman correlations between tonic HRV (Rest condition) and phasic HRV across participants were computed specifically for the HF HRV metric. Respiration analysis Similar to the cardiac signal analysis, we applied a similar procedure to analyze the respiratory signal. Respiratory time series were extracted for 11 min for each condition and each participant. The respiratory signal underwent preprocessing, including artifact cleaning (linear detrending followed by a fifth-order 2 Hz low-pass IIR But- terworth filter), and preparation for breath peak detection. Using Neurokit2 48 , we extracted features from the preprocessed respiratory signal. The mean respiratory rate and amplitude were computed, along with the inspiratory to expiratory time ratio or Phase Duration Ratio. To compute the BB interval (interval between breath peaks), we located the B peak corresponding to exhalation onsets. Similar to the RR interval for the cardiac signal, every BB time series was manually checked and corrected for abnormal detection, such as aberrant BB intervals. For the RRV analysis, we computed time domain features including the standard deviation of BB interval (SDBB) and the root mean square of the standard deviation of BB interval (RMSSD). These metrics provide insights into the variability of breath-to-breath intervals and reflect the dynamic nature of the respiratory system. Machine learning Cardiac and respiratory signal classifications were performed using a linear discriminant analysis model within a leave-one-subject-out cross-validation scheme. The goal was to assess the discriminative changes in HRV and RRV between the Rest, Imagination, and SICT conditions. In this approach, the model was trained on n − 1 subjects (25 individuals) and evaluated on the remaining sub- ject. This process was repeated until each subject had been left out, and performance metrics were averaged across these repetitions. Three distinct classifications were conducted: Rest vs. SICT, Rest vs. Imagination, and SICT vs. Imagination. The purpose of these classifications was to characterize the differences between each condition. Each classification was performed using a single feature at a time, utilizing the 27 features (16 cardiac metrics and 11 respiratory metrics) extracted from the data. The performance of the classification was assessed based on the decoding accuracy (DA), which represents the percentage of correctly classified observations in the test set. To determine the statistical significance of the obtained decoding accuracy, permutation tests were applied. This involved generating a null distribution of decoding accuracies by randomly shuffling class labels and run- ning multiple instances of the classification (n = 1000). By comparing the observed decoding accuracy to the null distribution, a statistical threshold was derived to assess significance. Maximum statistics were applied to control for multiple comparisons across all the metrics, following the methods described by 53–55 . HRVphasic (Imag) = HRV(Imag) − HRV(Rest) HRVphasic (SICT) = HRV(SICT) − HRV(Rest) Table 1. Heart Rate Variability (HRV) time and frequency domain metric, definition, and physiological mechanisms (adapted from Berntson et al., 1997; Laborde et al., 2017; Malik, 1996; Shaffer & Ginsberg, 2017). Name Definition Physiological mechanism Time domain SDNN Standard deviation of all RR intervals Cyclic components responsible for heart rate variability RMSSD Root mean square of successive differences Cardiac vagal activity pNN50 Percentage of successive normal sinus RR intervals more than 50 ms Cardiac vagal activity ApEn Approximate Entropy of RR intervals Frequency domain LF Low frequencies (0.04–0.15 Hz) Mix of cardiac sympathetic and vagal activ- ity, as well as baroreflex activity HF High-frequencies (0.15–0.4 Hz) Cardiac vagal activity for ms2 only LF/HF Low frequencies/high-frequencies ratio Unclear, suggested mix of cardiac sympa- thetic and vagal activity LFn Low frequencies divided by the total power HFn High frequencies divided by the total power LnHF Logarithmic transformed HF Cardiac vagal activity 5 Vol.:(0123456789) Scientific Reports | (2023) 13:15811 | https://doi.org/10.1038/s41598-023-42393-7 www.nature.com/scientificreports/ Results Self‑induced cognitive trance behaviors All participants reported being able to reach a trance state in the SICT condition, with a reported mean intensity of 6.69 (sd = 1.89; min = 3; max = 10). The mean time of induction was 3.03 min (min = 1; max = 5). Time and frequency domain of HRV We observed an increase in heart rate during SICT (mean = 81.07 ± 12.53, standard deviation), compared to Rest (69.25 ± 8.67; DA = 71.9%) and between SICT and Imag (71.2 ± 9.88; DA = 68%) while no difference was found between Rest and Imag (DA = 56%). RR mean interval was lower in SICT (764.82 ± 121.8) compared to Rest (884.16 ± 120.64; DA = 69.9%) and also compared to Imag (862.27 ± 128.38; DA = 71.9%) while no significant difference was found between Rest and Imag (DA = 52%) (Fig. 1). HRV time domain SDNN was higher in SICT (67.45 ± 28.08) compared to Rest (51.57 ± 17.37; DA = 69.9%) and also compared to Imag (47.43 ± 12.75; DA = 76%) while no difference was found between Rest and Imag (DA = 50%). No difference was found for RMSSD between SICT (34.53 ± 25.30) and Rest (37.40 ± 19.57; DA = 50%) nor between SICT and Imag (34.63 ± 18.06) (DA = 52%), nor between Rest and Imag (DA = 52%). No difference was found for pNN50 between SICT (13.45 ± 17.68) and Rest (16.59 ± 17.06) (DA = 41.9%), nor between SICT and Imag (14.94 ± 18.23) (DA = 47.9%), nor between Rest and Imag (DA = 52%). HRV frequency domain We found a decrease for HF in absolute power (ms2 ) between SICT (0.0020 ± 0.0020) and Rest (0.0079 ± 0.0020; DA = 70%) but no statistical difference between SICT and Imag (0.0070 ± 0.0069; DA = 64%) nor between Rest and Imag (DA = 50%). We reported a decrease for HFn between SICT (0.13 ± 0.11) and Rest (0.29 ± 0.13; DA = 76%) and between SICT and Imag (0.26 ± 0.17; DA = 69%) while no difference was found between Rest and Imag Figure 1. HRV time and frequency domain during Rest (green), Imagination (Imag, yellow) and self-induced cognitive trance (SICT, purple). SDNN standard deviation of RR interval, RMSSD the root mean square of standard deviation of RR interval, pNN50 the probability of RR interval greater than 50 ms, LF low frequency, HF high frequency, LFn low Frequency normalized, HFn high Frequency normalized, lnHF logarithmic transformed HF, ApEn approximate entropy (corrected for multiple comparisons, *p = 0.05; **p = 0.01). 6 Vol:.(1234567890) Scientific Reports | (2023) 13:15811 | https://doi.org/10.1038/s41598-023-42393-7 www.nature.com/scientificreports/ (DA = 54%). A decrease for lnHF between SICT (−6.80 ± 1.28) and Rest (−5.28 ± 1.32; DA = 71.9%) and between SICT and Imag (−5.45 ± 1.05; DA = 69.9%) were found while no difference was found between between Rest and Imag (DA = 59%). LH/HF ratio was found to be increased during SICT (5.59 ± 5.18) compared to Rest (2.39 ± 4.68; DA = 71.9%) and compared to Imag (3.16; ± 5.65; DA = 64%), as well as an increase in Rest compared to Imag (DA = 56%). We found a drop in ApEn when comparing SICT (1.04 ± 0.21) to Rest (1.25 ± 0.15; DA = 73.9%) or to Imag (1.28 ± 0.16; DA = 68%) while no difference in ApEn was found between Rest and Imag (DA = 50%). We did not find statistical difference for LF and LFn across conditions (Fig. 2, Supplementary Information). Respiration rate, amplitude, phase and variability Thoracic respiration rates showed a trend to decrease during SICT but no significant difference was found across all conditions, SICT (14.96 ± 6.54) compared to Rest (15.41 ± 4.38; DA = 58%), SICT compared to Imag (15.63 ± 4.80; DA = 62%) and Imag compared to Rest (DA = 46%) (Figs. 3 and 4). Respiratory amplitude increased during SICT (374.89 ± 225.53) compared to Rest (209.86 ± 134.01; DA = 68%) and between SICT compared to Imag (226.38 ± 136.33; DA = 68%) while no difference was found between Rest and Imag (DA = 54%). The phase ratio duration between inspiration and expiration was increased between SICT (1.04 ± 0.36) compared to Rest (0.74 ± 0.2; DA = 70%) and between SICT and Imag (0.83 ± 0.25; DA = 64%), while no difference was found between Rest and Imag (DA = 58%). Respiratory rate variability changed across conditions. An increase in standard deviation (SD) was found between SICT (2.60 × 10 3 ± 1.48 × 10 3 ) compared to Rest (1.32 × 10 3 ± 0.73 × 10 3 ; DA = 68%) and between Imag (1.2 × 10 3 ± 0.65 × 10 3 ) and SICT (DA = 66%), while no difference was found between Rest and Imag (DA = 48%). The root mean square of standard deviation (RMSDD) was higher during SICT (3.01 × 103 ± 1.77 × 103 ) compared to Rest (1.44 × 10 3 ± 0.92 × 10 3 ; DA = 70%) but also compared to Imag (1.35 × 10 3 ± 0.81 × 10 3 ; DA = 66%), while no difference was found between Rest and Imag (DA = 48%). Phasic and tonic vagal cardiac control We found that phasic HF during SICT decreased compared to tonic HF while in the Imag condition no difference were observed (Fig. 5). We also found significant spearman correlations across participants between phasic and tonic during SICT for HF (r = −0.92); while the same correlation between phasic and tonic HF during Imag was less strong (r = −0.62) (Fig. 5). Figure 2. Decoding accuracy for HRV time and frequency domain during Rest, Imagination (Imag) and self- induced cognitive trance (SICT) between two class problem by linear discriminant (LDA) classifier; *: 5% across maximum statistic distribution. 7 Vol.:(0123456789) Scientific Reports | (2023) 13:15811 | https://doi.org/10.1038/s41598-023-42393-7 www.nature.com/scientificreports/ Figure 3. Comparison of thoracic respiratory rate, amplitude, phase and variability during Rest (green), Imagination (Imag, yellow) and Self-Induced Cognitive Trance (SICT, purple). RRV respiratory rate variability, SD BB standard deviation of BB (breath to breath) interval. RMSSD the root mean square of standard deviation of BB interval, a.u arbitrary unit, I/E inspiration/expiration (corrected for multiple comparisons, *p = 0.05; **p = 0.01). Figure 4. Decoding accuracy for respiratory features during Rest, Imagination (Imag) and self-induced cognitive trance (SICT) between two class problem by linear discriminant classifier; *: 5% across maximum statistic distribution. 8 Vol:.(1234567890) Scientific Reports | (2023) 13:15811 | https://doi.org/10.1038/s41598-023-42393-7 www.nature.com/scientificreports/ Discussion The objective of this study was to examine the physiological mechanisms of autonomic nervous system underly- ing SICT, focusing specifically on heart rate and breathing modulation during this non-ordinary state of con- sciousness. We compared the cardio-respiratory changes between ordinary consciousness states (resting state and imagination) and SICT. Our findings revealed significant differences in HRV and respiratory variability during SICT compared to the control conditions. Specifically, SICT was associated with an increase in heart rate, an increase in HRV time domain, a decrease in HF in the frequency domain of HRV, and a decrease in the entropy of the cardiac signal (ApEn). These changes were also characterized by a phasic decrease in HF-HRV. Additionally, respiratory changes indicated an increase in respiratory amplitude and variability. Importantly, no significant differences were observed between the imagination and resting state conditions. These changes suggest a parasympathetic modulation of the autonomic nervous system, referred to as a vagolitic effect 27 , which implies an inhibition of the vagus nerve. The vagus nerve is the primary nerve of the parasympathetic nervous system, consisting of 80% afferent sensory fibers and 20% efferent motor fibers27,56 . The efferent part of the vagus nerve that affects the heart is known as CVC and can be directly monitored through HF 57 . CVC is associated with a wide range of positive outcomes related to executive functions, emotion, and health, indicating overall better self-regulation of the organism 57,58 . The observed withdrawal of HF during SICT is surprising, as it contrasts with a significant body of research demonstrating the detrimental effects of HF reduction. Decreased HF is considered predictive of poorer health outcomes, increased risk of cardiovascular disease and mortality 59 , chronic stress 60 , depression 61 , phobia 62 , schizophrenia and post traumatic syndrome disorders 63 . It has even been proposed that HF could serve as a psychopathology biomarker related to the autonomic nervous system 64 . Our data support a withdrawal of HF during SICT, indicating a modulation of the nervous system associated with a decrease in parasympathetic activity. This modulation is observed from the resting state and continues dur- ing the task condition of SICT. These results can be interpreted in light of the vagal tank theory20 , which suggests that self-regulation is based on vagal control at rest (also known as tonic CVC). This theory complements the Figure 5. Individual difference and correlation between phasic and tonic HF during SICT and Imag conditions. Upper panel shows Tonic (i.e., Rest) and Phasic HF (i.e., HF-SICT minus HF-Rest or HF-Imag minus HF-Rest). Lower panel shows Spearman correlation between phasic and tonic HF for both conditions (SICT and Imag); HF high frequency, *p = 0.05. 9 Vol.:(0123456789) Scientific Reports | (2023) 13:15811 | https://doi.org/10.1038/s41598-023-42393-7 www.nature.com/scientificreports/ neurovisceral integration theory57 by adding the mechanism of self-regulation, which involves how vagal control responds to specific events or tasks (reactivity, also called phasic CVC) and how the system returns to normal after the event (recovery). CVC withdrawal is associated with responding to environmental demands, including metabolically demanding states such as exercise, stress, attention, and information processing 27 . Now, one might wonder why there is also a withdrawal of CVC observed during the imagination condition. According to the vagal tank theory 20 , one would have expected this effect to be specific to SICT since only this condition involves an environmental demand. These results could be attributed to the cognitive load required to induce imagination. This finding is consistent with the scientific literature conducted within the framework of the vagal tank theory, which suggests that self-regulation during tasks requires increased executive cognitive demands, resulting in a lesser decrease or even an increase in vagal control 19,20 . Furthermore, in our study, no significant difference was found between the imagination and resting state conditions. This result suggests that the recall of a memory does not recruit the same physiological mechanisms or intensity at the autonomic nerv- ous system level compared to the experience of SICT itself. During SICT, we also observed a consistent and larger decrease in phasic CVC compared to the other two control conditions. Previous studies have found a decrease in CVC during the fight-or-flight stress response, which led to near complete vagal withdrawal 27,65 . Therefore, a possible interpretation of our results is that this parasympathetic withdrawal is associated with an acute stress response. We also found an increase in ampli- tude and RRV during SICT, which could also be associated with an acute stress response. However, there is an inconsistency related to the trend of decreasing respiratory rate, as both an increased rate and tidal volume are markers of the acute respiratory response to stress 66 . However, unlike a stress response, SICT is induced and can be stopped by will. Studies on inter-individual differences in physiological stress regulation have shown that a decrease in HRV is associated with an increase in stress-induced cortisol67 . A plausible hypothesis is that regular practice of autonomic nervous system modulation (i.e., increase or decrease of CVC) could lead to improved autonomic regulation, as already demonstrated by regular physical activity (decrease of CVC) or regular mindfulness meditation (increase of CVC)32,68 . Conse- quently, regular practice of SICT could enhance control over acute stress responses typically triggered by uncon- trolled external stimuli (i.e., the detection and interpretation of danger). Developing the ability to manage and control these stress response mechanisms through willpower could be beneficial for various clinical populations affected by psychopathologies associated with dysregulation of the physiological stress system, such as anxiety and post-traumatic stress disorder. The withdrawal of CVC could be a hyperactivation of the autonomic nervous system associated with a SICT- induced dissociation process, which is typical of other forms of non-ordinary states of consciousness8,9,69 . Certain forms of intense stress have been associated with a dissociation process in humans dissociative experience that is physiologically linked to the autonomic nervous system through the limbic system, which, along with the frontal cortex, forms the neural basis of the depersonalization model 70–72 . In the literature, non-ordinary states of consciousness can be classified at the physiological level as either an increase in parasympathetic tone, cor- responding to the relaxation model and also associated with hypnosis 73 or certain types of meditation 31,74,75 and ecstatic experience40 , or a decrease in parasympathetic tone, as found in SICT and certain types of meditation23,76 , and probably other types of trance 40 . The effects of SICT on the autonomic nervous system found in this study are not consistent with the scientific literature on other forms of non-ordinary states of consciousness, such as meditation, which typically induces an increase in parasympathetic tone 36 . However, when comparing autonomic nervous system changes under SICT with other forms of non-ordinary states of consciousness, it is important to consider the techniques used, the type of induction employed, the associated psychological state, and the reported phenomenology under these states. That said, some studies have shown varying effects of meditation on the autonomic nervous system 76 . For example, Amihai et al. (2014) compared two meditation styles (Theravada and Vajrayana) and found oppos- ing results in terms of HF power Theravada meditation, which includes practices involving sitting, focusing on breath, and monitoring attention, showed an increase in parasympathetic tone associated with the intensity of the practice. In contrast, Vajrayana meditation, which includes practices of high absorption level, dance, and emotional experience, showed the opposite effect, i.e., a parasympathetic withdrawal, with a decrease in HF power, similar to what we observed in SICT in our results 76 . Many similarities can be drawn between Vajrayana and SICT states, not only at the physiological level but also at the phenomenological or behavioral level (body signs), including rich multimodal activity, sensory-motor patterns, feelings and emotions, as well as dereifica- tion, self-decentering, and an altered sense of self 76 . Integrative studies comparing the effects of various forms of non-ordinary consciousness on autonomic nervous system activity would, therefore, be necessary to confirm their differences and similarities from a neurobiological point of view. Limitations and future studies Our study has several limitations that should be acknowledged. Firstly, the generalizability of our results is limited due to the small sample size and the imbalanced female-male ratio. Replication of our findings in a larger and more diverse population is necessary. Secondly, the naturalistic design of our study restricts the ability to make long-term predictions based on our findings. Future research should aim to investigate the long-term effects of SICT on autonomic nervous system activity in a more controlled setting. Thirdly, comparing long-term basal values with the phenomenological content and emotional content experienced during SICT would provide a more comprehensive understanding of the physiological impact of repetitive practices. Furthermore, integrating interoceptive and brain data would be valuable in elucidating autonomic self-reg- ulatory mechanisms. Future studies should explore the relationship between autonomic nervous system modu- lation and neuro-psychological interactions to better comprehend why a hyperarousal state of the autonomic 10 Vol:.(1234567890) Scientific Reports | (2023) 13:15811 | https://doi.org/10.1038/s41598-023-42393-7 www.nature.com/scientificreports/ nervous system can either contribute to or be indicative of psychological states such as dissociation or absorption. Assessing the activity of brain regions known to activate the physiological stress system, such as the amygdala, in association with SICT would be of interest. Additionally, evaluating the activity of the hypothalamic–pitui- tary–adrenal axis during SICT could further confirm the hypothesis of physiological stress system involvement. In summary, addressing these limitations and incorporating additional measures and analyses in future studies would enhance our understanding of the physiological and psychological aspects of SICT. Conclusion This study is the first to provide insights into the physiological component of the autonomic nervous system during SICT. Our findings reveal that SICT is characterized by a decrease in cardiac vagal control, leading to a hyperarousal state of the autonomic nervous system. This physiological response bears resemblance to the stress response observed in other forms of non-ordinary states of consciousness, such as meditation, and may involve the autoregulatory systems of the cortex and limbic regions. The induction of SICT appears to be associated with alterations in the amplitude, phase duration, and variability of the respiratory rhythm. Overall, this research emphasizes the central role of the autonomic nervous system in SICT as an embodied experience. It invites further research investigations into the regular practice of SICT and its potential implications on the underlying physiological mechanisms of stress regulation. Data availability The complete dataset is available upon reasonable request (avanhaudenhuyse@chuliege.be; ogosseries@uliege.be). Received: 31 May 2023; Accepted: 9 September 2023 References 1. Flor-Henry, P., Shapiro, Y. & Sombrun, C. Brain changes during a shamanic trance: Altered modes of consciousness, hemispheric laterality, and systemic psychobiology. Cogent. Psychol. 4, 1313 (2017). 2. Lutz, A., Jha, A. P., Dunne, J. D. & Saron, C. D. Investigating the phenomenological matrix of mindfulness-related practices from a neurocognitive perspective. Am. Psychol. 70, 632–658 (2015). 3. Dorjee, D. Defining contemplative science: The metacognitive self-regulatory capacity of the mind, context of meditation practice and modes of existential awareness. Front. Psychol. 7, 1–15 (2016). 4. Lifshitz, M., van Elk, M. & Luhrmann, T. M. Absorption and spiritual experience: A review of evidence and potential mechanisms. Conscious. Cogn. 73, 102760 (2019). 5. Tellegen, A. & Atkinson, G. Openness to absorbing and self-altering experiences (“absorption”), a trait related to hypnotic suscep- tibility. Hypn. Theory Res. Appl. 83, 223–232 (1974). 6. Demertzi, A., Vanhaudenhuyse, A., Noirhomme, Q., Faymonville, M. E. & Laureys, S. Hypnosis modulates behavioural measures and subjective ratings about external and internal awareness. J. Physiol. Paris 109, 173–179 (2015). 7. Lynn, S. J. & Green, J. P. The sociocognitive and dissociation theories of hypnosis: Toward a rapprochement. Int. J. Clin. Exp. Hypn. 59, 277–293 (2011). 8. Van Der Kruijs, S. J. M. et al. Neurophysiological correlates of dissociative symptoms. J. Neurol. Neurosurg. Psychiatry 85, 174–179 (2014). 9. Vanhaudenhuyse, A. et al. Can subjective ratings of absorption, Dissociation, and time perception during ‘neutral hypnosis” predict hypnotizability ? An exploratory study. Int. J. Clin. Exp. Hypn. 67, 28–38 (2019). 10. Varela, F. J. Neurophenomenology: a methodological remedy for the hard problem. J. Conscius. Stud. 3, 330 (1996). 11. Husserl, E. On the phenomenology of the consciousness of internal time. Trans. J. B. Brough Dordr. Kluwer 4, 117 (1917). 12. Papies, E. K., Pronk, T. M., Keesman, M. & Barsalou, L. W. The benefits of simply observing: Mindful attention modulates the link between motivation and behavior. J. Pers. Soc. Psychol. 108, 148–170 (2015). 13. Bernstein, A. et al. Decentering and related constructs: A critical review and metacognitive processes model. Perspect. Psychol. Sci. 10, 599–617 (2015). 14. Fletcher, L. & Hayes, S. C. Relational frame theory, acceptance and commitment therapy, and a functional analytic definition of mindfulness. J. Ration. Emotive Cogn. Behav. Ther. 23, 315–336 (2005). 15. Wehrwein, E. A., Orer, H. S. & Barman, S. M. Overview of the anatomy, physiology, and pharmacology of the autonomic nervous system. Compr. Physiol. 6, 1239–1278 (2016). 16. Quigley, K. S., Kanoski, S., Grill, W. M., Barrett, L. F. & Tsakiris, M. Functions of interoception: From energy regulation to experi- ence of the self. Trends Neurosci. 44, 29–38 (2021). 17. Barrett, L. F. & Simmons, W. K. Interoceptive predictions in the brain. Nat. Rev. Neurosci. 16, 419–429 (2015). 18. Hutchinson, J. B. & Barrett, L. F. The power of predictions: An emerging paradigm for psychological research. Curr. Dir. Psychol. Sci. 28, 280–291 (2019). 19. Laborde, S., Mosley, E. & Mertgen, A. A unifying conceptual framework of factors associated to cardiac vagal control. Heliyon 4, e01002 (2018). 20. Laborde, S., Mosley, E. & Mertgen, A. Vagal Tank theory: The three Rs of cardiac vagal control functioning—Resting, reactivity, and recovery. Front. Neurosci. 12, 1–14 (2018). 21. Kim, D.-K., Lee, K.-M., Kim, J., Whang, M.-C. & Kang, S.-W. Dynamic correlations between heart and brain rhythm during autogenic meditation. Front. Hum. Neurosci. 7, 1–8 (2013). 22. Yüksel, R., Ozcan, O. & Dane, S. The effects of hypnosis on heart rate variability. Int. J. Clin. Exp. Hypn. 61, 162–171 (2013). 23. Kox, M. et al. Voluntary activation of the sympathetic nervous system and attenuation of the innate immune response in humans. Proc. Natl. Acad. Sci. USA 111, 7379–7384 (2014). 24. Russo, M. A., Santarelli, D. M. & O’Rourke, D. The physiological effects of slow breathing in the healthy human. Breathe 13, 298–309 (2017). 25. Grossman, P. & Taylor, E. W. Toward understanding respiratory sinus arrhythmia: Relations to cardiac vagal tone, evolution and biobehavioral functions. Biol. Psychol. 74, 263–285 (2007). 26. Muzik, O., Reilly, K. T. & Diwadkar, V. A. “Brain over body”—A study on the willful regulation of autonomic function during cold exposure. Neuroimage 172, 632–641 (2018). 27. Porges, S. W. The polyvagal perspective. Biol. Psychol. 74, 116–143 (2007). 28. Sevoz-Couche, C. & Laborde, S. Heart rate variability and slow-paced breathing: When coherence meets resonance. Neurosci. Biobehav. Rev. 135, 104576 (2022). 11 Vol.:(0123456789) Scientific Reports | (2023) 13:15811 | https://doi.org/10.1038/s41598-023-42393-7 www.nature.com/scientificreports/ 29. Laborde, S. et al. Effects of voluntary slow breathing on heart rate and heart rate variability: A systematic review and a meta-analysis. Neurosci. Biobehav. Rev. 138, 104711 (2022). 30. Aysin, B. & Aysin, E. Effect of respiration in heart rate variability (HRV) analysis. Annu Int. Conf. IEEE Eng. Med. Biol. Proc. https:// doi.org/10.1109/IEMBS.2006.260773 (2006). 31. Soni, R. & Muniyandi, M. Breath rate variability: A novel measure to study the meditation effects. Int. J. Yoga 12, 45 (2019). 32. Tyagi, A. & Cohen, M. Yoga and heart rate variability: A comprehensive review of the literature. Int. J. Yoga 9, 97 (2016). 33. Aubert, A. E., Verheyden, B., Beckers, F., Tack, J. & Vandenberghe, J. Cardiac autonomic regulation under hypnosis assessed by heart rate variability: Spectral analysis and fractal complexity. Neuropsychobiology 60, 104–112 (2009). 34. Huhle, R. et al. Effects of awareness and nociception on heart rate variability during general anaesthesia. Physiol. Meas. 33, 207–217 (2012). 35. Arya, N. K., Singh, K., Malik, A. & Mehrotra, R. Effect of heartfulness cleaning and meditation on heart rate variability. Indian Heart J. 70, S50–S55 (2018). 36. Santarcangelo, E. L. et al. Hypnotizability modulates the cardiovascular correlates of subjective relaxation. Int. J. Clin. Exp. Hypn. 60, 383–396 (2012). 37. Diamond, S. G., Davis, O. C. & Howe, R. D. Heart-rate variability as a quantitative measure of hypnotic depth. Int. J. Clin. Exp. Hypn. 56, 1–18 (2008). 38. Lane, J. D., Seskevich, J. E. & Pieper, C. F. Brief meditation training can improve perceived stress and negative mood. Altern. Ther. Health Med. 13, 38–44 (2007). 39. Toneatto, T. & Nguyen, L. Does mindfulness meditation improve anxiety and mood symptoms? A review of the controlled research. Can. J. Psychiatry 52, 260–266 (2007). 40. Rouget, G. La Musique et la Transe. (1990). 41. Lafon, Y. et al. De la transe chamanique à la transe cognitive : Revue de la littérature et approche biopsychosociale. In Yannick. 98–116 (2018). 42. Grégoire, C., Sombrun, C., Gosseries, O. & Vanhaudenhuyse, A. La Transe Cognitive Auto-Induite: Caracteristiques et Applications Thérapeutique Potentielles. (2021). 43. Peters, L. Shamanism: Phenomenology of a spiritual discipline. J. Transpers. Psychol. 21, 115 (1989). 44. Rock, A. J. Shamanic Journeying to the Lower World. Vol. 25. 45–55 (2006). 45. Hove, M. J. et al. Brain network reconfiguration and perceptual decoupling during an absorptive state of consciousness. Cereb. Cortex 26, 3116–3124 (2016). 46. Gosseries, O. et al. Behavioural and brain responses in cognitive trance: A TMS-EEG case study. Clin. Neurophysiol. Lett. Ed. 48, 550–554 (2020). 47. Huels, E. R. et al. Neural correlates of the shamanic state of consciousness. Front. Hum. Neurosci. 15, 1–16 (2021). 48. Makowski, D. et al. NeuroKit2: A Python toolbox for neurophysiological signal processing. Behav. Res. Methods 53, 1689–1696 (2021). 49. Sacha, J. et al. How to strengthen or weaken the HRV dependence on heart rate—Description of the method and its perspectives. Int. J. Cardiol. 168, 1660–1663 (2013). 50. Gasior, J. S., Sacha, J., Jelen, P. J., Zielinski, J. & Przybylski, J. Heart rate and respiratory rate influence on heart rate variability repeatability: Effects of the correction for the prevailing heart rate. Front. Physiol. 7, 1–11 (2016). 51. Gasior, J. S. et al. Interaction between heart rate variability and heart rate in pediatric population. Front. Physiol. 6, 1–10 (2015). 52. Sacha, J. Interaction between heart rate and heart rate variability. Ann. Noninvasive Electrocardiol. 19, 207–216 (2014). 53. Combrisson, E. & Jerbi, K. Exceeding chance level by chance: The caveat of theoretical chance levels in brain signal classification and statistical assessment of decoding accuracy. J. Neurosci. Methods 250, 126–136 (2015). 54. Nichols, T. E. & Holmes, A. P. Nichols-2001.pdf. Vol. 25. 1–25 (2001). 55. Pantazis, D., Nichols, T. E., Baillet, S. & Leahy, R. M. A comparison of random field theory and permutation methods for the statistical analysis of MEG data. Neuroimage 25, 383–394 (2005). 56. Brodal, P. The Central Nervous System—Structure and Function. (2010). 57. Thayer, J. F., Hansen, A. L., Saus-Rose, E. & Johnsen, B. H. Heart rate variability, prefrontal neural function, and cognitive perfor- mance: The neurovisceral integration perspective on self-regulation, adaptation, and health. Ann. Behav. Med. 37, 141–153 (2009). 58. Thayer, J. F., Åhs, F., Fredrikson, M., Sollers, J. J. & Wager, T. D. A meta-analysis of heart rate variability and neuroimaging studies: Implications for heart rate variability as a marker of stress and health. Neurosci. Biobehav. Rev. 36, 747–756 (2012). 59. Thayer, J. F. & Lane, R. D. The role of vagal function in the risk for cardiovascular disease and mortality. Biol. Psychol. 74, 224–242 (2007). 60. Kuehl, L. K. et al. Two separable mechanisms are responsible for mental stress effects on high frequency heart rate variability: An intra-individual approach in a healthy and a diabetic sample. Int. J. Psychophysiol. 95, 299–303 (2015). 61. Schiweck, C., Piette, D., Berckmans, D., Claes, S. & Vrieze, E. Heart rate and high frequency heart rate variability during stress as biomarker for clinical depression. A systematic review. Psychol. Med. 49, 200–211 (2019). 62. Åhs, F., Sollers, J. J., Furmark, T., Fredrikson, M. & Thayer, J. F. High-frequency heart rate variability and cortico-striatal activity in men and women with social phobia. Neuroimage 47, 815–820 (2009). 63. Jung, W., Jang, K. I. & Lee, S. H. Heart and brain interaction of psychiatric illness: A review focused on heart rate variability, cogni- tive function, and quantitative electroencephalography. Clin. Psychopharmacol. Neurosci. 17, 459–474 (2019). 64. Beauchaine, T. P. & Thayer, J. F. Heart rate variability as a transdiagnostic biomarker of psychopathology. Int. J. Psychophysiol. 98, 338–350 (2015). 65. Beauchaine, T. P., Gatzke-Kopp, L. & Mead, H. K. Polyvagal Theory and developmental psychopathology: Emotion dysregulation and conduct problems from preschool to adolescence. Biol. Psychol. 74, 174–184 (2007). 66. Tipton, M. J., Harper, A., Paton, J. F. R. & Costello, J. T. The human ventilatory response to stress: rate or depth?. J. Physiol. 595, 5729–5752 (2017). 67. Pulopulos, M. M., Vanderhasselt, M. A. & De Raedt, R. Association between changes in heart rate variability during the anticipa- tion of a stressful situation and the stress-induced cortisol response. Psychoneuroendocrinology 94, 63–71 (2018). 68. Platisa, M. M., Mazic, S., Nestorovic, Z. & Gal, V. Complexity of heartbeat interval series in young healthy trained and untrained men. Physiol. Meas. 29, 439–450 (2008). 69. Peres, J. F., Moreira-Almeida, A., Caixeta, L., Leao, F. & Newberg, A. Neuroimaging during trance state: A contribution to the study of dissociation. PLoS ONE 7, 1–9 (2012). 70. Frewen, P. A. & Lanius, R. A. Toward a psychobiology of posttraumatic self-dysregulation: Reexperiencing, hyperarousal, dissocia- tion, and emotional numbing. Ann. N. Y. Acad. Sci. 1071, 110–124 (2006). 71. Sierra, M. et al. Autonomic response in depersonalization disorder. Arch. Gen. Psychiatry 59, 833–838 (2002). 72. Mula, M., Pini, S. & Cassano, G. B. The neurobiology and clinical significance of depersonalization in mood and anxiety disorders: A critical reappraisal. J. Affect. Disord. 99, 91–99 (2007). 73. Fernandez, A., Urwicz, L., Vuilleumier, P. & Berna, C. Impact of hypnosis on psychophysiological measures: A scoping literature review. Am. J. Clin. Hypn. 64, 36–52 (2022). 74. Krygier, J. R. et al. Mindfulness meditation, well-being, and heart rate variability: A preliminary investigation into the impact of intensive vipassana meditation. Int. J. Psychophysiol. 89, 305–313 (2013). 12 Vol:.(1234567890) Scientific Reports | (2023) 13:15811 | https://doi.org/10.1038/s41598-023-42393-7 www.nature.com/scientificreports/ 75. Brown, L. et al. The effects of mindfulness and meditation on vagally mediated heart rate variability: A meta-analysis. Psychosom. Med. 83, 631–640 (2021). 76. Amihai, I. & Kozhevnikov, M. Arousal vs. relaxation: A comparison of the neurophysiological and cognitive correlates of Vajrayana and Theravada meditative practices. PLoS ONE 9, 102990 (2014). Acknowledgements The study was supported by the University and University Hospital of Liege, the Belgian National Funds for Scientific Research (FRS-FNRS), the BIAL Foundation, the European Union’s Horizon 2020 Framework Pro- gramme for Research and Innovation under the Specific Grant Agreement No. 945539 (Human Brain Project SGA3), the MIS FNRS (F.4521.23), the fund Generet, the King Baudouin Foundation, the Télévie Foundation, the Belgium Foundation Against Cancer (Grants Number 2017064 and C/2020/1357), the Benoit Foundation (Bruxelles), the ERA-Net FLAG-ERA JTC2021 project ModelDXConsciousness (Human Brain Project Part- nering Project), the European Space Agency (ESA) and the Belgian Federal Science Policy Office (BELSPO) in the framework of the PRODEX Programme, the Public Utility Foundation ‘Université Européenne du Travail’, “Fondazione Europea di Ricerca Biomedica”, the Mind Science Foundation, the European Commission, the Fondation Leon Fredericq, the Mind-Care foundation, Wallonia as part of a program of the BioWin Health 619 Cluster framework. We are also thankful to all participants for agreeing to take part in this study. KJ is supported by funding from the Canada Research Chairs program (950-232368) and a Discovery Grant from the Natural Sciences and Engineering Research Council of Canada (2021-03426), a Strategic Research Clusters Program (2023-RS6-309472) from the Fonds de recherche du Quebec-Nature et technologies. OG is research associate and SL research director at FRS-FNRS. Author contributions Design of the study: O.G., A.V. Data collection: O.G., A.V., C.M., J.A. Data analysis: V.O., H.Y., J.K., O.G., A.V. Manuscript writing: V.O., O.G., A.V. Reviewing paper: C.M., J.A., V.O., Y.H., K.J., S.L., M.E.F., A.B., F.R., C.S. Competing interests CS is the founder of the TranceScience Research Insitute (Paris, France) and developed the SICT technique. Additional information Supplementary Information The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598- 023- 42393-7. Correspondence and requests for materials should be addressed to V.O., A.V. or O.G. Reprints and permissions information is available at www.nature.com/reprints. Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creativecommons. org/ licenses/ by/4. 0/. © The Author(s) 2023