- 1Laboratory of Physiology, Department of Translational Medicine, University of Eastern Piedmont, Novara, Italy
- 2Radiology Unit, Department of Translational Medicine, University of Eastern Piedmont, Novara, Italy
- 3Psychiatry Unit, Department of Translational Medicine, University of Eastern Piedmont, Novara, Italy
- 4Neurology Unit, Department of Translational Medicine, University of Eastern Piedmont, Novara, Italy
- 5Statistic Unit, Department of Translational Medicine, University of Eastern Piedmont, Novara, Italy
Background/Aim: Misophonia is a disorder characterized by reduced tolerance to specific sounds or stimuli known as “triggers,” which tend to evoke negative emotional, physiological, and behavioral responses. In this study, we aimed to better characterize participants with misophonia through the evaluation of the response of the autonomic nervous system to “trigger sounds,” a psychometric assessment, and the analysis of the neurological pathways.
Materials and methods: Participants included 11 adults presenting with misophonic disturbance and 44 sex-matched healthy controls (HCs). Following recently proposed diagnostic criteria, the participants listened to six “trigger sounds” and a “general annoyance” sound (baby crying) during a series of physiological tests. The effects were examined through functional magnetic resonance imaging (fMRI), the analysis of heart rate variability (HRV), and of galvanic skin conductance (GSC). The fMRI was performed on a 3T Scanner. The HRV was obtained through the analysis of electrocardiogram, whereas the GSC was examined through the positioning of silver-chloride electrodes on fingers. Furthermore, the psychometric assessment included questionnaires focused on misophonia, psychopathology, resilience, anger, and motivation.
Results: Participants with misophonia showed patterns of increased sympathetic activation in response to trigger sounds and a general annoyance sound, the low frequency (LF) component of HRV, the sympathetic index, and the number of significant GSC over the threshold, where the amplitude/phasic response of GSC was higher. The fMRI analysis provided evidence for the activation of the temporal cortex, the limbic area, the ventromedial prefrontal/premotor/cingulate cortex, and the cerebellum in participants with misophonia. In addition, the psychometric assessment seemed to differentiate misophonia as a construct independent from general psychopathology.
Conclusion: These results suggest the activation of a specific auditory-insula-limbic pathway at the basis of the sympathetic activation observed in participants with misophonia in response to “trigger and general annoyance sounds.” Further studies should disentangle the complex issue of whether misophonia represents a new clinical disorder or a non-pathological condition. These results could help to build diagnostic tests to recognize and better classify this disorder. The relevance of this question goes beyond purely theoretical issues, as in the first case, participants with misophonia should receive a diagnosis and a targeted treatment, while in the second case, they should not.
Introduction
The recent consensus work of experts has led to the definition of misophonia as a reduced tolerance to specific sounds or stimuli related to them (Ferrer-Torres and Giménez-Llort, 2022). These “trigger” stimuli can generate strong negative emotional, physiological, and behavioral reactions that are not commonly observed by the majority of people (Sweedo et al., 2022). For individuals with misophonia, it is difficult to distract themselves from trigger sounds, and they may experience a range of unpleasant consequences, such as suffering, distress, overall functioning impairment, and mental health problems. In individuals suffering from misophonia, symptoms should not be better explained by any co-occurring disorder, including psychiatric conditions or disorders, such as anxiety disorders, mood disorders, personality disorders, obsessive, compulsive related disorders, post-traumatic stress disorder, autism spectrum disorder, and attention deficit hyperactivity disorder.
Past studies showed that subjects with misophonia can exhibit a wide range of “triggers” (Daniels et al., 2020; Jager et al., 2020a), which can be different from one subject to another. Although each person can have his/her own “trigger” stimuli, there are some that, more than others, can better serve as “misophonic triggers.” In particular, sounds associated with oral functions like chewing, eating, and swallowing are the most often reported as misophonic triggers. Also, sounds produced by other people, such as pen clicking, keyboard typing, finger or foot tapping, and shuffling footsteps often serve as “triggers.”
In response to those “trigger” stimuli, individuals with misophonia may experience anger, irritation, anxiety, and aggressive impulses, as well as symptoms related to autonomic nervous system activation, such as motor tics and increased heart rate (Edelstein et al., 2013; Webber et al., 2014; Palumbo et al., 2018). In addition, patients often develop coping mechanisms, such as avoiding social situations in which the trigger stimuli might occur (Schröder et al., 2013; Taylor et al., 2014; Potgieter et al., 2019). To date, it has been clarified that misophonia is an affective auditory processing disorder (Edelstein et al., 2013; Erfanian et al., 2017; Kumar et al., 2017) associated with an increased number of connections, or strength of connections, between the limbic and sympathetic nervous systems, which can cause abnormal processing of sound stimuli. In participants with misophonia, increased functional connectivity within brain regions like the anterior insular cortex (AIC), the anterior cingulate cortex, the ventromedial prefrontal cortex, the posteromedial cortex (posterior cingulate and retrosplenial cortex), the hippocampus, and the amygdala (Kumar et al., 2017) could reflect abnormal salience attribution to misophonic stimuli (Schröder et al., 2019; Seeley, 2019). In particular, it is noteworthy that the ventromedial prefrontal cortex and posteromedial cortex are core parts of the default mode network (DMN) (Raichle et al., 2001), which are connected to AIC and are activated when subjects are engaged in internally directed thoughts and retrieval of memories (Huijbers et al., 2011).
Repeated exposure to the same cues will amplify the salience network activity. The causal mechanism of this phenomenon reflects a conditioned response in which the initially neutral stimulus is increasingly associated with intensified aversive emotions and augmented by increased vigilance (Schröder et al., 2019). Among those areas, the anterior cingulate cortex and insular activity have been linked to increased cardiovascular arousal (Critchley, 2005). This statement was confirmed by the study of Kumar et al. (2017), which showed that the augmented heart rate (HR) and galvanic skin response by “trigger sounds” were mediated by the activity of AIC in participants with misophonia.
Therefore, the nervous structures mentioned above could represent the anatomical basis of the symptoms related to the activation of the autonomic nervous system that characterizes misophonia and could explain the visceral responses associated with emotions (Schwartz et al., 1981).
Recently, increased connectivity between both the auditory and the visual cortex, and between those brain regions and the ventral premotor cortex, has been reported in response to trigger sounds in subjects with misophonia (Kumar et al., 2021). This is of interest as the ventral premotor cortex is responsible for orofacial movement (Grabski et al., 2012). These findings support a model of misophonia based on “hyper-mirroring” according to which sounds would be the “medium” through which actions of other people are excessively mirrored.
Although knowledge about misophonia has increased, a complete analysis of the neuropsychiatric features, of the autonomic nervous system activation and of the neuronal pathways in participants with misophonia (Sweedo et al., 2022) in response to “trigger” sounds, is still lacking. In particular, previous studies have focused either on only one of the aspects mentioned above or have only partially investigated the activation of the autonomic nervous system by recording a few parameters related to the orthosympathetic/parasympathetic nervous system. Furthermore, the sample size calculation was not performed.
In the present study, we aimed to better characterize misophonia through the integrated evaluation of the physiologic, psychiatric, and neurological correlates in response to a trigger sound protocol in participants with misophonia. In particular, our primary endpoint was the analysis of the orthosympathetic/parasympathetic balance during the application of six “trigger sounds” and one “general annoyance” sound. To reach our primary endpoint, the sample size was calculated on the basis of previously reported tests related to the activation of the orthosympathetic nervous system in participants with misophonia. The data obtained in the present study were correlated with the psychometric assessment and the anatomic pathways activated during the trigger sounds protocol, which was investigated through functional magnetic resonance imaging (fMRI).
Materials and methods
We performed a double-blind case-control study by comparing subjects complaining of misophonic disturbances (participants with misophonia, N = 11) to gender-matched healthy controls (HCs) (N = 44) on the following variables: the fMRI, the analysis of heart rate variability (HRV), galvanic skin conductance (GSC), and psychiatric assessment. The study was approved by the local Ethical Committee, Azienda Ospedaliero-Universitaria Maggiore della Carità, Novara (CE 81/18). Participant recruitment took place from April 2018 to January 2020; subjects (or their legal representatives) were asked to sign a written informed consent and were treated according to Good Clinical Practice principles (Declaration of Helsinki: 2013).
The participants with misophonia were recruited from the community through the online misophonia support group (Misofonia Italia in Facebook). The HCs were recruited through local advertisements. The recruitment of the participants was carried out through an interview conducted according to the model proposed by Sanchez and da Silva (Schröder et al., 2013; Sanchez and da Silva, 2018). The selection was executed by an experienced psychiatrist.
Other potentially comorbid medical conditions were investigated as well.
Inclusion criteria were as follows: age ≥14 years, presence of at least one misophonic symptom, and written informed consent by each subject or parents.
Exclusion criteria were as follows: diagnosis of cardiac or neurologic disease, diabetes, intellectual disability/dementia/cognitive impairment, autoimmune/inflammatory diseases, pregnancy, use of medications, such as α or β blockers, diuretics, calcium channels blockers, smoking habits (>1 cigarette day), use of alcohol/psychoactive substances, unwillingness to participate and/or to sign the written informed consent. Participants were recruited based on the evidence of near-normal hearing as documented from previous visits relating to occupational medicine.
The participants with misophonia and HCs were instructed to abstain from caffeinated beverages for 24 h and from moderate or strenuous physical activity for 48 h before the analyses.
Trigger sounds protocol
Before beginning the study, the subjects were given a detailed explanation of the experimental procedures. They were instructed to fast for at least 3 h before the beginning of the experiment and all the evaluations were conducted from 9 a.m. to 1.30 p.m. All the participants underwent a rest period lasting 15 min before the start of the procedures (psychometric assessment, HRV, GSC, and fMRI analysis), in a quiet room with a controlled temperature between 26 and 27°C. We have chosen this temperature because it is the one that is reported to be the most comfortable for patients subjected to HRV analysis (Liu et al., 2008). In addition, the results of the GSC measurement are not affected by the selected temperature range (Doberenz et al., 2011). All evaluations were performed in a blinded condition.
In each subject, six different stimulations with sounds that have been widely reported as “triggers” (Edelstein et al., 2013; Kumar et al., 2017; Sweedo et al., 2022) were administered in the following sequence to avoid any bias when interpreting the statistical results: crunchy sound, nails tapping, chewing sound, fast breathing, typing, and click pen sound. Baby crying was used as the “general annoyance” sound, too. The audio clips were selected from YouTube and were composed by means of Audacity software (free and open-source digital audio editor and recording application software licensed under GPL-2.0) which allowed us to remove the background noise and adjust the volume to 70 dB HL. In particular, the sounds were adjusted by selecting the individual parts to get a normalization, in order to have a uniform amplitude. This was done by normalizing the peak width between 0.1 and 0.5 dB HL. As for the measurement of the output volume in headphones, Audacity has its own dB measurement system. The soundtracks were administered to the participants with misophonia and HCs by Beats Solo2 Headphones, with a sensitivity of 115 dB/mW, and an impedance of 45 ohms.
Moreover, each sound did not exceed the 70 dB HL limit, so as to make the experimental conditions as similar as possible to the real context.
The pattern of stimulation for each sound was as follows: 5 s silence, 30 s sound, 10 s silence, 30 s sound, 10 s silence, 30 s sound, 10 s silence, 30 s sound, and 5 s silence (total period of stimulation for each sound: 2 min and 40 s; five periods of silence and four sound applications; Figure 1).
Figure 1. Flowchart about the trigger sounds protocol. The entire sounds stimulation protocol was made of seven patterns of stimulation (six “trigger sounds” and one “annoyance sound”). All patterns were included in a single track. We measured the baseline GSC and HRV variables during each period of stimulation for comparison with the same values at baseline. GSC, galvanic skin conductance; HRV, heart rate variability.
This pattern of stimulation was preceded and followed by 2 min and 20 s, to allow the return to the baseline before the following stimulation, respectively.
The trigger sounds protocol was similar, in terms of duration and repetition of each sound/silence period, for the psychiatric assessment and the analysis of HRV and GSC. The only difference was during fMRI, where the silence period was 30 s instead of 10 s; this was due to the technical specifications of the fMRI equipment. Thus, the overall time of stimulation was about 42 min.
Psychometric assessment
Each participant was assessed with a protocol of psychometric tests, including both clinician-rated and self-rated tools, as detailed below. The psychometric assessment was performed by an experienced psychiatrist, trained in the use of the assessment tools described below.
We chose to include several of the available measures for misophonia, as most of them are not fully validated and show intrinsic limitations (for instance, they do not measure the actual magnitude of response to triggers), which is the reason why there have been recent attempts to develop new tools in this field of research (Dibb et al., 2021).
As possible associations have been suggested between misophonia and psychiatric disorders and/or symptoms, including anxiety, depression, and obsessive-compulsive tendencies, in our psychometric assessment, we included an overall measure of psychopathology (Symptom Checklist 90-R) and specific measures for the symptoms detailed above. Some of these measures are clinician-rated to avoid the possible biases of self-rated questionnaires. Last, we included a measure of resilience, which can be defined as the process of adapting well in the face of adversity and stress (American Psychological Association, 2014).
While we clearly expected to find higher scores on the misophonia scales in participants with misophonia when compared with controls, two scenarios could be hypothesized for the psychopathological measures: no difference between participants with misophonia and controls, or higher scores in participants with misophonia compared to controls, which would suggest a possible overlap of the misophonia construct with psychopathological symptoms. Regarding resilience, we were interested in assessing possible differences between the two groups, groups with possibly lower resilience ratings suggesting poorer resources in participants with misophonia in adapting to perceived stressors.
Amsterdam Misophonia Scale
This scale is an adaptation of the Yale-Brown Obsessive Compulsive Scale (Y-BOCS) to address misophonic symptoms. It has six items, scored from 0 (no symptom) to 4 (extreme). Scores 0–4 suggest subclinical misophonic symptoms, 5–9 mild, 10–14 moderate, 15–19 severe, and 20–24 extreme (Schröder et al., 2013).
Misophonia Questionnaire
The Misophonia Questionnaire (MQ) is a three-part self-rated questionnaire that evaluates the presence of misophonic symptoms, the resulting emotions and behaviors, and the overall severity of sound sensibility. Misophonia symptoms and resulting emotions and behaviors are assessed with 19 items which can be scored from 0 (not at all true) to 4 (always true). Ratings from the first two parts are summed together to produce a total score ranging from 0 to 72. Misophonia sound sensitivity severity is rated on a scale from 1 (minimal) to 15 (very severe) with a score greater than or equal to seven indicating clinical misophonia (McKay et al., 2018). Initial validation of the MQ indicates good internal consistency for the symptom scale (α = 0.86), the emotions and behaviors scale (α = 0.86), and total score (α = 0.89) (Wu et al., 2014).
Misophonia Activation Scale
This is a clinician-rated scale focusing on physical and emotional reactions to the trigger sounds in misophonic subjects. Participants are presented with 11 levels of responses to the known and personal misophonic triggers, with higher levels of response reflecting the severity of the disorder. The score ranges from 0 (no reaction to specific sound) to 10 (physical aggression, both self and others-directed) (Fitzmaurice, 2010).
Misophonia Assessment Questionnaire
The Misophonia Assessment Questionnaire (MAQ) includes 21 questions about the time spent on thoughts and feelings related to misophonic cues, evaluated on a four-points Likert-type scale ranging from 0 (not at all) to 3 (always true). Scores ranging from 0 to 11 indicate subclinical misophonic symptoms, 12–24: mild, 25–37 moderate, 38–50 severe, and 51–63 extreme (Johnson, 2014).
Yale-Brown Obsessive Compulsive Scale
Clinician-rated, 10-item scale, which assesses the severity of obsessions and compulsions in the week prior to the test. Total scores range from 0 to 40 (Goodman et al., 1989).
Hamilton Anxiety Rating Scale
This clinician-rated scale consists of 14 items, measuring both the mental anxiety (psychological distress) and somatic anxiety (physical symptoms related to anxiety). The total score ranges from 0 to 56, where <17 indicates mild severity, 18–24 mild to moderate severity, and 25–30 moderate to severe symptoms (Hamilton, 1959).
Hamilton Depression Rating Scale
The Hamilton Depression Rating Scale (HAM-D) is a clinician-rated 24-items test. The evaluation for most items is the result of the integration between the objective observation of signs and subjective exposure of symptoms, although the severity criterion mainly refers to the former. Scores 0–9 indicate subclinical depression, 10–13 mild depression, 14–17 moderate depression, and greater than17 indicate severe depression (Hamilton, 1960).
State-Trait Anxiety Inventory 1 and 2
It is a 40-item self-administered test for the assessment of state and trait anxiety. Each item is rated from 1 to 4 (1 = not at all, 4 = severe), and no specific cut-offs exist. The higher the score, the higher is the anxiety. Internal consistency coefficients for the scale ranged from 0.86 to 0.95; test–retest reliability coefficients ranged from 0.65 to 0.75 over a 2-month interval. Test–retest coefficients for this measure in the present study ranged from 0.69 to 0.89 (Spielberger, 1989).
Resilience Scale for Adults
It is a 33-item self-administered scale evaluating intra- or inter-relational stress-preventing factors (positive self-perception, positive future perception, social competence, structured style, family cohesion, and social resources). The higher the total score, the greater is the subject’s resilience. The Resilience Scale for Adults (RSA) is reliable, with good internal consistency demonstrated by Cronbach alpha values ranging from 0.79 to 0.88 in various studies, while among the six factors, it ranges from 0.67, for the structured style, to 0.81, for the perception of Self (Friborg et al., 2005).
Symptom Checklist 90 Revised
The Symptom Checklist 90 Revised (SCL-90-R) is a multidimensional self-report measure that assesses the severity of current psychological symptoms and distress. It assesses nine symptom dimensions: somatization, obsessive–compulsive, interpersonal sensitivity, depression, anxiety, hostility, phobic anxiety, paranoid ideation, and psychoticism. It also includes three global indices of psychological distress: global severity index (number of symptoms endorsed and intensity of distress), positive symptom distress index (average level of distress for those items that were endorsed; exaggerating or attenuating response style), and positive symptoms total (total symptoms endorsed/breadth of distress). All the Italian version subscales show a good internal coherence, with α-values between 0.70 and 0.96 (Derogatis, 1992).
Analysis of heart rate variability and galvanic skin conductance
In order to achieve our primary endpoint, which was related to the evaluation of the orthosympathetic/parasympathetic balance, we performed analyses of HRV and GSC, which are widely adopted methods for the assessment of the autonomic nervous system (Mackersie and Calderon-Moultrie, 2016). In particular, we expected to find increased parameters of HRV and GSC related to orthosympathetic activation in the participants with misophonia and increased parameters of HRV and GSC related to parasympathetic activation in the HCs.
Heart rate variability and GSC were monitored before (baseline) and during the trigger sounds protocol at the Laboratory of Physiology. In particular, we have extrapolated the numerical values of the entire registrations through special programs, as reported below. Thereafter, values were taken from the baseline and from the end of the trigger sounds protocol and were used to calculate an average value (Figures 2–4) and to make the graphs. All this was done for each sound and for each participant. The HRV was obtained through automatized analysis of electrocardiogram (ECG) by Kubios.
Figure 2. Effects of the trigger sounds protocol on HRV variables in participants with misophonia and participants without misophonia [healthy controls (HCs)] vs. specific baseline (set as 1). The results show increased orthosympathetic variables (LF, LF/HF) in participants with misophonia and increased parasympathetic variables (HRV triangular index, RMSSD) in participants without misophonia in comparison with values registered at baseline. It is also to note that the parasympathetic variables (HF, HRV, and RMSSD) were decreased in participants without misophonia RMS in comparison with values registered at baseline LF, low frequency (A); HF, high frequency (B); LF/HF (C), low frequency/high frequency ratio; HRV (D), heart rate variability triangular index; RMSSD (E), the root mean square of successive differences between normal heart beats. The results are the mean ± SE. The parenthesis indicate significance between groups, as specified in panel (F). In panel (F), the explanation for various groups’ representation is reported. The statistical analysis between M and HCs was performed through the Mann–Whitney test, whereas the statistical analysis between M/HCs and baseline was performed through the Wilcoxon test.
Figure 3. Effects of the trigger sounds protocol on HRV variables in participants with misophonia and participants without misophonia [Healthy controls (HCs)] vs. specific baseline (set as 1). The results show increased orthosympathetic variables (SNS index, Stress index, and VLF) in participants with misophonia in comparison with values registered at the baseline. It is also to note that the parasympathetic variables (PNS index and mean RR) were decreased in participants without misophonia in comparison with values registered at baseline. SNS (A), sympathetic; PNS (B), parasympathetic; stress index (C), stress index; RR (D), interval the elapsed time between ECG R waves; VLF (E), very low frequency; SDNN (F), standard deviation of RR intervals. The results are the mean ± SE. The parentheses indicate significance between groups, as specified in Figure 2F. Various groups are represented as in Figure 2. The statistical analysis between M and HCs was performed through the Mann–Whitney test, whereas the statistical analysis between M/HCs and baseline was performed through the Wilcoxon test.
Figure 4. Effects of the trigger sounds protocol on GSC variables in participants with misophonia and participants without misophonia [Healthy controls (HCs)] vs. specific baseline (set as 1). The results show increased orthosympathetic variables (CDA.nSCR, CDA.AmpSum, CDA.SCR, CDA.iSCR, CDA.PhasicMax, and CDA.Tonic) in participants with misophonia in comparison with values registered at the baseline (set as 1). In participants without misophonia, some of those variables related to orthosympathetic activation (CDA.AmpSum and CDA.PhasicMax) were lower in comparison with values registered at the baseline. CDA, continuous decomposition analysis; SCR, skin conductance response; CDA.nSCR (A), SCRs within response window (wrw); CDA.AmpSum (B), sum of SCR-amplitudes of significant SCRs; CDA.SCR (C), continuous decomposition analysis skin conductance response; CDA.iSCR (D), integral of SCR over the 10-s non-overlapping time; CDA.PhasixMax (E), maximum value of phasic activity wrw; CDA.Tonic (F), the mean tonic activity wrw. The results are the mean ± SE. The parentheses indicate significance between groups, as specified in Figure 2F. Various groups are represented as in Figure 2. The statistical analysis between M and HCs was performed through the Mann–Whitney test, whereas the statistical analysis between M/HCs and baseline was performed through the Wilcoxon test.
The data from ECG and GSC were recorded together on sitting subjects; furthermore, the subjects were seated in front of a white wall without any possible distraction and were asked to fixate on a point on the wall.
Heart rate variability analysis
Heart rate variability parameters were obtained by the computerized analysis of 6-lead ECG performed through Easy ECG pocket and Software Easy View Plus Stress, Ates Medical Device, Colognola Ai Colli, Verona, Italy. Kubios HRV program version 3.1.0 was used to analyze the collected data in both the time domain and frequency domain, as previously described (Tarvainen et al., 2014).
In particular, in the time domain, the following variables have been measured: mean RR (means of RR intervals, at RR intervals, or the elapsed time between ECG R waves), standard deviation of RR intervals (SDNN, for vagal activity) (Gernot, 2017; Forte et al., 2019), and the root mean square of successive differences between normal heart beats (RMSSD, for vagal activity). In the frequency domain, the following variables have been measured: high frequency (HF, for vagal activity), low frequency (LF) (LF, sympathetic activity, or a mix between sympathetic and vagal influences), LF/HF (sympathetic activity), very low frequency (VLF, for sympatho-vagal balance) and the HRV triangular index (for vagal activity).
In addition, sympathetic (SNS), parasympathetic (PNS) index, and the stress index have been examined through Kubios analysis of the registrations.
Galvanic skin conductance analysis
In order to collect GSC, a pair of Ag-AgCl electrodes was attached to the palmar surface of the middle and ring fingers of the participant’s non-dominant hand. Prior to attachment, participants’ hands were cleaned with an alcohol wipe and a skin conductance gel was applied to each electrode.
Galvanic skin conductance was recorded with a Shimmer3 GSR+ unit wireless device (Shimmer Research Ltd., Dublin, Ireland), with two 8 mm Ag chloride electrodes (GSR electrodes shimmer sensins), closed by velcro. Between electrodes, a constant potential of 0.5 V was maintained.
The recorder was connected via Bluetooth to a portable PC and the recorded data were shown on Consensys Pro (official Shimmer acquisition software) interface PC window, over a range 0–100 micro S (μS). The gain parameter was set at 5 μs/V and A/D resolution was set at 16 bits in order to acquire data over this time range. The sample acquisition rate was 30 samples/s (Gatti et al., 2018).
Analysis of skin conductance data was performed with Matlab version R2015a (The Mathworks, Natick, MA, United States) using Ledalab v3.4.7 and custom analysis programs. The raw skin conductance data were downsampled to 8 Hz and visually examined for the presence of motion artifacts, as indicated by HF fluctuations in the signal amplitude. The relatively few artifacts that were identified were replaced using linear interpolation. The GSC analysis was performed by using a continuous decomposition analysis (CDA) method in order to capture only the variations of the activity, under stimulation, that reflected an effective difference from baseline in each subject (Benedek and Kaernbach, 2010). In this case, period analysis was done only at every 10 s stimulation peak, and the results consisted of the average of GSC during the entire stimulation period.
The CDA was performed to separate the tonic component [skin conductance level (SCL)] and the phasic component [skin conductance response (SCR)] of the signal, using an amplitude criterion of 0.04 μS for defining SCRs, which reveals the activation of the sympathetic nervous system (Benedek and Kaernbach, 2010; Staib et al., 2015; Clark et al., 2018; Posada-Quintero et al., 2020).
Measured variables were as follows: the number of significant (above-threshold) SCRs within the response window (CDA.nSCR), the sum of SCR amplitudes of significant SCRs within the response window (reconvolved from the corresponding phasic driver-peaks) (CDA.AmpSum in μS), the mean tonic activity within the response window (of decomposed tonic component) (CDA.Tonic in μS), the maximum value of phasic activity within the response window (CDA.PhasicMax in μS), and the integral of SCR over the 10-s non-overlapping time (CDA.iSCR in μS × s), which was calculated to represent the overall SCR in a certain time period. The CDA.iSCR is thought to capture the cumulative effect of the signals while avoiding any influences by the usually arbitrary decision of the thresholds for peak detection and event definition (Zhang et al., 2018; Kim et al., 2019).
Functional magnetic resonance imaging analysis
The objective of the fMRI analysis was to describe the neural correlates of activation in participants with misophonia and to identify any abnormally activated cortical areas compared to healthy subjects. We expected to observe an activation of areas within the central nervous system involved in the control of the orthosympathetic/parasympathetic balance and in the processing of emotions in participants with misophonia during the trigger sounds protocol. Brodmann’s classification system was used for the definition of cortical areas.
Resonance magnetic imaging specificity and protocol
A 3-Tesla Magnetic Resonance scanner was used (Ingenia, Philips, Bergen, Norway) with a multichannel head coil (16 channels) and a NordicAktiva (NordicNeuroLab) fMRI system for stimulus administration.
The fMRI protocol was a continuous scanning protocol with two conditions, sound and silence. BOLD sequences, T1-weighted, and T2-weighted sequences were acquired for functional and anatomical reconstructions.
The parameters of T1-weighted sequence were as follows: FOV 130 × 130 × 120, voxel size = 2 × 2 × 4 mm, slice thickness = 4 mm, MPRAGE sequence, TR = 3000 ms, TE = 35 ms, flip angle = 90°, and acquisition matrix = 96 × 94. The parameters of BOLD sequences were as follows: FOV 182 × 240 × 256, voxel size = 1 × 1 × 1, slice thickness = 4 mm, TR = shortest, TE = shortest, flip angle = 8°, and acquisition matrix = 320 × 300.
The sounds were administered according to the trigger sounds protocol; during the fMRI continuous scanning protocol, for each sound, a sequence was repeated, which included four periods of sound presentation (30 s duration each) and 5 periods of silence (30 s duration each) where no sound was emitted, as described above. The total duration of this protocol was 41 min and 18 s (including centering and morphological T1 sequences).
Patient preparation
Each participant, before being positioned inside the MRI scanner, was instructed by the radiologist about the execution of the examination, the background sounds, the timing, and the protocol in use. No information was provided on the nature of the sound.
Before starting the acquisition, head coils and soundproofing headphones were used for the subjects (external noise reduction of about 30 dB). The chosen volume amounted to about 70 dB HL, which was set up before the starting fMRI through the machine software. This sound intensity reflects the average intensity value found in everyday life. Therefore, it represents a sound intensity that participants with misophonia may commonly encounter.
At the end of the fMRI scan, oral feedback on the perception of sounds was requested. All participants confirmed that they heard and recognized the sounds.
Post-processing analysis
Once the data were acquired, data were converted into DICOM format with the MRIcroGL software and then reconstructed with the SPM12 Software (2020) (SPM12, updated at October 2014). Within- and between-group comparisons were made in results obtained through the fMRI. The statistical performed parameter mapping (SPM) analysis was the same for all subjects and verified by an experienced operator external to the study.
Functional images were realigned to correct for motion, spatially transformed to standard stereotaxic space (based on the Montreal Neurological Institute coordinate system) and smoothed with an 8-mm full-width at half-maximum Gaussian kernel to decrease spatial noise prior to statistical analysis. Rotational motion in degrees (pitch, roll, and yaw), and translational movement in millimeters (x,y,z), were calculated based on the SPM12 parameters for motion correction of the functional images in each subject. No participants had movement greater than 2.5-mm translation or 3 degrees of rotation; therefore, none were excluded from further analysis. One-sample t-tests were performed to reveal statistically significant signal differences generated by trigger sounds and the general annoyance sound in participants with misophonia and HCs. The family-wise error rate corrected to p < 0.05 was used for all analysis methods (Nandy and Cordes, 2004).
During the fMRI model estimation process, reconstruction was performed for p-values = 0.1 (that were not statistically significant) without the family-wise error rate correction, as a control for an actual BOLD signal recording. All participants showed the presence of a signal; therefore, none were excluded from further analysis.
All reported areas having a statistically significant p-value and a cluster size greater than 2 were considered as “activated areas” or “activation.” No masks were applied during the post-processing analysis.
Statistical analyses
The calculation of the sample size was performed based on the primary endpoint, which in our study was defined as an activation of the autonomic nervous system. To do this, we considered a quantitative parameter related to the activity of the orthosympathetic system, that is GSC. Literature data reported a difference in the GSC between HCs and participants with misophonia amounting to 0.15 ± 0.4 μS (Edelstein et al., 2013). Using a statistical power of 80% and an alpha error of 0.05, we obtained a total of 55 subjects divided between 11 participants with misophonia and 44 HCs (controls). Misophonic and control groups were compared through descriptive statistics. Categorical data are reported as a percentage and absolute frequencies. Continuous data are summarized as mean ± standard deviation (SD) or standard error (SE). Wilcoxon-type and Mann–Whitney tests were performed for continuous variables and the Pearson chi-square test or Fisher-exact test, for categorical variables.
A generalized linear mixed effect (GLME) model was used to examine and compare response trends over time (during the trigger sounds protocol vs. baseline), compare group averages (HCs vs. participants with misophonia), and compare the measurement time within each group (interaction term). A random effect (intercept) on the subject’s identification code term has been considered to account for correlations within repeated measurements across time. The time effect with group membership interaction was also considered as a fixed-effect factor to adjust the estimates, together for sex and age. Separate univariable models were estimated. The p-values were adjusted for multiple endpoints by considering the Benjamini–Hochberg correction. The model estimates together with the SE and p-values are reported. The HRV and GSC results are presented in Figures 2–4 created using GraphPad Prism 6.
P-values < 0.05 were considered statistically significant. Statistical analyses were conducted using R 3.5.2 with rms packages (Core and Team, 2019) and lme4 (Bates et al., 2015).
Results
Data from participants with misophonia and HCs are reported in Table 1. The two groups were significantly different with regard to age, with the misophonic subjects being older than HCs.
Psychometric assessment
Table 2 shows significant group differences for the psychometric assessment between misophonic subjects and HCs. The most striking result emerging from the psychometric assessment is that, as expected according to the study sampling procedure, subjects in the misophonia group scored significantly higher on all the questionnaires specifically assessing the misophonia construct [Amsterdam Misophonia Scale (A-MISO-S), Misophonia Activation Scale (MAS-1), MQ, and MAQ]. In contrast, we failed to find group differences in the questionnaires assessing the overall psychopathology (Symptom Checklist-90-R), obsessive-compulsive, depressive, and anxiety symptoms [Y-BOCS, HAM-D, Hamilton Anxiety Rating Scale (HAM-A), and State-Trait Anxiety Inventory].
Specifically, according to the mean scores on the A-MISO-S and MAQ, participants in the misophonia group had mild levels of misophonia, while those in the HCs scored within normal limits. The mean MQ severity score in the misophonia group, nonetheless, was 5.20, lower than the cutoff (score > 7) indicating possible misophonia.
Another interesting result is the finding of lower resilience scores in misophonic subjects, suggesting that compared to HCs, they possess less resilience when facing stressors.
Effects of the trigger sounds protocol on heart rate variability and galvanic skin conductance
As reported in section “Materials and methods,” we performed this analysis in order to address our primary endpoint, which was to determine changes in the orthosympathetic/parasympathetic balance in the participants with misophonia vs. participants without misophonia.
First, we measured HRV and GSC variables in both participants with and without misophonia before the start of the trigger sounds protocol in order to collect baseline values of all participants. These data were used for comparison with the results obtained within each group of participants after the trigger sounds protocol. The changes between measurements of HRV and GSC variables obtained during the trigger sounds protocol vs. baseline values (set as 1) were calculated as a percentage. As shown in Table 3, the baseline analysis of HRV and GSC parameters evidenced significant differences between the participants with misophonia and HCs. Hence, all examined variables except the maximum value of phasic activity within the response window, LF, and LF/high frequency were higher in the participants with misophonia than HCs. In addition, high frequency, parasympathetic index, and VLF were lower in the formers. These findings were in agreement with our test hypothesis that the orthosympathetic tone (basal orthosympathetic activity) would be dominant in participants with misophonia, and a parasympathetic tone (basal parasympathetic activity) would be dominant in HCs in the baseline condition.
These baseline values were then used to compare the effects of the trigger sounds protocol in the two groups of subjects. It is noteworthy that in comparison with the baseline, in HCs, the trigger sounds protocol caused a significant increase of the HRV triangular index and a decrease of sympathetic index (as% vs. baseline: 11.7 ± 2.7 and −83.2 ± 31, respectively, p < 0.05; Figures 2, 3), the sum of SCR-amplitudes of significant SCRs within the response window, and the maximum value of phasic activity within the response window (as% vs. baseline: −16.36 ± 3.2 and −23.7 ± 2.5, respectively, p < 0.05; Figure 4). Also, the mean square difference between successive RR intervals showed an increase, which was, however, at the limit of significance (p = 0.05; Figure 3).
The analysis of these variables in participants with misophonia showed opposite results (Figures 2, 3). Moreover, the trigger sounds protocol resulted in a strong increase in all GSC parameters in misophonics compared to HCs (as percent vs. baseline, the number of significant, above-threshold, SCRs within the response window: 159.9 ± 12.1; the sum of SCR-amplitudes of significant SCRs within the response window: 141.3 ± 14.5; SCRs within the response window: 230.8 ± 33.1; the integral of SCR over the 10-s non-overlapping time: 101.6 ± 13.8; the maximum value of phasic activity within the response window: 87.8 ± 11.7; the mean tonic activity within the response window: 23.13 ± 1.8 vs. the number of significant, above-threshold, SCRs within the response window: −4.5 ± 6.8; the sum of SCR-amplitudes of significant SCRs within the response window:−16.3 ± 3.2; SCRs within the response window: −4.5 ± 64.9; the integral of SCR over the 10-s non-overlapping time: −9.3 ± 2.6; the maximum value of phasic activity within the response window: −23.7 ± 2.5; the mean tonic activity within the response window: −1.9 ± 1.2; p < 0.05; Figure 3), LF, LF/HF, sympathetic index, stress index and VLF (as% vs. baseline, LF: 64.8 ± 5.5; LF/HF: 271.7 ± 36.6; sympathetic index: 327.7 ± 119.2; stress index: 43.5 ± 5.3; VLF: 323.8 ± 77.5 vs. LF: −4.3 ± 1.6; LF/HF: −0.22 ± 2.6; sympathetic index: −83.2 ± 31; stress index: 12.77 ± 3.4; VLF: 0.22 ± 4.1; p < 0.05; Figures 2, 3). In addition, a reduction of all other HRV parameters was observed in participants with misophonia vs. HCs (Figures 2, 3).
An example of GSC analysis taken from one participant with misophonia is shown in Figure 5.
Figure 5. Example of GSC analysis performed during the trigger sounds protocol in a participant with misophonia. From top to bottom: Skin conductance raw data; Inter impulse activity used for continuous deconvolution analysis; Continuous deconvolution analysis-tonic; Continuous deconvolution analysis-phasic; Total reconstruction with differences between phasic and tonic.
A significant interaction emerged for all parameters indicating that the change in values between the baseline and the trigger sounds protocol showed different patterns for participants with misophonia and HCs. In particular, the values for sympathetic index, stress index, VLF, LF, and LF/HF increased over time among participants with misophonia, whereas, in the same subjects, the values for parasympathetic index, HF, and HRV decreased. When treated as a covariate for modeling group differences, “age” was only a significant factor in two of the N models tested: LF and HF. The GSC values showed increasing trends over time among participants with misophonia (Tables 4, 5).
On the whole, the results obtained provide evidence for a higher orthosympathetic and lower parasympathetic tone at baseline registrations in participants with misophonia than HCs. The state of these tones was then strengthened by the trigger sounds protocol.
Effects of the trigger sounds protocol on functional magnetic resonance imaging
Functional magnetic resonance imaging was used to characterize the anatomical pathways involved in physiological and emotional responses to the trigger sounds protocol. In particular, it was hypothesized that greater activation in brain regions involving the limbic system, the temporal cortex, and the cerebellum would be observed due to the relationship between these areas and the orthosympathetic/parasympathetic balance. The principal findings are summarized in Table 6.
Table 6. Functional magnetic resonance imaging (fMRI) activations in participants with misophonia and without misophonia.
The fMRI analysis confirmed our test hypotheses about the involvement of the above neuronal areas in the participants with misophonia during the trigger sounds protocol. Hence, in all cases, a strong activation of the temporal superior gyrus (BA 22; 100%), the temporal cortex (BA 21; 100%), and the auditory cortex (BA 41 and 42; 100%) was found. In about 70% of cases, the premotor cortex (BA 6) and the cerebellum were found to be activated, as well.
Furthermore, 8/11 (73%) participants with misophonia showed activation signals in at least one of the following: the hippocampus, the ventromedial prefrontal cortex, the insula, and the cingulate cortex.
Figure 6A shows the atypical activation of multiple brain regions in a participant with misophonia after hearing the “eating crunching” sounds. In particular, the auditory cortex (BA 41 and 42), the cerebellum, the insula, the premotor cortex (BA 6), and the frontal cortex (BA 4) were activated during the trigger sounds protocol. Another example of an fMRI image taken in one participant with misophonia is shown in Figure 6B. In that subject, the presentation of “breathing” caused the stimulation of the auditory cortex (BA 41 and 42), the cerebellum, and the premotor areas (BA 6). It is to note that those findings were evidenced after the application of the correction factor “family-wise error rate” in SPM analysis.
Figure 6. Brain areas activated during the trigger sounds protocol in two participants with misophonia in the three planes. The trigger sounds were “eating crunchy” in panel (A), and “breathing” in panel (B). In panel (A), details on left insula. Also, auditory cortex (BA 41 and 42; light blue arrow), the cerebellum (yellow arrow), the insula (green arrow), the premotor cortex (BA 6; red arrow), and the frontal cortex (BA 4; pink arrow) were activated during the trigger sounds protocol. In panel (B), the areas identified were the auditory areas (BA 41 and 42; light blue arrow), the cerebellum (yellow arrow), and the premotor cortex (BA 6; red arrow).
Concerning the HCs, 36/44 (82%) subjects demonstrated the activation of the auditory cortex (BA 41 and 42), while in 8 (18%) subjects, we could not find any activated area. In these eight subjects, further post-processing analysis was performed on the images, showing the presence of stimulus response for SPM with p = 0.1.
In Figure 7, examples of fMRI in two HCs are shown; in those subjects, “eating crunching” and “breathing” caused the activation of the primary auditory cortex (BA 42 and 41), only. Also, in this case, as specified above, the correction factor “family-wise error rate” in SPM analysis had been applied.
Figure 7. Activation of the auditory cortex (BA 41 and 42; light blue arrow) in two participants without misophonia. The stimulus administered was “eating crunchy” in panel (A), and “breathing” in panel (B).
Thus, the fMRI analysis highlighted the activation of brain regions involved both in the orthosympathetic response and in the processing of trigger stimuli only in participants with misophonia.
Discussion
The results of this study obtained by combining the fMRI analysis with physiological and psychiatric evaluations showed specific auditory-insula-limbic patterns of activation associated with increased sympathetic tone in participants with misophonia. Those findings are unlikely due to specific psychopathologic features, as shown by the absence of differences between the two groups of participants (those with misophonia and HCs) in the questionnaires assessing obsessive-compulsive, depressive, and anxiety symptoms.
As previously described, misophonia is a disorder characterized by reduced tolerance to specific sounds or stimuli known as “triggers,” which tend to evoke negative emotional, physiological, and behavioral responses (Kumar et al., 2017; Sweedo et al., 2022).
Recently collected data showed a prevalence of participants with misophonia that reaches up to 49% in the general population (Naylor et al., 2021) and about 52% in subjects with obsessive- compulsive disorder (Siepsiak et al., 2020). In addition to feelings of anxiety, distress, and anger, physical manifestations including tightness or pain in the entire body, increased muscular tone, dyspnea, tachycardia, and hypertension have been described (Cavanna and Seri, 2015).
Moreover, those experiences may lead, in some subjects, to a severe decline in daily functioning and the development of behavioral health problems.
Although misophonia now has a consensus definition (Sweedo et al., 2022), and some information is available on the utility of psychometric assessment and involvement of the autonomic nervous system, more in-depth investigations relating to the activation of the orthosympathetic and parasympathetic systems and to the concomitant activation of brain regions involved in emotional responses could be useful for a better understanding of this disorder.
Heart rate variability and galvanic skin conductance effects induced by the trigger sounds protocol in participants with misophonia
We utilized both HRV and GSC to analyze the autonomic nervous system balance in misophonic participants and HCs (Benedek and Kaernbach, 2010; Tarvainen et al., 2014; Seri, 2015; Staib et al., 2015; Gernot, 2017; Clark et al., 2018; Zhang et al., 2018; Forte et al., 2019; Kim et al., 2019; Posada-Quintero et al., 2020). In particular, we analyzed both the orthosympathetic and parasympathetic drive using HRV and analyzed the orthosympathetic profile using GSC, which is related to the cholinergic-dependent sympathetic stimulation of the cutaneous sweat glands.
The baseline registration of HRV and GSC variables showed interesting results in participants with misophonia. In these participants, HRV variables were more distinctly related to an increased orthosympathetic tone, such as LF and LF/HF and almost all GSC variables. Meanwhile, markers of parasympathetic tone given by HRV, including HF, parasympathetic index, and LF, were higher in HCs than in participants with misophonia. These findings add important information regarding the state of the autonomic nervous system in resting conditions in participants with misophonia vs. HCs and are useful for outlining their orthosympathetic/parasympathetic background.
Participants with misophonia showed increased levels for almost all orthosympathetic HRV measures and decreased levels for almost all parasympathetic measures. The activation of the orthosympathetic system was also confirmed by the GSC analysis. In HCs, it is noteworthy that trigger sounds resulted in a reduction of the sympathetic index, the sum of SCR-amplitudes of significant SCRs, and the maximum value of phasic activity. These results suggest that the trigger sounds protocol caused an inhibition of the orthosympathetic system in HCs unlike what was observed in the participants with misophonia.
Since the main bias of this study was the difference in age between participants with misophonia and HCs, we added a correction for age to the statistical analysis of HRV and GSC data. While we could not confirm our previous observations, since for LF and HF, the statistical significance in results was lost after the correction for age, the general patterns appeared to be similar to those observed without performing the correction for age.
While our findings were similar to those obtained by Kumar et al. (2017), there are several key differences in our experimental protocol that must be considered. For example, Kumar et al. (2017), used two trigger sounds and one neutral sound, whereas we concentrated on six trigger sounds and a general annoyance sound. Furthermore, changes in the autonomic nervous system balance were examined by analyzing many GSC and HRV variables, which evaluated the orthosympathetic/parasympathetic activity both at the baseline and during auditory stimulation. Finally, it should be noted that the sample size in this study was calculated through appropriate statistical tests based on the variation of a parameter related to the activation of the autonomic nervous system (Edelstein et al., 2013). In this way, we could answer the primary hypothesis of the study, which referred to the presence of an altered orthosympathetic/parasympathetic balance in participants with misophonia. In the study by Kumar et al. (2017), instead, participants with misophonia were matched without any power size calculation with participants without misophonia. In addition, in that study, results focused on fMRI analyses. Only one GSC measurement was performed (galvanic skin response), and one HRV measurement (heart rate); no psychometric assessment was provided.
In the study by Edelstein et al. (2013), participants with misophonia were recruited after an interview and the emotional responses to different “trigger sounds” were evaluated. In the second part of that study, a comparison was executed between six participants with misophonia and five HCs in terms of GSC response to auditory–visual stimulations. A similar kind of experimental protocol was executed by Schröder et al. (2019), who conducted psychometric assessments and fMRI analysis in 21 participants with misophonia and 23 HCs. However, a change in the HR was the only measure used to evaluate the autonomic nervous system activation.
As reported above, our findings about the activation of the orthosympathetic nervous system in participants with misophonia are in agreement with previous observations. However, here we conducted a more in-depth evaluation of autonomic nervous system involvement in response to a trigger sounds protocol, by combining it with psychometric and imaging evaluation.
Neuroanatomical pathways activated by the trigger sounds protocol in participants with misophonia
In participants with misophonia, we observed activations of Brodmann areas 21, 22, 41, and 42, which could be attributed to listening, detection, and understanding of a perceived sound (Pickles, 2015). In addition, in participants with misophonia, the trigger sounds protocol caused activation in the hippocampus, the cingulate cortex, the ventromedial prefrontal cortex, and the insula. Overall, those findings are in agreement with the data by Kumar et al. (2017), who showed that, in participants with misophonia, the trigger sounds were associated with abnormal functional connectivity between the anterior insula cortex, the ventromedial prefrontal cortex, the hippocampus, and the amygdala. It is to note that all of those brain regions are involved with the perception of interoceptive signals and emotion processing (McLachlan, 2009). Moreover, they play a role in the modulation of the autonomic nervous system and in the integration of physiological signals and the dynamic representation of emotional states. In particular, the findings of the involvement of the superior temporal cortex (BA 22), the area with the largest change in the activation in participants with misophonia, could be argued to be related to an increased sensitivity to the trigger sounds in those subjects. In this way, the trigger sounds protocol could have induced auditory attention and caused a greater response to the stimuli, resulting in signals being labeled as emotionally relevant. The absence of activation of the amygdala in participants with misophonia could be attributable to the possible inhibition exerted by the ventromedial prefrontal cortex, which would be involved in the regulation of emotions by the inhibition of neuronal areas involved in this process, such as the amygdala (Motzkin et al., 2015).
Alternatively, it is possible that we did not find any involvement of amygdala because this region is more associated with fear than representing a primary emotional site, as reported by Schröder et al. (2019).
The HCs group showed the activation of Brodmann areas 41 and 42, only, with a mean cluster size of 130 voxel. This could be explained by the trigger sounds protocol that, when compared to the background noise of the fMRI continuous protocol, was of low intensity. Moreover, in eight HCs, there was no activation. Those eight HCs stated during the brief interview following the fMRI that they perceived and recognized the sounds during the examination. These were subjects whose fMRI images were acquired on different days and did not belong to a single session. In order to exclude mistakes in subject sampling, reconstruction analyses were followed for SPM values with p-values above the significant threshold (p = 0.1). It is to note that this post-processing reconstruction showed activations of Brodmann areas 41 and 42. This activation, in eight HCs, also persisted for values of p = 0.05. However, the application of the family-wise error rate resulted in its elimination.
Psychiatric assessment: Comparison of participants with misophonia and healthy controls
The results of the psychiatric evaluation indicate that misophonia could be a construct independent of the general psychopathology (as assessed with the Symptom Checklist-90 R), anxiety (HAM-A, State-Trait Anxiety Inventory), depression (HAM-D), and obsessive-compulsive features (Y-BOCS).
The currently available literature has described possible associations of misophonia with a variety of psychiatric symptoms, such as traits of obsessive-compulsive personality disorder, mood disorders, attention-deficit (hyperactivity) disorder, autism spectrum conditions (Jager et al., 2020a), post-traumatic stress disorder, and anorexia (Erfanian et al., 2019). Nonetheless, data about the frequency of misophonia in patients suffering from psychiatric disorders are lacking (Siepsiak et al., 2020). Recent research supported the possible relationship of misophonia symptoms with clinician-rated symptoms of personality disorders, but not other psychiatric disorders, even though anxiety was found to partially mediate this relationship (Cassiello-Robbins et al., 2021).
There are also suggestions that misophonia could be a discrete psychiatric disorder, with corresponding implications for treatment (Schröder et al., 2013). Features, such as neuroticism, impulsivity, and difficulties with emotion regulation have been suggested as important risk factors and treatment targets for adults with misophonia (Cassiello-Robbins et al., 2020). Recently, both short-term and long-term efficacy of CBT for misophonia have been reported (Jager et al., 2020b).
In support of our baseline hypothesis, we found higher scores on the misophonia scales in participants with misophonia when compared to HCs. Nonetheless, it has to be underscored that participants in the misophonia group did not score above the cutoff established for misophonia in all scales; actually, they had mild misophonia according to the mean scores on the A-MISO-S and MAQ, but scored lower than the cutoff on the MQ. While this might lead one to argue whether the differences we found between the two groups could be clearly attributed to misophonia (Schröder et al., 2017), we also have to underscore that, as described in the “Materials and methods” section, our assessment included more than one questionnaire, the involvement of a psychiatrist, and proper screening in a face-to-face interview where information about co-morbidity was gathered.
As far as psychopathology is concerned, our first hypothesis is supported, as our results did not evidence the presence of psychopathologic symptoms in participants with misophonia (as assessed with the Symptom Checklist-90 R, HAM-A, State-Trait Anxiety Inventory, HAM-D, and Y-BOCS). Thus, the higher scores found in the misophonia measures seem independent of the overall psychopathology, leading to the possibility that misophonia may be either a discrete pathological condition or a non-pathological variant of human nature linked to a particular way of processing trigger sounds.
Also, misophonia could be involved in the creation of the so-called “global emotional moment.” This is meant as the set of homeostatic, environmental, hedonic, motivational, social, and cognitive activities (Bud, 2009), which contribute to one’s feelings and represent the sentient self at one moment in time.
Finally, the finding of lower resilience scores in participants with misophonia may suggest poorer resources in this group of subjects when facing stressors.
Connections between neuroanatomic pathways, autonomic nervous system response, and psychiatric assessment in participants with misophonia
The primary auditory cortex lies in the transverse temporal gyri of Heschl that are juxtaposed to the insula, which is considered a key brain area in the homeostasis of visceral information processing and interoception. In particular, the insula is involved in the control of the autonomic nervous system and in the integration of physiological signals, as well as the dynamic representation of emotional states to create the “global emotional moment” (McGeoch and Rouw, 2020). Negative emotional experiences have been reported to induce anterior insular activation, including disgusting, frightening, happy, sad, or sexual images (Uddin et al., 2017).
The cingulate cortex, which is connected to the amygdala, is responsible for the processing of emotions and for the regulation of associated endocrine and autonomic responses. This region is also involved in reward-related processing of endocrine and autonomic responses to them due to the fact that this area is involved in reward-related processing. The cingulate cortex can be considered a connecting center of emotions, sensations, and action. Due to its links with the hippocampus and amygdala, the cingulate cortex could also have a role in the consolidation of long-term memories and the processing of emotionally relevant stimuli. The integration of signals originating from all the above areas would, thus, result in the modulation of the autonomic nervous system drive and in processing emotional states. In predisposed subjects, an auditory-insular synesthesia model could account for the onset of psychiatric symptoms and clinical manifestation of changes in the orthosympathetic/parasympathetic balance. In this way, a dynamic process of altered neurological activation could turn the specific auditory stimuli into trigger sounds to induce strong emotional responses (Nagai et al., 2007; Rolls, 2019). In this context, our data would also support the concept of wellbeing that has been linked to the balance between the two sides of the autonomic nervous system. In particular, it has been hypothesized that conditions of chronic sympathetic hyperactivity and parasympathetic hypoactivity would be associated with reduced emotional wellbeing and a variety of mental disorders and vice versa (Thayer and Brosschot, 2005; Strigo and Craig, 2016).
The findings of increased baseline sympathetic variables in participants with misophonia and parasympathetic variables in HCs measured with HRV and GSC would confirm the above issues regarding the reduced “emotional wellbeing reserve” and increased predisposition to undergo changes in wellbeing in response to stressful auditory conditions.
With regard to resilience, which was found to be reduced in participants with misophonia, and in general emotional states, numerous studies have associated this trait with the activation of the mesocorticolimbic area, such as the hippocampus and the prefrontal cortex (Richter et al., 2019). Similarly, Tan et al. (2018) correlated heartbeat interoception and anxious state and found the activation of the insula to be related to this emotional state. Cholinergic signaling in the hippocampus has also been described to regulate the social stress resilience and anxiety- and depression-like behavior. Finally, this study found the activation of the cingulate cortex activation in misophonia, which would be involved with homeostatic motivations that guide an adaptive behavior. Taken together, the insula and limbic areas could represent the neuroanatomical basis linking orthosympathetic system activation, emotion, and feelings (Strigo and Craig, 2016).
A novel finding from this study involved motor control, such as Brodmann area 6 and the cerebellum. The premotor cortex (BA 6), in particular, would be activated “during motor imagery,” that is, when the subject visually imagines a movement or imagines the sensations, he would experience during that movement. Thus, it could be hypothesized that in participants with misophonia, the trigger sounds protocol would activate the neural pathways implicated in the preparation/execution of escape (Kumar et al., 2021).
Conclusion
In conclusion, the results obtained from our study underline the existence of a specific response pattern within the auditory cortex-insula-limbic areas which, in the presence of “trigger sounds,” would be activated in predisposed subjects, such as participants with misophonia. The recruitment of those neuronal patterns would, in turn, alter the autonomic nervous system balance in favor of the orthosympathetic drive, influencing the emotional wellbeing.
Considering the limitations of the study, future work should investigate different activation patterns resulting from various trigger sounds and the general annoyance sound. Multimodal stimulation using video clips and sounds that simulate everyday life may also reveal patterns similar to those described in this study. Moreover, future comparisons should include a “no annoyance control condition.” The selection of participants should also be based on the recent consensus definition of misophonia (Sweedo et al., 2022), which was not available at the time of this study, and the enlarged sample size to determine specific effects of gender. Finally, any bias related to differences in age could be avoided by specific selection of participants, which would allow for an age correction (this was not performed in the current study).
Data availability statement
The original contributions presented in this study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Ethics statement
The study was approved by the local Ethical Committee, Azienda Ospedaliero-Universitaria Maggiore della Carità, Novara (CE 81/18). Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author contributions
EG, AS, CG, AC, and PZ: conception of the work. DD, RC, BG, DN, and AF: data collection. AS, DD, BG, DN, DA, DF, and AC: data analysis and interpretation. EG, AS, CG, RC, DF, AC, and PZ: drafting the manuscript. All authors critically revised the manuscript, contributed to the article, and approved the submitted version.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
American Psychological Association (2014). The Road to Resilience. Washington, DC: American Psychological Association.
Bates, D., Mächler, M., Bolker, B., and Walker, S. (2015). Fitting Linear Mixed-Effects Models Using Lme4. J. Stat. Soft. 67, 1–48. doi: 10.18637/jss.v067.i01
Benedek, M., and Kaernbach, C. (2010). Decomposition of skin conductance data by means of nonnegative deconvolution. Psychophysiology 47, 647–658. doi: 10.1111/j.1469-8986.2009.00972.x
Bud, A. D. (2009). How do you feel–now? The anterior insula and human awareness. Nat. Rev. Neurosci. 10, 59–70. doi: 10.1038/nrn2555
Cassiello-Robbins, C., Anand, D., McMahon, K., Brout, J., Kelley, L., and Rosenthal, M. Z. (2021). A Preliminary Investigation of the Association Between Misophonia and Symptoms of Psychopathology and Personality Disorders. Front. Psychol. 11:519681. doi: 10.3389/fpsyg.2020.519681
Cassiello-Robbins, C., Anand, D., McMahon, K., Guetta, R., Trumbull, J., Kelley, L., et al. (2020). The Mediating Role of Emotion Regulation Within the Relationship Between Neuroticism and Misophonia: A Preliminary Investigation. Front. Psychiatry 11:847. doi: 10.3389/fpsyt.2020.00847
Cavanna, A. E., and Seri, S. (2015). Misophonia: Current perspectives Neuropsychiatr. Dis. Treat. 11, 2117–2123. doi: 10.2147/NDT.S81438
Clark, D. J., Chatterjee, S. A., McGuirk, T. E., Porges, E. C., Fox, E. J., and Balasubramanian, C. K. (2018). Sympathetic nervous system activity measured by skin conductance quantifies the challenge of walking adaptability tasks after stroke. Gait. Post. 60, 148–153. doi: 10.1016/j.gaitpost.2017.11.025
Core, R., and Team, R. (2019). A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing.
Critchley, H. D. (2005). Neural mechanisms of autonomic, affective, and cognitive integration. J. Comp. Neurol. 493, 154–166. doi: 10.1002/cne.20749
Daniels, E. C., Rodriguez, I. A., and Zabelina, D. L. (2020). Severity of misophonia symptoms is associated with worse cognitive control when exposed to misophonia trigger sounds. PLoS One 15:e0227118. doi: 10.1371/journal.pone.0227118
Derogatis, L. R. (1992). SCL-90-R: Administration, Scoring of Procedures Manual-II for the R (revised) Version and Other Instruments of the Psychopathology Rating Scale Series. Towson: Clinical Psychometric Research Incorporated.
Dibb, B., Golding, S. E., and Dozier, T. H. (2021). The development and validation of the Misophonia response scale. J. Psychosom Res. 149:110587. doi: 10.1016/j.jpsychores.2021.110587
Doberenz, S., Roth, W. T., Maslowski, N., Wollburg, E., and Kima, S. (2011). Methodological Considerations in Ambulatory Skin Conductance Monitoring. Int. J. Psychophysiol. 80, 87–95. doi: 10.1016/j.ijpsycho.2011.02.002
Edelstein, M., Brang, D., Rouw, R., and Ramachandran, V. S. (2013). Misophonia: Physiological investigations and case descriptions. Front. Hum. Neurosci. 7:296. doi: 10.3389/fnhum.2013.00296
Erfanian, M., Brout, J. J., Edelstein, M., Kumar, S., Mannino, M., Miller, L. J., et al. (2017). Investigating misophonia: A review of the literature, clinical implications and research agenda reflecting current neuroscience and emotion research perspectives. Eur. Psychiatry 41S:S681. doi: 10.1016/j.eurpsy.2017.01.1180
Erfanian, M., Kartsonaki, C., and Keshavarz, A. (2019). Misophonia and comorbid psychiatric symptoms: A preliminary study of clinical findings. Nord. J. Psychiatry 73, 219–228. doi: 10.1080/08039488.2019.1609086
Ferrer-Torres, A., and Giménez-Llort, L. (2022). Misophonia: A Systematic Review of Current and Future Trends in This Emerging Clinical Field. Int. J. Environ. Res. Public Health 19:6790. doi: 10.3390/ijerph19116790
Fitzmaurice, G. (2010). The Misophonia Activation Scale (MAS-1). Available online at: http://www.misophonia-uk.org/the-misophonia-activation-scale.html [accessed Nov 17, 2021].
Forte, G., Favieri, F., and Casagrande, M. (2019). Heart Rate Variability and Cognitive Function: A Systematic Review. Front. Neurosci. 13:710. doi: 10.3389/fnins.2019.00710
Friborg, O., Barlaug, D., Martinussen, M., Rosenvinge, J. H., and Hjemdal, O. (2005). Resilience in relation to personality and intelligence. Int. J. Methods Psychiatr. Res. 14, 29–40.
Gatti, E., Calzolari, E., Maggioni, E., and Obrist, M. (2018). Emotional ratings and skin conductance response to visual, auditory and haptic stimuli Sci. Data. 26:180120. doi: 10.1038/sdata.2018.120
Gernot, E. (2017). Hidden Signals—The History and Methods of Heart Rate Variability. Front. Public Health 5:265. doi: 10.3389/fpubh.2017.00265
Goodman, W. K., Price, L. H., Rasmussen, S. A., Mazure, C., Fleischmann, R. L., Hill, C. L., et al. (1989). The Yale-Brown Obsessive Compulsive Scale. I. Development, use, and reliability. Arch. Gen. Psychiatry 46, 1006–1011. doi: 10.1001/archpsyc.1989.01810110048007
Grabski, K., Lamalle, L., and Sato, M. (2012). Somatosensory-motor adaptation of orofacial actions in posterior parietal and ventral premotor cortices. PLoS One 7:e49117. doi: 10.1371/journal.pone.0049117
Hamilton, M. (1959). The assessment of anxiety states by rating. Br. J. Med. Psychol. 32, 50–55. doi: 10.1111/j.2044-8341.1959.tb00467.x
Hamilton, M. (1960). A rating scale for depression. J. Neurol. Neurosurg. Psychiatry. 23, 56–62. doi: 10.1136/jnnp.23.1.56
Huijbers, W., Pennartz, C. M., Cabeza, R., and Daselaar, S. M. (2011). The hippocampus is coupled with the default network during memory retrieval but not during memory encoding. PLoS One 6:e17463. doi: 10.1371/journal.pone.0017463
Jager, I., de Koning, P., Bost, T., Denys, D., and Vulink, N. (2020a). Misophonia: Phenomenology, comorbidity and demographics in a large sample. PLoS One 15:e0231390. doi: 10.1371/journal.pone.0231390
Jager, I. J., Vulink, N. C. C., Bergfeld, I. O., van Loon, A. J. J. M., and Denys, D. A. J. P. (2020b). Cognitive behavioral therapy for misophonia: A randomized clinical trial. Depress Anxiety 38, 708–718. doi: 10.1002/da.23127
Johnson, M. (2014). Misophonia Assessment Questionaire (MAQ). Available online at: https://misophoniatreatment.com/wp-content/uploads/2016/02/Binder_all_forms.pdf.
Kim, A. Y., Jang, E. H., Choi, K. W., Jeon, H. J., Byun, S., Sim, J. Y., et al. (2019). Skin conductance responses in Major Depressive Disorder (MDD) under mental arithmetic stress. PLoS One 14:e0213140. doi: 10.1371/journal.pone.0213140
Kumar, S., Dheerendra, P., Erfanian, M., Benzaquén, E., Sedley, W., Gander, P. E., et al. (2021). The Motor Basis for Misophonia. J. Neurosci. 41, 5762–5770. doi: 10.1523/JNEUROSCI.0261-21.2021
Kumar, S., Tansley-Hancock, O., Sedley, W., Winston, J. S., Callaghan, M. F., Allen, M., et al. (2017). The Brain Basis for Misophonia. Curr. Biol. 27, 527–533.
Liu, W., Lian, Z., and Liu, Y. (2008). Heart rate variability at different thermal comfort levels. Eur. J. Appl. Physiol. 103, 361–366. doi: 10.1007/s00421-008-0718-6
Mackersie, C. L., and Calderon-Moultrie, N. (2016). Autonomic Nervous System Reactivity During Speech Repetition Tasks: Heart Rate Variability and Skin Conductance. Ear. Hear. 37:118S–125S. doi: 10.1097/AUD.0000000000000305
McGeoch, P. D., and Rouw, R. (2020). How everyday sounds can trigger strong emotions: ASMR, misophonia and the feeling of wellbeing. Bioessays 42:e2000099. doi: 10.1002/bies.202000099
McKay, D., Kim, S. K., Mancusi, L., Storch, E. A., and Spankovich, C. (2018). Profile analysis of psychological symptoms associated with misophonia: A community sample. Behav. Ther. 49, 286–294. doi: 10.1016/j.beth.2017.07.002
McLachlan, R. S. (2009). A brief review of the anatomy and physiology of the limbic system. Can. J. Neurol. Sci. 2:S84–S87.
Motzkin, J. C., Philippi, C. L., Wolf, R. C., Baskaya, M. K., and Koenigs, M. (2015). Ventromedial prefrontal cortex is critical for the regulation of amygdala activity in humans. Biol. Psychiatry 77, 276–284. doi: 10.1016/j.biopsych.2014.02.014
Nagai, M., Kishi, K., and Kato, S. (2007). Insular cortex and neuropsychiatric disorders: A review of recent literature. Eur. Psychiatry 22, 387–394. doi: 10.1016/j.eurpsy.2007.02.006
Nandy, R., and Cordes, D. (2004). Improving the spatial specificity of canonical correlation analysis in fMRI. Magn. Reson. Med. 52, 947–952. doi: 10.1002/mrm.20234
Naylor, J., Caimino, C., Scutt, P., Hoare, D. J., and Baguley, D. M. (2021). The Prevalence and Severity of Misophonia in a UK Undergraduate Medical Student Population and Validation of the Amsterdam Misophonia Scale. Psychiatr. Q. 92, 609–619. doi: 10.1007/s11126-020-09825-3
Palumbo, D. B., Alsalman, O., De Ridder, D., Song, J. J., and Vanneste, S. (2018). Misophonia and Potential Underlying Mechanisms: A Perspective. Front. Psychol. 9:953. doi: 10.3389/fpsyg.2018.00953
Pickles, J. O. (2015). Auditory pathways: Anatomy and physiology. Handb. Clin. Neurol. 129, 3–25. doi: 10.1016/B978-0-444-62630-1.00001-9
Posada-Quintero, H. F., Kong, Y., Nguyen, K., Tran, C., Beardslee, L., Chen, L., et al. (2020). Using electrodermal activity to validate multilevel pain stimulation in healthy volunteers evoked by thermal grills. Am. J. Physiol. Regul. Integr. Comp. Physiol. 319:R366–R375. doi: 10.1152/ajpregu.00102.2020
Potgieter, I., MacDonald, C., Partridge, L., Cima, R., Sheldrake, J., and Hoare, D. J. (2019). Misophonia: A scoping review of research. J. Clin. Psychol. 75, 1203–1218. doi: 10.1002/jclp.22771
Raichle, M. E., MacLeod, A. M., Snyder, A. Z., Powers, W. J., Gusnard, D. A., and Shulman, G. L. (2001). A default mode of brain function. Proc. Natl. Acad. Sci. U.S.A. 98, 676–682.
Richter, A., Krämer, B., Diekhof, E. K., and Gruber, O. (2019). Resilience to adversity is associated with increased activity and connectivity in the VTA and hippocampus. Neuroimage. Clin. 23:101920. doi: 10.1016/j.nicl.2019.101920
Rolls, E. T. (2019). The cingulate cortex and limbic systems for action, emotion, and memory. Handb. Clin. Neurol. 166, 23–37. doi: 10.1016/B978-0-444-64196-0.00002-9
Sanchez, T. G., and da Silva, F. E. (2018). Familial misophonia or selective sound sensitivity syndrome: Evidence for autosomal dominant inheritance? Braz. J. Otorhinolaryngol. 84, 553–559. doi: 10.1016/j.bjorl.2017.06.014
Schröder, A., van Wingen, G., Eijsker, N., San Giorgi, R., Vulink, N. C., Turbyne, C., et al. (2019). Misophonia is associated with altered brain activity in the auditory cortex and salience network. Sci. Rep. 9:7542. doi: 10.1038/s41598-019-44084-8
Schröder, A., van Wingen, G., Vulink, N. C., and Denys, D. (2017). Commentary: The Brain Basis for Misophonia. Front. Behav. Neurosci. 11:111. doi: 10.3389/fnbeh.2017.00111
Schröder, A., Vulink, N., and Denys, D. (2013). Misophonia: Diagnostic criteria for a new psychiatric disorder. PLoS One 8:e54706. doi: 10.1371/journal.pone.0054706
Schwartz, G. E., Weinberger, D. A., and Singer, J. A. (1981). Cardiovascular differentiation of happiness, sadness, anger, and fear following imagery and exercise. Psychosom. Med. 43, 343–364. doi: 10.1097/00006842-198108000-00007
Seeley, W. (2019). The salience network: A neural system for perceiving and responding to homeostatic demands. J. Neurosci. 39, 9878–9882. doi: 10.1523/JNEUROSCI.1138-17.2019
Siepsiak, M., Sobczak, A. M., Bohaterewicz, B., Cichocki, L., and Dragan, W. L. (2020). Prevalence of Misophonia and Correlates of Its Symptoms among Inpatients with Depression. Int. J. Environ. Res. Public Health 15:5464. doi: 10.3390/ijerph17155464
Spielberger, C. D. (1989). State-Trait Anxiety Inventory: Bibliography, 2nd Edn. Palo Alto, CA: Consulting Psychologists Press.
SPM12 Software (2020). Statistical Parametric Mapping, updated at 1st October 2014. Available online at: https://www.fil.ion.ucl.ac.uk/spm/software/spm12/
Staib, M., Castegnetti, G., and Bach, D. R. (2015). Optimising a model-based approach to inferring fear learning from skin conductance responses. J. Neurosci. Methods 255, 131–138. doi: 10.1016/j.jneumeth.2015.08.009
Strigo, I. A., and Craig, A. D. (2016). Interoception, homeostatic emotions and sympathovagal balance. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 371:20160010. doi: 10.1098/rstb.2016.0010
Sweedo, S. E., Baguley, D. M., Denys, D., Dixon, L. J., Erfanian, M., Fioretti, A., et al. (2022). Consensus Definition of Misophonia: A Delphi Study. Front. Neurosci. 16:841816. doi: 10.3389/fnins.2022.841816
Tan, Y., Wei, D., Zhang, M., Yang, Y., Jelinèiæ, V., and Qiu, J. J. (2018). The role of mid-insula in the relationship between cardiac interoceptive attention and anxiety: Evidence from an fMRI study. Sci. Rep. 8:17280. doi: 10.1038/s41598-018-35635-6
Tarvainen, M. P., Niskanen, J. P., Lipponen, J. A., Ranta-Aho, P. O., and Karjalainen, P. A. (2014). Kubios HRV–heart rate variability analysis software. Comput. Methods Programs. Biomed. 113, 210–220. doi: 10.1016/j.cmpb.2013.07.024
Taylor, S., Conelea, C. A., McKay, D., Crowe, K. B., and Abramowitz, J. S. (2014). Sensory intolerance: Latent structure and psychopathologic correlates. Compr. Psychiatry 55, 1279–1284. doi: 10.1016/j.comppsych.2014.03.007
Thayer, J. F., and Brosschot, J. F. (2005). Psychosomatics and psychopathology: Looking up and down from the brain. Psychoneuroendocrinology 30, 1050–1058. doi: 10.1016/j.psyneuen.2005.04.014
Uddin, L. Q., Nomi, J. S., Hebert-Seropian, B., Ghaziri, J., and Boucher, O. (2017). Structure and function of the human insula J. Clin. Neurophysiol. 34, 300–306. doi: 10.1097/WNP.0000000000000377
Webber, T. A., Johnson, P. L., and Storch, E. A. (2014). Pediatric misophonia with comorbid obsessive-compulsive spectrum disorders. Gen. Hosp. Psychiatry 36:231.e1-2. doi: 10.1016/j.genhosppsych.2013.10.018
Wu, M. S., Lewin, A. B., Murphy, T. K., and Storch, E. A. (2014). Misophonia: Incidence, phenomenology, and clinical correlates in an undergraduate student sample. J. Clin. Psychol. 70, 994–1007. doi: 10.1002/jclp.22098
Keywords: anatomic pathways, autonomic system, functional magnetic resonance, psychometric assessment, trigger sounds
Citation: Grossini E, Stecco A, Gramaglia C, De Zanet D, Cantello R, Gori B, Negroni D, Azzolina D, Ferrante D, Feggi A, Carriero A and Zeppegno P (2022) Misophonia: Analysis of the neuroanatomic patterns at the basis of psychiatric symptoms and changes of the orthosympathetic/ parasympathetic balance. Front. Neurosci. 16:827998. doi: 10.3389/fnins.2022.827998
Received: 02 December 2021; Accepted: 14 July 2022;
Published: 11 August 2022.
Edited by:
Julia Campbell, The University of Texas at Austin, United StatesReviewed by:
Phillip Gilley, University of Colorado Boulder, United StatesSharon Hunter, University of Colorado Anschutz Medical Campus, United States
Copyright © 2022 Grossini, Stecco, Gramaglia, De Zanet, Cantello, Gori, Negroni, Azzolina, Ferrante, Feggi, Carriero and Zeppegno. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Elena Grossini ZWxlbmEuZ3Jvc3NpbmlAbWVkLnVuaXVwby5pdA==
†These authors have contributed equally to this work and share first authorship
‡These authors have contributed equally to this work and share first authorship