The habenula is an epithalamic structure that functionally links the forebrain with the midbrain structures that are involved in the release of dopamine (i.e., the substantia nigra pars compacta and the ventral tegmental area) and serotonin (i.e., raphe nucleus) (Wang and Aghajanian, 1977; Herkenham and Nauta, 1979; Hikosaka et al., 2008; Hong et al., 2011; Proulx et al., 2014; Hu et al., 2020). As a region that could influence both the dopaminergic and serotonergic systems, the habenula is thought to play a key role in not only sleep and wakefulness but also in regulating various emotional and cognitive functions. Animal studies showed that activities in lateral habenula (LHb) increased during the processing of aversive events such as omission of predicted rewards, and stimuli provoking anxiety, stress, pain and fear (Matsumoto and Hikosaka, 2007; Hikosaka, 2010; Yamaguchi et al., 2013; Hu et al., 2020).

Hyperexcitability and dysfunction of the LHb have been implicated in the development of psychiatric disorders including depressive disorder and bipolar disorders (Fakhoury, 2017; Yang et al., 2018b). In rodents, LHb firing rate and metabolism is elevated in parallel with depressive-like phenotypes such as reduction in locomotor and rearing behaviors (Caldecott-Hazard et al., 1988). Habenula activities also increase during acquisition and recall of conditioned fear (González-Pardo et al., 2012). In human participants, high-resolution magnetic resonance imaging has revealed smaller habenula volume in patients with depressive and bipolar disorders (Savitz et al., 2011a). Dysfunction of the LHb has also been involved in different cognitive disorders, such as schizophrenia (Shepard et al., 2006) and addiction (Velasquez et al., 2014). More direct evidence of the involvement of the LHb in psychiatric disorders in humans comes from deep brain stimulation (DBS) of the LHb that has potential therapeutic effects in treatment-resistant depression, bipolar disorder, and schizophrenia (Sartorius et al., 2010; Zhang et al., 2019; Wang et al., 2020). However, how habenula activities are modulated during the processing of emotional information in humans is still poorly understood.

The processing of emotional information is crucial for an individual’s mental health and has a substantial influence on social interactions and different cognitive processes. Dysfunction and dysregulation of emotion-related brain circuits may precipitate mood disorders (Phillips et al., 2003b). Investigating the neural activities in response to emotional stimuli in the cortical-habenula network is crucial to our understanding of emotional information processing in the brain. This might also shed light on how to modulate habenula in the treatment of psychiatric disorders. In this study, we utilize the unique opportunity offered by DBS surgery targeting habenula as a potential treatment for psychiatric disorders. We measured local field potentials (LFPs) from the habenula area using the electrodes implanted for DBS in patients during a passive emotional picture-viewing task (Figure 1; Materials and methods). Whole-brain magnetoencephalography (MEG) was simultaneously recorded. This allowed us to investigate changes in the habenula neural activity and its functional connectivity with cortical areas induced by the stimuli of different emotional valence. The high temporal resolution of the LFP and MEG measurements also allowed us to evaluate how local activities and cross-region connectivity change over time in the processing of emotional stimuli. Previous studies on rodent models of depression showed that, during the depression-like state in rodents, LHb neuron firing increased with the mean firing rate at the theta band (Li et al., 2011) and LHb neurons fire in bursts and phase locked to local theta band field potentials (Yang et al., 2018a). Therefore, we hypothesize that theta band activity in the habenula LFPs in humans would increase in response to negative emotional stimuli.

Figure 1

Download asset Open asset

This study has showed that neural activities in the theta/alpha frequency band within the habenula and prefrontal cortical regions, as well as the connectivity between these structures in the same frequency band, are modulated in an emotional picture-viewing task in human participants. Compared with positive emotional stimuli, negative emotional stimuli were associated with higher transient increase in theta/alpha activity in both habenula and bilateral frontal cortex with a short latency (from 100 to 500 ms) after stimulus onset. Furthermore, higher theta/alpha coherence between habenula and right PFC was observed at 800–1300 ms after the stimulus onset, which was correlated with another increase in theta power in the habenula with a long latency (from 2700 to 3300 ms) after stimulus onset. These changes correlated with the emotional valence but not with the stimulus arousal of the presented figures. These activity changes at different time windows may reflect the different neuropsychological processes underlying emotion perception including identification and appraisal of emotional material, production of affective states, and autonomic response regulation and recovery (Phillips et al., 2003a). The later effects of increased theta activities in the habenula when the stimuli disappeared were also supported by other literature showing that there can be prolonged effects of negative stimuli in the neural structure involved in emotional processing (Haas et al., 2008; Puccetti et al., 2021). In particular, greater sustained patterns of brain activity in the medial PFC when responding to blocks of negative facial expressions were associated with higher scores of neuroticism across participants (Haas et al., 2008). Slower amygdala recovery from negative images also predicts greater trait neuroticism, lower levels of likability of a set of social stimuli (neutral faces), and declined day-to-day psychological well-being (Schuyler et al., 2014; Puccetti et al., 2021). This is the first study, to our knowledge, implicating increased theta band activities in the habenula-PFC network in negative emotions in human patients.

All data generated or analysed during this study are included in the manuscript and supporting files.

1. 1. Avants BB
   2. Epstein CL
   3. Grossman M
   4. Gee JC

   (2008) Symmetric diffeomorphic image registration with cross-correlation: evaluating automated labeling of elderly and neurodegenerative brain

   Medical Image Analysis 12:26–41.

   https://doi.org/10.1016/j.media.2007.06.004

   • PubMed
   • Google Scholar
2. 1. Bai L
   2. Ma H
   3. Huang Y-X

   (2005) The development of native chinese affective picture system-a pretest in 46 college students

   Chinese Mental Health J 11:719–722.

   https://doi.org/10.3321/J.ISSN:1000-6729.2005.11.001

   • Google Scholar
3. 1. Baldermann JC
   2. Melzer C
   3. Zapf A
   4. Kohl S
   5. Timmermann L
   6. Tittgemeyer M
   7. Huys D
   8. Visser-Vandewalle V
   9. Kühn AA
   10. Horn A
   11. Kuhn J

   (2019) Connectivity profile predictive of effective deep brain stimulation in Obsessive-Compulsive disorder

   Biological Psychiatry 85:735–743.

   https://doi.org/10.1016/j.biopsych.2018.12.019

   • PubMed
   • Google Scholar
4. 1. Beraha E
   2. Eggers J
   3. Hindi Attar C
   4. Gutwinski S
   5. Schlagenhauf F
   6. Stoy M
   7. Sterzer P
   8. Kienast T
   9. Heinz A
   10. Bermpohl F

   (2012) Hemispheric asymmetry for affective stimulus processing in healthy subjects--a fMRI study

   PLOS ONE 7:e46931.

   https://doi.org/10.1371/journal.pone.0046931

   • PubMed
   • Google Scholar
5. 1. Boraud T
   2. Bezard E
   3. Bioulac B
   4. Gross C

   (1996) High frequency stimulation of the internal globus pallidus (GPi) simultaneously improves parkinsonian symptoms and reduces the firing frequency of GPi neurons in the MPTP-treated monkey

   Neuroscience Letters 215:17–20.

   https://doi.org/10.1016/S0304-3940(96)12943-8

   • PubMed
   • Google Scholar
6. 1. Bromberg-Martin ES
   2. Hikosaka O

   (2011) Lateral habenula neurons signal errors in the prediction of reward information

   Nature Neuroscience 14:1209–1216.

   https://doi.org/10.1038/nn.2902

   • PubMed
   • Google Scholar
7. 1. Caldecott-Hazard S
   2. Mazziotta J
   3. Phelps M

   (1988) Cerebral correlates of depressed behavior in rats, visualized using 14C-2-deoxyglucose autoradiography

   The Journal of Neuroscience 8:1951–1961.

   https://doi.org/10.1523/JNEUROSCI.08-06-01951.1988

   • PubMed
   • Google Scholar
8. 1. Canli T
   2. Desmond JE
   3. Zhao Z
   4. Glover G
   5. Gabrieli JD

   (1998) Hemispheric asymmetry for emotional stimuli detected with fMRI

   NeuroReport 9:3233–3239.

   https://doi.org/10.1097/00001756-199810050-00019

   • PubMed
   • Google Scholar
9. 1. Chiken S
   2. Nambu A

   (2016) Mechanism of deep brain stimulation: inhibition, excitation, or disruption?

   The Neuroscientist : A Review Journal Bringing Neurobiology, Neurology and Psychiatry 22:313–322.

   https://doi.org/10.1177/1073858415581986

   • PubMed
   • Google Scholar
10. Software

    1. Crosse MJ
    2. Foxe JJ
    3. Molholm S

    (2020) PERMUTOOLS: A MATLAB Package for Multivariate Permutation Testing, version 9f4c6c8

    Github.

    https://github.com/mickcrosse/PERMUTOOLS
11. 1. Cui Y
    2. Yang Y
    3. Ni Z
    4. Dong Y
    5. Cai G
    6. Foncelle A
    7. Ma S
    8. Sang K
    9. Tang S
    10. Li Y
    11. Shen Y
    12. Berry H
    13. Wu S
    14. Hu H

    (2018) Astroglial Kir4.1 in the lateral habenula drives neuronal bursts in depression

    Nature 554:323–327.

    https://doi.org/10.1038/nature25752

    • PubMed
    • Google Scholar
12. 1. Davidson RJ

    (1992) Anterior cerebral asymmetry and the nature of emotion

    Brain and Cognition 20:125–151.

    https://doi.org/10.1016/0278-2626(92)90065-T

    • PubMed
    • Google Scholar
13. 1. Donoghue T
    2. Haller M
    3. Peterson EJ
    4. Varma P
    5. Sebastian P
    6. Gao R
    7. Noto T
    8. Lara AH
    9. Wallis JD
    10. Knight RT
    11. Shestyuk A
    12. Voytek B

    (2020) Parameterizing neural power spectra into periodic and aperiodic components

    Nature Neuroscience 23:1655–1665.

    https://doi.org/10.1038/s41593-020-00744-x

    • PubMed
    • Google Scholar
14. 1. Fakhoury M

    (2017) The habenula in psychiatric disorders: more than three decades of translational investigation

    Neuroscience & Biobehavioral Reviews 83:721–735.

    https://doi.org/10.1016/j.neubiorev.2017.02.010

    • PubMed
    • Google Scholar
15. 1. Fonov V
    2. Evans AC
    3. Botteron K
    4. Almli CR
    5. McKinstry RC
    6. Collins DL

    (2011) Unbiased average age-appropriate atlases for pediatric studies

    NeuroImage 54:313–327.

    https://doi.org/10.1016/j.neuroimage.2010.07.033

    • Google Scholar
16. 1. Friedman A
    2. Lax E
    3. Dikshtein Y
    4. Abraham L
    5. Flaumenhaft Y
    6. Sudai E
    7. Ben-Tzion M
    8. Yadid G

    (2011) Electrical stimulation of the lateral habenula produces an inhibitory effect on sucrose self-administration

    Neuropharmacology 60:381–387.

    https://doi.org/10.1016/j.neuropharm.2010.10.006

    • Google Scholar
17. 1. Gainotti G

    (2012) Unconscious processing of emotions and the right hemisphere

    Neuropsychologia 50:205–218.

    https://doi.org/10.1016/j.neuropsychologia.2011.12.005

    • Google Scholar
18. 1. González-Pardo H
    2. Conejo NM
    3. Lana G
    4. Arias JL

    (2012) Different brain networks underlying the acquisition and expression of contextual fear conditioning: a metabolic mapping study

    Neuroscience 202:234–242.

    https://doi.org/10.1016/j.neuroscience.2011.11.064

    • Google Scholar
19. 1. Gross J
    2. Kujala J
    3. Hamalainen M
    4. Timmermann L
    5. Schnitzler A
    6. Salmelin R

    (2001) Dynamic imaging of coherent sources: Studying neural interactions in the human brain

    PNAS 98:694–699.

    https://doi.org/10.1073/pnas.98.2.694

    • Google Scholar
20. 1. Haas BW
    2. Constable RT
    3. Canli T

    (2008) Stop the sadness: Neuroticism is associated with sustained medial prefrontal cortex response to emotional facial expressions

    NeuroImage 42:385–392.

    https://doi.org/10.1016/j.neuroimage.2008.04.027

    • Google Scholar
21. Preprint

    1. Haller M
    2. Donoghue T
    3. Peterson E
    4. Varma P
    5. Sebastian P
    6. Gao R

    (2018) Parameterizing neural power spectra

    bioRxiv.

    https://doi.org/10.1101/299859

    • Google Scholar
22. 1. Herkenham M
    2. Nauta WJH

    (1979) Efferent connections of the habenular nuclei in the rat

    The Journal of Comparative Neurology 187:19–47.

    https://doi.org/10.1002/cne.901870103

    • Google Scholar
23. 1. Hikosaka O
    2. Sesack SR
    3. Lecourtier L
    4. Shepard PD

    (2008) Habenula: crossroad between the basal ganglia and the limbic system

    Journal of Neuroscience 28:11825–11829.

    https://doi.org/10.1523/JNEUROSCI.3463-08.2008

    • PubMed
    • Google Scholar
24. 1. Hikosaka O

    (2010) The habenula: from stress evasion to value-based decision-making

    Nature Reviews Neuroscience 11:503–513.

    https://doi.org/10.1038/nrn2866

    • Google Scholar
25. 1. Hirschmann J
    2. Özkurt TE
    3. Butz M
    4. Homburger M
    5. Elben S
    6. Hartmann CJ
    7. Vesper J
    8. Wojtecki L
    9. Schnitzler A

    (2011) Distinct oscillatory STN-cortical loops revealed by simultaneous MEG and local field potential recordings in patients with Parkinson's disease

    NeuroImage 55:1159–1168.

    https://doi.org/10.1016/j.neuroimage.2010.11.063

    • PubMed
    • Google Scholar
26. 1. Hong S
    2. Jhou TC
    3. Smith M
    4. Saleem KS
    5. Hikosaka O

    (2011) Negative reward signals from the lateral habenula to dopamine neurons are mediated by rostromedial tegmental nucleus in primates

    Journal of Neuroscience 31:11457–11471.

    https://doi.org/10.1523/JNEUROSCI.1384-11.2011

    • Google Scholar
27. 1. Horn A
    2. Reich M
    3. Vorwerk J
    4. Li N
    5. Wenzel G
    6. Fang Q
    7. Schmitz-Hübsch T
    8. Nickl R
    9. Kupsch A
    10. Volkmann J
    11. Kühn AA
    12. Fox MD

    (2017) Connectivity predicts deep brain stimulation outcome in parkinson disease

    Annals of Neurology 82:67–78.

    https://doi.org/10.1002/ana.24974

    • PubMed
    • Google Scholar
28. 1. Horn A
    2. Li N
    3. Dembek TA
    4. Kappel A
    5. Boulay C
    6. Ewert S
    7. Tietze A
    8. Husch A
    9. Perera T
    10. Neumann WJ
    11. Reisert M
    12. Si H
    13. Oostenveld R
    14. Rorden C
    15. Yeh FC
    16. Fang Q
    17. Herrington TM
    18. Vorwerk J
    19. Kühn AA

    (2019) Lead-DBS v2: towards a comprehensive pipeline for deep brain stimulation imaging

    NeuroImage 184:293–316.

    https://doi.org/10.1016/j.neuroimage.2018.08.068

    • PubMed
    • Google Scholar
29. 1. Horn A
    2. Kühn AA

    (2015) Lead-DBS: a toolbox for deep brain stimulation electrode localizations and visualizations

    NeuroImage 107:127–135.

    https://doi.org/10.1016/j.neuroimage.2014.12.002

    • PubMed
    • Google Scholar
30. 1. Hu H
    2. Cui Y
    3. Yang Y

    (2020) Circuits and functions of the lateral habenula in health and in disease

    Nature Reviews Neuroscience 21:277–295.

    https://doi.org/10.1038/s41583-020-0292-4

    • Google Scholar
31. 1. Huebl J
    2. Spitzer B
    3. Brücke C
    4. Schönecker T
    5. Kupsch A
    6. Alesch F
    7. Schneider G-H
    8. Kühn AA

    (2014) Oscillatory subthalamic nucleus activity is modulated by dopamine during emotional processing in Parkinson's disease

    Cortex 60:69–81.

    https://doi.org/10.1016/j.cortex.2014.02.019

    • Google Scholar
32. 1. Huebl J
    2. Brücke C
    3. Merkl A
    4. Bajbouj M
    5. Schneider G-H
    6. Kühn AA

    (2016) Processing of emotional stimuli is reflected by modulations of beta band activity in the subgenual anterior cingulate cortex in patients with treatment resistant depression

    Social Cognitive and Affective Neuroscience 11:1290–1298.

    https://doi.org/10.1093/scan/nsw038

    • Google Scholar
33. 1. Husch A
    2. V. Petersen M
    3. Gemmar P
    4. Goncalves J
    5. Hertel F

    (2018) PaCER - A fully automated method for electrode trajectory and contact reconstruction in deep brain stimulation

    NeuroImage: Clinical 17:80–89.

    https://doi.org/10.1016/j.nicl.2017.10.004

    • Google Scholar
34. 1. Jakobs M
    2. Pitzer C
    3. Sartorius A
    4. Unterberg A
    5. Kiening K

    (2019) Acute 5-Hz deep brain stimulation of the lateral habenula is associated with depressive-like behavior in male wild-type Wistar rats

    Brain Research 1721:146283.

    https://doi.org/10.1016/j.brainres.2019.06.002

    • PubMed
    • Google Scholar
35. 1. Ji H
    2. Shepard PD

    (2007) Lateral Habenula Stimulation Inhibits Rat Midbrain Dopamine Neurons through a GABAA Receptor-Mediated Mechanism

    Journal of Neuroscience 27:6923–6930.

    https://doi.org/10.1523/JNEUROSCI.0958-07.2007

    • Google Scholar
36. 1. Kuhn AA
    2. Hariz MI
    3. Silberstein P
    4. Tisch S
    5. Kupsch A
    6. Schneider G-H
    7. Limousin-Dowsey P
    8. Yarrow K
    9. Brown P

    (2005) Activation of the subthalamic region during emotional processing in Parkinson disease

    Neurology 65:707–713.

    https://doi.org/10.1212/01.wnl.0000174438.78399.bc

    • Google Scholar
37. 1. Kuhn AA
    2. Kempf F
    3. Brucke C
    4. Gaynor Doyle L
    5. Martinez-Torres I
    6. Pogosyan A
    7. Trottenberg T
    8. Kupsch A
    9. Schneider G-H
    10. Hariz MI
    11. Vandenberghe W
    12. Nuttin B
    13. Brown P

    (2008) High-Frequency Stimulation of the Subthalamic Nucleus Suppresses Oscillatory Activity in Patients with Parkinson's Disease in Parallel with Improvement in Motor Performance

    Journal of Neuroscience 28:6165–6173.

    https://doi.org/10.1523/JNEUROSCI.0282-08.2008

    • Google Scholar
38. 1. Lawson RP
    2. Drevets WC
    3. Roiser JP

    (2013) Defining the habenula in human neuroimaging studies

    NeuroImage 64:722–727.

    https://doi.org/10.1016/j.neuroimage.2012.08.076

    • Google Scholar
39. 1. Lawson RP
    2. Nord CL
    3. Seymour B
    4. Thomas DL
    5. Dayan P
    6. Pilling S
    7. Roiser JP

    (2017) Disrupted habenula function in major depression

    Molecular Psychiatry 22:202–208.

    https://doi.org/10.1038/mp.2016.81

    • Google Scholar
40. 1. Lecca S
    2. Pelosi A
    3. Tchenio A
    4. Moutkine I
    5. Lujan R
    6. Hervé D
    7. Mameli M

    (2016) Rescue of GABAB and GIRK function in the lateral habenula by protein phosphatase 2A inhibition ameliorates depression-like phenotypes in mice

    Nature Medicine 22:254–261.

    https://doi.org/10.1038/nm.4037

    • Google Scholar
41. 1. Li B
    2. Piriz J
    3. Mirrione M
    4. Chung C
    5. Proulx CD
    6. Schulz D
    7. Henn F
    8. Malinow R

    (2011) Synaptic potentiation onto habenula neurons in the learned helplessness model of depression

    Nature 470:535–539.

    https://doi.org/10.1038/nature09742

    • Google Scholar
42. 1. Li K
    2. Zhou T
    3. Liao L
    4. Yang Z
    5. Wong C
    6. Henn F
    7. Malinow R
    8. Yates JR
    9. Hu H

    (2013) βCaMKII in Lateral Habenula Mediates Core Symptoms of Depression

    Science 341:1016–1020.

    https://doi.org/10.1126/science.1240729

    • Google Scholar
43. 1. Litvak V
    2. Eusebio A
    3. Jha A
    4. Oostenveld R
    5. Barnes GR
    6. Penny WD
    7. Zrinzo L
    8. Hariz MI
    9. Limousin P
    10. Friston KJ
    11. Brown P

    (2010) Optimized beamforming for simultaneous MEG and intracranial local field potential recordings in deep brain stimulation patients

    NeuroImage 50:1578–1588.

    https://doi.org/10.1016/j.neuroimage.2009.12.115

    • Google Scholar
44. 1. Litvak V
    2. Jha A
    3. Eusebio A
    4. Oostenveld R
    5. Foltynie T
    6. Limousin P
    7. Zrinzo L
    8. Hariz MI
    9. Friston K
    10. Brown P

    (2011) Resting oscillatory cortico-subthalamic connectivity in patients with Parkinson’s disease

    Brain 134:359–374.

    https://doi.org/10.1093/brain/awq332

    • Google Scholar
45. 1. Maris E
    2. Oostenveld R

    (2007) Nonparametric statistical testing of EEG- and MEG-data

    Journal of Neuroscience Methods 164:177–190.

    https://doi.org/10.1016/j.jneumeth.2007.03.024

    • Google Scholar
46. 1. Marmor O
    2. Valsky D
    3. Joshua M
    4. Bick AS
    5. Arkadir D
    6. Tamir I
    7. Bergman H
    8. Israel Z
    9. Eitan R

    (2017) Local vs. volume conductance activity of field potentials in the human subthalamic nucleus

    Journal of Neurophysiology 117:2140–2151.

    https://doi.org/10.1152/jn.00756.2016

    • Google Scholar
47. 1. Matsumoto M
    2. Hikosaka O

    (2007) Lateral habenula as a source of negative reward signals in dopamine neurons

    Nature 447:1111–1115.

    https://doi.org/10.1038/nature05860

    • Google Scholar
48. 1. Matsumoto M
    2. Hikosaka O

    (2009) Representation of negative motivational value in the primate lateral habenula

    Nature Neuroscience 12:77–84.

    https://doi.org/10.1038/nn.2233

    • Google Scholar
49. 1. Morris JS
    2. Smith KA
    3. Cowen PJ
    4. Friston KJ
    5. Dolan RJ

    (1999) Covariation of Activity in Habenula and Dorsal Raphé Nuclei Following Tryptophan Depletion

    NeuroImage 10:163–172.

    https://doi.org/10.1006/nimg.1999.0455

    • Google Scholar
50. 1. Neumann W-J
    2. Jha A
    3. Bock A
    4. Huebl J
    5. Horn A
    6. Schneider G-H
    7. Sander TH
    8. Litvak V
    9. Kühn AA

    (2015) Cortico-pallidal oscillatory connectivity in patients with dystonia

    Brain 138:1894–1906.

    https://doi.org/10.1093/brain/awv109

    • Google Scholar
51. 1. Oostenveld R
    2. Fries P
    3. Maris E
    4. Schoffelen JM

    (2011) FieldTrip: open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data

    Computational Intelligence and Neuroscience 2011:1–9.

    https://doi.org/10.1155/2011/156869

    • PubMed
    • Google Scholar
52. 1. Oswal A
    2. Beudel M
    3. Zrinzo L
    4. Limousin P
    5. Hariz M
    6. Foltynie T
    7. Litvak V
    8. Brown P

    (2016) Deep brain stimulation modulates synchrony within spatially and spectrally distinct resting state networks in Parkinson’s disease

    Brain 139:1482–1496.

    https://doi.org/10.1093/brain/aww048

    • Google Scholar
53. 1. Pascual-Leone A
    2. Catala MD
    3. Pascual AP-L
    4. Pascual-Leone Pascual A

    (1996) Lateralized effect of rapid-rate transcranial magnetic stimulation of the prefrontal cortex on mood

    Neurology 46:499–502.

    https://doi.org/10.1212/WNL.46.2.499

    • Google Scholar
54. 1. Phillips ML
    2. Drevets WC
    3. Rauch SL
    4. Lane R

    (2003a) Neurobiology of emotion perception I: the neural basis of normal emotion perception

    Biological Psychiatry 54:504–514.

    https://doi.org/10.1016/S0006-3223(03)00168-9

    • Google Scholar
55. 1. Phillips ML
    2. Drevets WC
    3. Rauch SL
    4. Lane R

    (2003b) Neurobiology of emotion perception II: implications for major psychiatric disorders

    Biological Psychiatry 54:515–528.

    https://doi.org/10.1016/S0006-3223(03)00171-9

    • Google Scholar
56. 1. Proulx CD
    2. Hikosaka O
    3. Malinow R

    (2014) Reward processing by the lateral habenula in normal and depressive behaviors

    Nature Neuroscience 17:1146–1152.

    https://doi.org/10.1038/nn.3779

    • Google Scholar
57. 1. Puccetti NA
    2. Schaefer SM
    3. van Reekum CM
    4. Ong AD
    5. Almeida DM
    6. Ryff CD
    7. Davidson RJ
    8. Heller AS

    (2021) Linking amygdala persistence to Real-World emotional experience and psychological Well-Being

    The Journal of Neuroscience 41:3721–3730.

    https://doi.org/10.1523/JNEUROSCI.1637-20.2021

    • PubMed
    • Google Scholar
58. 1. Riva-Posse P
    2. Choi KS
    3. Holtzheimer PE
    4. McIntyre CC
    5. Gross RE
    6. Chaturvedi A
    7. Crowell AL
    8. Garlow SJ
    9. Rajendra JK
    10. Mayberg HS

    (2014) Defining critical white matter pathways mediating successful subcallosal cingulate deep brain stimulation for treatment-resistant depression

    Biological Psychiatry 76:963–969.

    https://doi.org/10.1016/j.biopsych.2014.03.029

    • Google Scholar
59. 1. Sartorius A
    2. Kiening KL
    3. Kirsch P
    4. von Gall CC
    5. Haberkorn U
    6. Unterberg AW
    7. Henn FA
    8. Meyer-Lindenberg A

    (2010) Remission of major depression under deep brain stimulation of the lateral habenula in a therapy-refractory patient

    Biological Psychiatry 67:e9–e11.

    https://doi.org/10.1016/j.biopsych.2009.08.027

    • Google Scholar
60. 1. Savitz JB
    2. Nugent AC
    3. Bogers W
    4. Roiser JP
    5. Bain EE
    6. Neumeister A
    7. Zarate CA
    8. Manji HK
    9. Cannon DM
    10. Marrett S
    11. Henn F
    12. Charney DS
    13. Drevets WC

    (2011a) Habenula Volume in Bipolar Disorder and Major Depressive Disorder: A High-Resolution Magnetic Resonance Imaging Study

    Biological Psychiatry 69:336–343.

    https://doi.org/10.1016/j.biopsych.2010.09.027

    • Google Scholar
61. 1. Savitz JB
    2. Bonne O
    3. Nugent AC
    4. Vythilingam M
    5. Bogers W
    6. Charney DS
    7. Drevets WC

    (2011b) Habenula volume in post-traumatic stress disorder measured with high-resolution MRI

    Biology of Mood & Anxiety Disorders 1:7.

    https://doi.org/10.1186/2045-5380-1-7

    • Google Scholar
62. 1. Schuyler BS
    2. Kral TRA
    3. Jacquart J
    4. Burghy CA
    5. Weng HY
    6. Perlman DM
    7. Bachhuber DRW
    8. Rosenkranz MA
    9. MacCoon DG
    10. van Reekum CM
    11. Lutz A
    12. Davidson RJ

    (2014) Temporal dynamics of emotional responding: amygdala recovery predicts emotional traits

    Social Cognitive and Affective Neuroscience 9:176–181.

    https://doi.org/10.1093/scan/nss131

    • Google Scholar
63. 1. Shepard PD
    2. Holcomb HH
    3. Gold JM

    (2006) Schizophrenia in translation: the presence of absence: habenular regulation of dopamine neurons and the encoding of negative outcomes

    Schizophrenia Bulletin 32:417–421.

    https://doi.org/10.1093/schbul/sbj083

    • Google Scholar
64. 1. Shumake J
    2. Ilango A
    3. Scheich H
    4. Wetzel W
    5. Ohl FW

    (2010) Differential neuromodulation of acquisition and retrieval of avoidance learning by the lateral habenula and ventral tegmental area

    Journal of Neuroscience 30:5876–5883.

    https://doi.org/10.1523/JNEUROSCI.3604-09.2010

    • Google Scholar
65. 1. Starkstein SE
    2. Robinson RG
    3. Price TR

    (1987) Comparison of cortical and subcortical lesions in the production of poststroke mood disorders

    Brain 110:1045–1059.

    https://doi.org/10.1093/brain/110.4.1045

    • PubMed
    • Google Scholar
66. 1. Taulu S
    2. Simola J

    (2006) Spatiotemporal signal space separation method for rejecting nearby interference in MEG measurements

    Physics in Medicine and Biology 51:1759–1768.

    https://doi.org/10.1088/0031-9155/51/7/008

    • Google Scholar
67. 1. Vadovičová K

    (2014) Affective and cognitive prefrontal cortex projections to the lateral habenula in humans

    Frontiers in Human Neuroscience 8:819.

    https://doi.org/10.3389/fnhum.2014.00819

    • PubMed
    • Google Scholar
68. 1. Velasquez KM
    2. Molfese DL
    3. Salas R

    (2014) The role of the habenula in drug addiction

    Frontiers in Human Neuroscience 8:174.

    https://doi.org/10.3389/fnhum.2014.00174

    • PubMed
    • Google Scholar
69. 1. Wang D
    2. Li Y
    3. Feng Q
    4. Guo Q
    5. Zhou J
    6. Luo M

    (2017) Learning shapes the aversion and reward responses of lateral habenula neurons

    eLife 6:e23045.

    https://doi.org/10.7554/eLife.23045

    • PubMed
    • Google Scholar
70. 1. Wang Y
    2. Zhang C
    3. Zhang Y
    4. Gong H
    5. Li J
    6. Jin H
    7. Li D
    8. Liu D
    9. Sun B

    (2020) Habenula deep brain stimulation for intractable schizophrenia: a pilot study

    Neurosurgical Focus 49:E9.

    https://doi.org/10.3171/2020.4.FOCUS20174

    • Google Scholar
71. 1. Wang R
    2. Aghajanian G

    (1977) Physiological evidence for habenula as major link between forebrain and midbrain raphe

    Science 197:89–91.

    https://doi.org/10.1126/science.194312

    • Google Scholar
72. 1. Welter M-L
    2. Houeto J-L
    3. Bonnet A-M
    4. Bejjani P-B
    5. Mesnage V
    6. Dormont D
    7. Navarro S
    8. Cornu P
    9. Agid Y
    10. Pidoux B

    (2004) Effects of high-frequency stimulation on subthalamic neuronal activity in parkinsonian patients

    Archives of Neurology 61:89–96.

    https://doi.org/10.1001/archneur.61.1.89

    • Google Scholar
73. 1. Wyczesany M
    2. Capotosto P
    3. Zappasodi F
    4. Prete G

    (2018) Hemispheric asymmetries and emotions: Evidence from effective connectivity

    Neuropsychologia 121:98–105.

    https://doi.org/10.1016/j.neuropsychologia.2018.10.007

    • Google Scholar
74. 1. Yamaguchi T
    2. Danjo T
    3. Pastan I
    4. Hikida T
    5. Nakanishi S

    (2013) Distinct Roles of Segregated Transmission of the Septo-Habenular Pathway in Anxiety and Fear

    Neuron 78:537–544.

    https://doi.org/10.1016/j.neuron.2013.02.035

    • Google Scholar
75. 1. Yang Y
    2. Cui Y
    3. Sang K
    4. Dong Y
    5. Ni Z
    6. Ma S
    7. Hu H

    (2018a) Ketamine blocks bursting in the lateral habenula to rapidly relieve depression

    Nature 554:317–322.

    https://doi.org/10.1038/nature25509

    • PubMed
    • Google Scholar
76. 1. Yang Y
    2. Wang H
    3. Hu J
    4. Hu H

    (2018b) Lateral habenula in the pathophysiology of depression

    Current Opinion in Neurobiology 48:90–96.

    https://doi.org/10.1016/j.conb.2017.10.024

    • Google Scholar
77. 1. Zhang C
    2. Kim S-G
    3. Li D
    4. Zhang Y
    5. Li Y
    6. Husch A
    7. Hertel F
    8. Yan F
    9. Voon V
    10. Sun B

    (2019) Habenula deep brain stimulation for refractory bipolar disorder

    Brain Stimulation 12:1298–1300.

    https://doi.org/10.1016/j.brs.2019.05.010

    • Google Scholar

Acceptance summary:

The findings provide novel insight into the network dynamics associated with the processing of emotional stimuli and in particular the role of the habenula. The key finding was a transient increase in habenula theta/alpha activity for negative compared to positive stimuli. Furthermore, there was a later increase in habenula-prefrontal oscillatory coupling in the same band. These are important data, as there are few reports of direct recordings from the habenula together with the MEG in humans performing cognitive tasks.

Decision letter after peer review:

Thank you for submitting your article "Increased theta/alpha synchrony in the habenula-prefrontal network with negative emotional stimuli in human patients" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, including Ole Jensen as the Reviewing Editor and Reviewer #1, and the evaluation has been overseen by Floris de Lange as the Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity: Jan Hirschmann (Reviewer #2); Diego Vidaurre (Reviewer #3).

The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.

Essential Revisions:

Please pay special attention to

1) Concerns on volume conduction

2) The issue of evoked versus oscillatory responses

3) Concerns on low-level features of the visual stimuli driving the evoked effects

4) Robustness of findings in the light of saccadic artefacts and statistical assessment

Reviewer #2 (Recommendations for the authors):

Were the stimuli matched for low-level properties such as contrast and luminance? If not, I would suggest grouping them by e.g. luminance and comparing these groups with the same methods to test this potential confound.

The authors speak of a modulation of theta oscillations by emotional stimuli, yet the dominant structures in the time-frequency spectra of the habenula look more like an evoked response (wide band, low latency, inverted "V" shape). It might be worthwhile inspecting the stimulus-locked temporal average of the habenula signal and the time-frequency spectrum of that average to better understand the nature of this response.

In the discussion, it is mentioned that clinical scores and electrophysiological effects did not correlate. This important information should appear in the abstract and results, too.

It seems that the authors have used Medtronic extensions for externalization. Unless the custom-made non-magnetic version is used (please specify), these cause tremendous artifacts in the MEG signal. In this case, the authors should provide information on the degree of artifact contamination and the effectiveness of data cleaning.

In Figure 2A it seems that some electrodes were not placed in the habenula. Were these excluded? If not, why not?

Regarding my last comment in the public review about the somewhat unexpected timing of the coherence increase: Could this timing be an "artifact" of averaging over all MEG channels? Maybe this all-channel average mixes frontal and occipital coherence changes which themselves occur at more intuitive time points?

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

Thank you for resubmitting your work entitled "Increased theta/alpha synchrony in the habenula-prefrontal network with negative emotional stimuli in human patients" for further consideration by eLife. Your revised article has been evaluated by Floris de Lange (Senior Editor) and a Reviewing Editor.

Reviewer #1:

The authors have adequately addressed most of my concerns. I do however remain concerned about whether the early transient activity in the theta/alpha band can be considered oscillatory given how short the effect is. I do appreciate that the authors now demonstrate that the activity is not only phase-locked to the stimulus, but that does not make it oscillatory per se. I suggest that the authors refer to this effect as transient theta/alpha band activity in the Abstract and though out the manuscript. I do also understand the coherence effect is later and more sustained and more like to reflect neuronal oscillations.

I realized it is not made clear how the correlations are calculated (Pearson I assume?). If so, are there concerns on outliers? Please clarify.

Reviewer #2:

My thanks to the authors for providing detailed and good answers to my comments. I have just two things left to say:

Regarding my comment 2: I am aware that the reported correlation includes all conditions, but based on the red dots in Figure 6C, it does not seem that the relationship is particularly strong for the "negative" condition. This can be checked by performing separate correlations for all three conditions. More importantly, I would ask the authors to repeat their power-coherence correlation for other time windows to check whether the correlation is specific to the time windows of interest. If it turns out the correlation is neither condition- nor time-specific, I would think that it reflects general similarity between power and coherence as well as temporal autocorrelation in both measures.

Regarding my comment 6: With the frontal channel selection we now see a difference between conditions in coherence emerging much earlier in the trial. To me, this shows that it is better use exactly the same channels/locations when presenting power and coherence.

Reviewer #3:

It's not my intention to be difficult, but it doesn't seem that the authors have addressed my concerns, or done any additional analyses to add statistical confidence on results that are based on a very low sample (inevitably, given the nature of the data).

To give an example, Figure 6 is exactly as before (except that now the number of subjects is noted in the caption). The authors have not followed my suggestion of performing nonparametric permutation testing on the correlation reported in Figure 6C, and it remains an open question whether this is driven by the blue dot at the bottom of the plot (potentially an outlier). Unfortunately, that other coherence measures (eg occipital – habenula) were not significant is not a valid argument; i.e. that does not give any information about whether this particular relationship is truly significant.

Essential Revisions:

Please pay special attention to

1) Concerns on volume conduction

Please see detailed response to Reviewer 1 Public Review Recommendation 1. In summary, we have tried different ways to confirm the location of the electrode and the contacts selected for analysis. We also tried to use bipolar recordings from the neighboring macroelectrodes in order to minimize activities from other areas due to volume conduction. However, considering the size of the habenula, the size of the macroelectrode contacts, and the resolution of the MRI, we still cannot completely exclude the possibility that our recordings can be contaminated by activities from other areas due to volume conduction. Therefore, we have now added in the following in the ‘Limitation’ of the Discussion:

“Another caveat we would like to acknowledge is that the human habenula is a small region. Existing data from structural MRI scans reported combined habenula (the sum of the left and right hemispheres) volumes of ~ 30–36 mm³ (Savitz et al., 2011a; Savitz et al., 2011b meaning that each habenula has the size of 2~3 mm in each dimension, which may be even smaller than the standard functional MRI voxel size (Lawson et al., 2013). The size of the habenula is also small relative to the standard DBS electrodes (as shown in Figure 2A). The electrodes used in this study (Medtronic 3389) have electrode diameter of 1.27 mm with each contact length of 1.5 mm, and contact spacing of 0.5 mm. We have tried different ways to confirm the location of the electrode and to select the contacts that is within or closest to the habenula: 1.) the MRI was co-registered with a CT image (General Electric, Waukesha, WI, USA) with the Leksell stereotactic frame to obtain the coordinate values of the tip of the electrode; 2.) Post-operative CT was co-registered to pre-operative T1 MRI using a two-stage linear registration using Lead-DBS software. We used bipolar signals constructed from neighbouring macroelectrode recordings, which have been shown to detect locally generated potentials from subthalamic nucleus and especially when the macroelectrodes are inside the subthalamic nucleus (Marmor et al., 2017). Considering that not all contacts for bipolar LFP construction are in the habenula in this study, as shown in Figure 2, we cannot exclude the possibility that the activities we measured are contaminated by activities from neighbouring areas through volume conduction. In particular, the human habenula is surrounded by thalamus and adjacent to the posterior end of the medial dorsal thalamus, so we may have captured activities from the medial dorsal thalamus. However, we also showed that those bipolar LFPs from contacts in the habenula tend to have a peak in the theta/alpha band in the power spectra density (PSD); whereas recordings from contacts outside the habenula tend to have extra peak in beta frequency band in the PSD. This supports the habenula origin of the emotional valence related changes in the theta/alpha activities reported here.”

References:

Savitz JB, Bonne O, Nugent AC, Vythilingam M, Bogers W, Charney DS, et al. Habenula volume in post-traumatic stress disorder measured with high-resolution MRI. Biology of Mood & Anxiety Disorders 2011a; 1(1): 7.

Savitz JB, Nugent AC, Bogers W, Roiser JP, Bain EE, Neumeister A, et al. Habenula volume in bipolar disorder and major depressive disorder: a high-resolution magnetic resonance imaging study. Biological Psychiatry 2011b; 69(4): 336-43.

Lawson RP, Drevets WC, Roiser JP. Defining the habenula in human neuroimaging studies. NeuroImage 2013; 64: 722-7.

Marmor O, Valsky D, Joshua M, Bick AS, Arkadir D, Tamir I, et al. Local vs. volume conductance activity of field potentials in the human subthalamic nucleus. Journal of Neurophysiology 2017; 117(6): 2140-51.

2) The issue of evoked versus oscillatory responses

Please see detailed response to Reviewer 1 Public Review Recommendation 2. In summary, we have now quantified the evoked related potential (ERPs) for different emotional valence conditions separately for each subject, but we did not found emotional valence related modulation in the amplitude or the latency of the ERPs. In addition, we have also calculated induced oscillations only (after removing the evoked activity). And we found similar difference in the theta/alpha band activities caused by the emotional valence condition in the induced activities as we showed in the main results. So we argue that the difference in the early time window we reported in the manuscript is not dominated by evoked activity. The figures of the new analyses and following have now been added in the main text:

“In addition, we tested whether stimuli-related habenula LFP modulations primarily reflect a modulation of oscillations, which is not phase-locked to stimulus onset, or, alternatively, if they are attributed to evoked event-related potential (ERP). We quantified the ERP for each emotional valence condition for each habenula. There was no significant difference in ERP latency or amplitude caused by different emotional valence stimuli (Figure 3—figure supplement 1). In addition, when only considering the non phase-locked activity by removing the ERP from the time series before frequency-time decomposition, the emotional valence effect (presented in Figure 3—figure supplement 2) is very similar to those shown in Figure 3. These additional analyses demonstrated that the emotional valence effect in the LFP signal is more likely to be driven by non-phase-locked (induced only) activity.”

3) Concerns on low-level features of the visual stimuli driving the evoked effects

Please see detailed response to Reviewer 2 Recommendations for the authors 2. Unfortunately, low-level properties, such as contrast brightness and complexity of the figures are not measured or reported in the Chinese Affective Figure system. Only very dark or bright pictures were excluded from the paradigm. We have now added this information in the method (Line 367 – 369), and have also acknowledged this potential confounding in the discussion. However, we also expect that the effect of low-level properties in the figures, such as luminance or color, may have more prominent effect in the visual cortex measured from the parietal and occipital regions as shown in previous studies Yener et al., 2009; Haigh et al., 2018; Peter et al., 2019; Perry et al., 2020.

References:

Yener GG, Guntekin B, Tulay E, Basar E. A comparative analysis of sensory visual evoked oscillations with visual cognitive event related oscillations in Alzheimer's disease. Neuroscience Letters 2009; 462(3): 193-7.

Haigh SM, Cooper NR, Wilkins AJ. Chromaticity separation and the alpha response. Neuropsychologia 2018; 108: 1-5.

Peter A, Uran C, Klon-Lipok J, Roese R, van Stijn S, Barnes W, et al. Surface color and predictability determine contextual modulation of V1 firing and γ oscillations. eLife 2019; 8.

Perry G, Taylor NW, Bothwell PCH, Milbourn CC, Powell G, Singh KD. The γ response to colour hue in humans: Evidence from MEG. PloS One 2020; 15(12): e0243237.

4) Robustness of findings in the light of saccadic artefacts and statistical assessment

We have now addressed the concerns about the saccadic artefacts, please see the detailed response to the Reviewer 1 Public Review Recommendation 4.

We have now addressed the concerns about the robustness of the finding due to the average trial number per condition per participant, please see the response to the Reviewer 2 (Public Review Detailed Comment 1).

We have now addressed the concerns about the signal-to-noise ratio and ‘artefacts due to the extension cable for LFP recordings’, please see the response to Reviewer 2 (Recommendations for the authors, Point #4).

We have now addressed the concerns about the robustness of the statistics, please see the response to the Reviewer 3 Recommendation.

Reviewer #2 (Recommendations for the authors):

Were the stimuli matched for low-level properties such as contrast and luminance? If not, I would suggest grouping them by e.g. luminance and comparing these groups with the same methods to test this potential confound.

We thank the reviewer for the comment. Unfortunately, low-level properties, such as contrast brightness and complexity of the figures are not measured or reported in the Chinese Affective Figure system. Only very dark or bright pictures were excluded from the paradigm. We have now added this information in the method: “As low-level properties of the figures, such as contrast brightness and complexity of the figures are not measured or reported in the CAPS, only very dark or bright pictures were excluded from the paradigm.” We have also acknowledged this potential confounding in the discussion. However, we also expect that the effect of low-level properties in the figures, such as luminance or color, may have more prominent effect in the visual cortex measured from the parietal and occipital regions as shown in previous studies Yener et al., 2009; Haigh et al., 2018; Peter et al., 2019; Perry et al., 2020.

References:

Yener GG, Guntekin B, Tulay E, Basar E. A comparative analysis of sensory visual evoked oscillations with visual cognitive event related oscillations in Alzheimer's disease. Neuroscience Letters 2009; 462(3): 193-7.

Haigh SM, Cooper NR, Wilkins AJ. Chromaticity separation and the alpha response. Neuropsychologia 2018; 108: 1-5.

Peter A, Uran C, Klon-Lipok J, Roese R, van Stijn S, Barnes W, et al. Surface color and predictability determine contextual modulation of V1 firing and γ oscillations. eLife 2019; 8.

Perry G, Taylor NW, Bothwell PCH, Milbourn CC, Powell G, Singh KD. The γ response to colour hue in humans: Evidence from MEG. PloS One 2020; 15(12): e0243237.

The authors speak of a modulation of theta oscillations by emotional stimuli, yet the dominant structures in the time-frequency spectra of the habenula look more like an evoked response (wide band, low latency, inverted "V" shape). It might be worthwhile inspecting the stimulus-locked temporal average of the habenula signal and the time-frequency spectrum of that average to better understand the nature of this response.

We agree with the reviewer that the activity increase in the first time window with short latency after the stimuli onset is very transient and not band-limited. This raise the question that whether this is oscillatory or a transient evoked activity. Please see the response to Reviewer 1 Recommendation 2 for the details of further analysis. In summary, we have now quantified the ERP but found no effect of emotional valence on the amplitude or latency of the ERP. We have also looked at the non-phase locked activities removing the influence of the ERP from the total power spectra, we didn’t observe any effect of the stimulus emotional valence in the pure induced (non-phase-locked) power spectra. Therefore, we argue that the theta/alpha increase with negative emotional stimuli we observed in both habenula and prefrontal cortex 0-500 ms after stimuli onset are not dominated by visual or other ERP.

In the discussion, it is mentioned that clinical scores and electrophysiological effects did not correlate. This important information should appear in the abstract and results, too.

Following the reviewer’s suggestion, we have now added the sentence in the results: “Furthermore, we also investigated the relationship between the neural characteristics we observed and the clinical symptoms. However, none of the electrophysiological effects we observed correlated with clinical scores of depression (the Beck Depression Inventory score or Hamilton Depression Rating Scale score) measured before the surgery across patients after correcting for multiple correction.” However due to the strict word limit in the abstract, we were not able to include this in the abstract.

It seems that the authors have used Medtronic extensions for externalization. Unless the custom-made non-magnetic version is used (please specify), these cause tremendous artifacts in the MEG signal. In this case, the authors should provide information on the degree of artifact contamination and the effectiveness of data cleaning.

Medtronic extensions for externalization was used in this study, NOT the custom-made non-magnetic version. The extension wires can reduce signal-to-noise ratio in the MEG recording. Therefore the spatiotemporal Signal Space Separation (tSSS) method (Taulu and Simola, 2006) implemented in the MaxFilter software (Elekta Oy, Helsinki, Finland) has been applied in this study to suppress strong magnetic artifacts caused by extension wires. This method has been proved to work well in de-noising the magnetic artifacts and movement artifacts in MEG data in our previous studies (Cao et al., 2019; Cao et al., 2020). In addition, the beamforming method proposed by several studies (Litvak et al., 2010; Hirschmann et al., 2011; Litvak et al., 2011) has been used in this study. In Litvak et al., 2010, the artifacts caused by DBS extension wires was described in detail and the beamforming was demonstrated to effectively suppress artifacts and thereby enable both localization of cortical sources coherent with the deep brain nucleus. We have now added more details and references about the data cleaning and the beamforming method in the main text. With the beamforming method, we did observe the standard movement-related modulation in the β frequency band in the motor cortex, shown in Figure 5 - figure supplement 1. This suggests that the beamforming method did work well to suppress the artefacts. The figure on movement-related modulation in the motor cortex in the MEG signals have now been added as a supplementary figure to demonstrate the effect of the beamforming.

References:

Taulu S, Simola J. Spatiotemporal signal space separation method for rejecting nearby interference in MEG measurements. Physics in Medicine and Biology 2006; 51(7): 1759-68.

Cao C, Huang P, Wang T, Zhan S, Liu W, Pan Y, et al. Cortico-subthalamic Coherence in a Patient With Dystonia Induced by Chorea-Acanthocytosis: A Case Report. Frontiers in Human Neuroscience 2019; 13: 163.

Cao C, Li D, Zhan S, Zhang C, Sun B, Litvak V. L-dopa treatment increases oscillatory power in the motor cortex of Parkinson's disease patients. NeuroImage Clinical 2020; 26: 102255.

Litvak V, Eusebio A, Jha A, Oostenveld R, Barnes GR, Penny WD, et al. Optimized beamforming for simultaneous MEG and intracranial local field potential recordings in deep brain stimulation patients. NeuroImage 2010; 50(4): 1578-88.

Litvak V, Jha A, Eusebio A, Oostenveld R, Foltynie T, Limousin P, et al. Resting oscillatory cortico-subthalamic connectivity in patients with Parkinson's disease. Brain 2011; 134(Pt 2): 359-74.

Hirschmann J, Ozkurt TE, Butz M, Homburger M, Elben S, Hartmann CJ, et al. Distinct oscillatory STN-cortical loops revealed by simultaneous MEG and local field potential recordings in patients with Parkinson's disease. NeuroImage 2011; 55(3): 1159-68.

In Figure 2A it seems that some electrodes were not placed in the habenula. Were these excluded? If not, why not?

We didn’t exclude any electrode shown in Figure 2A because we showed that those electrodes from which only alpha/theta peaks were detected are well placed in the habenula area based on the Lead-DBS software. However, we noticed that those electrodes from which large beta band peak were present in the power spectra density and seem to be outside the habenula, we can still detect another peak of smaller amplitude in the theta/alpha band. The human habenula is a small region. Existing data from structural MRI scans reported combined habenula (the sum of the left and right hemisphere) volumes of ~ 30–36 mm³ (Savitz et al., 2011a; Savitz et al., 2011b) which means each habenula has the size of 2~3 mm in each dimension, which may be even smaller than the standard functional MRI voxel size (Lawson et al. 2013). The Lead-DBS software gives a good indication on the location of the electrode relative to a standardised atlas of habenula: THOMAS Atlas (Su et al., 2019). But we feel that it might not be valid to use Lead-DBS results to provide the binary selection criteria for such a small structure, as the electrode outside the habenula as shown in Figure 2 also have contacts which are at most 1 or 2 pixels away from the habenula. We have now also acknowledged that we cannot exclude the possibility that the activities we measured are contaminated by activities from neighbouring areas through volume conduction. In particular, the human habenula is surrounded by thalamus and adjacent to the posterior end of the medial dorsal thalamus, so the β band activities captured by some electrodes may come from the medial dorsal thalamus. We have now added this as part of the ‘Limitation’. Please see detailed response to Reviewer 1 Comment 1.

References:

Savitz JB, Bonne O, Nugent AC, Vythilingam M, Bogers W, Charney DS, et al. Habenula volume in post-traumatic stress disorder measured with high-resolution MRI. Biology of Mood & Anxiety Disorders 2011a; 1(1): 7.

Savitz JB, Nugent AC, Bogers W, Roiser JP, Bain EE, Neumeister A, et al. Habenula volume in bipolar disorder and major depressive disorder: a high-resolution magnetic resonance imaging study. Biological Psychiatry 2011b; 69(4): 336-43.

Lawson RP, Drevets WC, Roiser JP. Defining the habenula in human neuroimaging studies. NeuroImage 2013; 64: 722-7.

Su JH, Thomas FT, Kasoff WS, Tourdias T, Choi EY, Rutt BK, et al. Thalamus Optimized Multi Atlas Segmentation (THOMAS): fast, fully automated segmentation of thalamic nuclei from structural MRI. NeuroImage 2019; 194: 272-82.

Regarding my last comment in the public review about the somewhat unexpected timing of the coherence increase: Could this timing be an "artifact" of averaging over all MEG channels? Maybe this all-channel average mixes frontal and occipital coherence changes which themselves occur at more intuitive time points?

We have now looked at the cross-trial coherence changes for frontal cortex-habenula and occipital cortex-habenula separately. When only MEG channels from the frontal cortex was considered, the coherence with habenula showed significant increase around 1 s which is very similar (Figure S7 related to Figure 6A) as we reported in the main manuscript. But the coherence first peaked at around 500 ms after stimuli onset, which is still after the initial activities increase we observed in both the prefrontal cortex and the habenula. In contrast, we didn’t observe any changes in the coherence between occipital channels and habenula. So the timing with significant difference in the prefrontal cortex-habenula area coherence is not an ‘artifact’ of averaging all MEG channels. The following has now been added in the main text: “We also performed the same analysis for cross-trial cortical-habenula coherence averaged across prefrontal channels and occipital channels separately. The emotional valence effect on the coherence was only observed in the frontal channels not in the occipital channels, as shown in Figure 6—figure supplement 1.”

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

Reviewer #1:

The authors have adequately addressed most of my concerns. I do however remain concerned about whether the early transient activity in the theta/alpha band can be considered oscillatory given how short the effect is. I do appreciate that the authors now demonstrate that the activity is not only phase-locked to the stimulus, but that does not make it oscillatory per se. I suggest that the authors refer to this effect as transient theta/alpha band activity in the Abstract and though out the manuscript. I do also understand the coherence effect is later and more sustained and more like to reflect neuronal oscillations.

We thank the reviewer for the comment. We have now made the revisions accordingly. The term “oscillatory activity” induced by stimulus has been changed to “transient activity” consistently through the abstract and main text, especially when referring to the modulation from 100 to 500 ms after the stimuli onset. In addition, we have now added the following text in the Results: “It should be noted that there was an increase in a broad band activity at short latency (from 100 to 500 ms) after the stimuli onset (Figure 3A-C). This raises the question as to whether the emotional valence related modulation observed in Figure 3D, especially the cluster at short latency, reflects a modulation of oscillations, which is not phase-locked to stimulus onset, or, alternatively, are they attributable to an evoked event-related potential (ERP).’. The following sentence has also been included: ‘These additional analyses demonstrated that the emotional valence effect in the LFP signal is more likely to be driven by non-phase-locked (induced only) activity, even though the possibility of the contribution from transient evoked potentials still cannot be completed excluded. Therefore, we refer to the activities in the habenula LFPs that are modulated by emotional valence at short latency after stimulus onset as ‘activity’ rather than ‘oscillations’.”

I realized it is not made clear how the correlations are calculated (Pearson I assume?). If so, are there concerns on outliers? Please clarify.

In the Figure 6C, the correlation was estimated by linear mixed-effect modelling method. Regarding the concerns on outliers, the nonparametric permutation test that is robust to outliers performed on the Pearson’s correlation using the PERMUTOOLS toolbox in MATLAB (https://github.com/mickcrosse/PERMUTOOLS) has now been used to evaluate correlation between the theta power and coherence. The new test has showed significant correlation between PFC-habenula theta coherence and habenula theta power as well (R = 0.3224 (95% confidence interval: [0.0422 0.5557]), p = 0.03). In addition, to address the comments from other reviewers, we have now applied the permutation based correlation test on data from different emotional valence conditions separately, and also applied the correlation test to the PFC-habenula coherence at other time windows. The following has now been added in the main text:

Methods: “The nonparametric permutation test that is robust to outliers was performed on the Pearson’s correlations using the PERMUTOOLS toolbox in MATLAB (https://github.com/mickcrosse/PERMUTOOLS) to further evaluate correlations when considering all data from different subjects together. The estimated R value with 95% confidence interval, and estimated p values of the test were reported.”

Results: “The nonparametric permutation test that is robust to outliers (https://github.com/mickcrosse/PERMUTOOLS) was also used to evaluate the correlation between the theta power and coherence when data from all participants and all emotional conditions were considered together. This confirmed significant correlation as well (R = 0.3224 (95% confidence interval: [0.0422 0.5557]), p = 0.03). On the other hand, when data from different emotional conditions were considered separately, none of the separate correlations between theta coherence at 800–1300 ms and habenula theta power at 2700–3300 ms were significant: (R = 0.3405 (95% confidential interval: [-0.1867 0.7154]), p = 0.2020 for Neutral; R = 0.3846 (95% confidential interval: [-0.1373 0.7394]), p = 0.1474 for Positive; R = -0.1655 (95% confidential interval: [-0.6111 0.3597]), p = 0.5474 for Negative). In addition, we tested whether this coherence-power correlation was specific to the time window identified based on Figure 6A. To do so, we quantified the correlation between the habenula theta power at 2700–3300 ms and the habenula-PFC theta coherence at -200–300 ms, 300–800 ms, 1300–1800 ms, and 1800–2300ms separately. None of the habenula-PFC coherences at other time windows correlated with habenula theta at 2700–3300 ms. We acknowledge that the effect shown in Figure 6C is weak and would not survive correction for multiple comparison. However, the selection of time window for the test shown in Figure 6C was based on the previous test shown in Figure 6A, not based on multiple tests.”

Reviewer #2:

My thanks to the authors for providing detailed and good answers to my comments. I have just two things left to say:

Regarding my comment 2: I am aware that the reported correlation includes all conditions, but based on the red dots in Figure 6C, it does not seem that the relationship is particularly strong for the "negative" condition. This can be checked by performing separate correlations for all three conditions. More importantly, I would ask the authors to repeat their power-coherence correlation for other time windows to check whether the correlation is specific to the time windows of interest. If it turns out the correlation is neither condition- nor time-specific, I would think that it reflects general similarity between power and coherence as well as temporal autocorrelation in both measures.

We thank the reviewer for the comments. We have now applied the permutation based correlation test on data from different emotional valence conditions separately, and also applied the correlation test to the habenula theta power and PFC-habenula coherence at other time windows. When the PFC-habenula theta band coherence at 800–1300 ms was considered, none of the separate correlations between PFC-habenula coherence and habenula theta power at 2700–3300 ms for different emotional valence conditions were significant (R² = 0.1156, p = 0.1969 for neutral valence condition; R² = 0.1475, p = 0.1413 for positive valence condition; R² = 0.0274, p = 0.5402 for negative valence condition). Results were similar when nonparametric tests (https://github.com/mickcrosse/PERMUTOOLS) based on permutation were used: the test confirmed significant correlation when data from all conditions were considered together (R = 0.3224 (95% confidence interval: [0.0422 0.5557]), p = 0.03), but this was not the case for correlation when data from different emotional valence conditions were considered separately: R = 0.3405 (95% confidential interval: [-0.1867 0.7154]), p = 0.2020 for neutral; R = 0.3846 (95% confidential interval: [-0.1373 0.7394]), p = 0.1474 for positive; R = -0.1655 (95% confidential interval: [-0.6111 0.3597]), p = 0.5474 for negative condition.

We also investigated whether the PFC-habenula coherence at other time windows, such as -200–300 ms, 300–800 ms, 1300–1800 ms, and 1800–2300ms is correlated with habenula theta power at the later time window (around 3000 ms). The Linear mixed-effect modelling analysis showed that the power-coherence correlations for other time windows are not significant (Figure S1), indicating the correlation between power and coherence is time-specific.

We have now added these results to the main text:

Methods: “The nonparametric permutation test that is robust to outliers was performed on the Pearson’s correlations using the PERMUTOOLS toolbox in MATLAB (https://github.com/mickcrosse/PERMUTOOLS) to further evaluate correlations when considering all data from different subjects together. The estimated R value with 95% confidence interval, and estimated p values of the test were reported.”

Results: “The nonparametric permutation test that is robust to outliers (https://github.com/mickcrosse/PERMUTOOLS) was also used to evaluate the correlation between the theta power and coherence when data from all participants and all emotional conditions were considered together. This confirmed significant correlation as well (R = 0.3224 (95% confidence interval: [0.0422 0.5557]), p = 0.03). On the other hand, when data from different emotional conditions were considered separately, none of the separate correlations between theta coherence at 800–1300 ms and habenula theta power at 2700–3300 ms were significant: (R = 0.3405 (95% confidential interval: [-0.1867 0.7154]), p = 0.2020 for Neutral; R = 0.3846 (95% confidential interval: [-0.1373 0.7394]), p = 0.1474 for Positive; R = -0.1655 (95% confidential interval: [-0.6111 0.3597]), p = 0.5474 for Negative). In addition, we tested whether this coherence-power correlation was specific to the time window identified based on Figure 6A. To do so, we quantified the correlation between the habenula theta power at 2700–3300 ms and the habenula-PFC theta coherence at -200–300 ms, 300–800 ms, 1300–1800 ms, and 1800–2300ms separately. None of the habenula-PFC coherences at other time windows correlated with habenula theta at 2700–3300 ms. We acknowledge that the effect shown in Figure 6C is weak and would not survive correction for multiple comparison. However, the selection of time window for the test shown in Figure 6C was based on the previous test shown in Figure 6A, not based on multiple tests.”

Regarding my comment 6: With the frontal channel selection we now see a difference between conditions in coherence emerging much earlier in the trial. To me, this shows that it is better use exactly the same channels/locations when presenting power and coherence.

The earlier time window (from 500 to 600 ms) with significant difference between frontal MEG-habenula LFP coherence for conditions was very brief. In addition, this time window was not significant when the other coupling method (imaginary coherence) was used to evaluate the connectivity between frontal MEG and habenula (Figure S2). The time window in which the PFC-habenula connectivity was consistently modulated by emotional valence was only around 1000 ms (8000-1200 ms), but not at 500-600 ms. Therefore, we still focused on this time window for the other analysis shown in Figure 6.

Reviewer #3:

It's not my intention to be difficult, but it doesn't seem that the authors have addressed my concerns, or done any additional analyses to add statistical confidence on results that are based on a very low sample (inevitably, given the nature of the data).

We are sorry for not addressing the question on statistical confidence in the previous version. We appreciate the concern of the reviewer given the small sample size we have. We have now added 95% confidence intervals into the main text in the ‘Results’. For example, we have now used a non-parametric bootstrap test for the post-hoc analysis shown in Figure 3 G and H, and have now shown the confidence interval of the test. We have now made the following changes in the main text: “The average powers in the determined frequency band and time window identified by the cluster-based permutation method between different valence conditions were further compared using post-hoc paired-sample permutation t-tests with the function in the PERMUTOOLS toolbox in MATLAB (https://github.com/mickcrosse/PERMUTOOLS). The t value and associated confidence interval of the difference, as well as the p value of the paired-sample permutation t-test are reported.”. The statistical results were updated in Figure 3.

In addition, we have now included the confidence interval in the effects of the multi-level modelling. For example, Page 9 Line 176: ‘k = 0.2434 ± 0.1031 (95% confidence interval [0.0358 0.4509])’ and the Table 2.

We have also added 95% confidential intervals for the results shown in Figure 7. The 95% CI for the partial correlation coefficient between valence rating and frontal theta/alpha power, between valence rating and coherence, between valence rating and habenula theta power, controlling for arousal rating are [-0.6946, -0.1491], [-0.7738, -0.2429], [-0.7519, -0.2353], respectively. The 95% confidence intervals for the partial correlation coefficient between arousal rating and frontal theta/alpha power, between arousal rating and coherence, between arousal rating and habenula theta power, controlling for valence rating are [-0.5517, 0.4273], [-0.3253, 0.4634], [-0.3748, 0.3926], respectively.

To give an example, Figure 6 is exactly as before (except that now the number of subjects is noted in the caption). The authors have not followed my suggestion of performing nonparametric permutation testing on the correlation reported in Figure 6C, and it remains an open question whether this is driven by the blue dot at the bottom of the plot (potentially an outlier). Unfortunately, that other coherence measures (eg occipital – habenula) were not significant is not a valid argument; i.e. that does not give any information about whether this particular relationship is truly significant.

Following the reviewer’s suggestion, the nonparametric permutation test that is robust to outliers has now been used based on the Pearson’s correlation between the theta power and coherence using the PERMUTOOLS toolbox in MATLAB (https://github.com/mickcrosse/PERMUTOOLS), and confirms significant coherence (R = 0.3224 (95% confidence interval: [0.0422 0.5557]), p = 0.03). In addition, the linear mixed-effect modelling analysis after removing the dot at the bottom of the plot showed significant correlation between the coherence and the theta power (R² = 0.0832, p = 0.0447).

We have now added the following to the main text:

Methods: ‘The nonparametric permutation test that is robust to outliers was performed on the Pearson’s correlations using the PERMUTOOLS toolbox in MATLAB (https://github.com/mickcrosse/PERMUTOOLS) to further evaluate correlations when considering all data from different subjects together. The estimated R value with 95% confidence interval, and estimated p values of the test were reported.’

Results: ‘The nonparametric permutation test that is robust to outliers (https://github.com/mickcrosse/PERMUTOOLS) was also used to evaluate the correlation between the theta power and coherence when data from all participants and all emotional conditions were considered together. This confirmed significant correlation as well (R = 0.3224 (95% confidence interval: [0.0422 0.5557]), p = 0.03). On the other hand, when data from different emotional conditions were considered separately, none of the separate correlations between theta coherence at 800–1300 ms and habenula theta power at 2700–3300 ms were significant: (R = 0.3405 (95% confidential interval: [-0.1867 0.7154]), p = 0.2020 for Neutral; R = 0.3846 (95% confidential interval: [-0.1373 0.7394]), p = 0.1474 for Positive; R = -0.1655 (95% confidential interval: [-0.6111 0.3597]), p = 0.5474 for Negative). In addition, we tested whether this coherence-power correlation was specific to the time window identified based on Figure 6A. To do so, we quantified the correlation between the habenula theta power at 2700–3300 ms and the habenula-PFC theta coherence at -200–300 ms, 300–800 ms, 1300–1800 ms, and 1800–2300ms separately. None of the habenula-PFC coherences at other time windows correlated with habenula theta at 2700–3300 ms. We acknowledge that the effect shown in Figure 6C is weak and would not survive correction for multiple comparison. However, the selection of time window for the test shown in Figure 6C was based on the previous test shown in Figure 6A, not based on multiple tests.’

Author details

1. Yongzhi Huang

   1. Academy of Medical Engineering and Translational Medicine, Tianjin University, Tianjin, China
   2. Nuffield Department of Surgical Sciences, University of Oxford, Oxford, United Kingdom

   Contribution

   Conceptualization, Data curation, Formal analysis, Methodology, Writing - original draft, Writing - review and editing

   Contributed equally with

   Bomin Sun

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-2503-1589

2. Bomin Sun

   Department of Neurosurgery, Affiliated Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China

   Contribution

   Resources, Funding acquisition, Investigation, Project administration, Writing - review and editing

   Contributed equally with

   Yongzhi Huang

   Competing interests

   No competing interests declared

3. Jean Debarros

   Medical Research Council (MRC) Brain Network Dynamics Unit at the University of Oxford, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, United Kingdom

   Contribution

   Formal analysis, Writing - review and editing

   Competing interests

   No competing interests declared

4. Chao Zhang

   Department of Neurosurgery, Affiliated Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China

   Contribution

   Investigation, Writing - review and editing

   Competing interests

   No competing interests declared

5. Shikun Zhan

   Department of Neurosurgery, Affiliated Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China

   Contribution

   Investigation, Writing - review and editing

   Competing interests

   No competing interests declared

6. Dianyou Li

   Department of Neurosurgery, Affiliated Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China

   Contribution

   Investigation, Writing - review and editing

   Competing interests

   No competing interests declared

7. Chencheng Zhang

   Department of Neurosurgery, Affiliated Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China

   Contribution

   Investigation, Writing - review and editing

   Competing interests

   No competing interests declared

8. Tao Wang

   Department of Neurosurgery, Affiliated Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China

   Contribution

   Investigation, Writing - review and editing

   Competing interests

   No competing interests declared

9. Peng Huang

   Department of Neurosurgery, Affiliated Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China

   Contribution

   Investigation, Writing - review and editing

   Competing interests

   No competing interests declared

10. Yijie Lai

    Department of Neurosurgery, Affiliated Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China

    Contribution

    Formal analysis, Writing - review and editing

    Competing interests

    No competing interests declared

11. Peter Brown

    Medical Research Council (MRC) Brain Network Dynamics Unit at the University of Oxford, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, United Kingdom

    Contribution

    Formal analysis, Writing - review and editing

    Competing interests

    No competing interests declared

    "This ORCID iD identifies the author of this article:" 0000-0002-5201-3044

12. Chunyan Cao

    Department of Neurosurgery, Affiliated Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China

    Contribution

    Conceptualization, Data curation, Formal analysis, Supervision, Funding acquisition, Investigation, Project administration, Writing - review and editing

    For correspondence

    ccy40646@rjh.com.cn

    Competing interests

    No competing interests declared

13. Huiling Tan

    Medical Research Council (MRC) Brain Network Dynamics Unit at the University of Oxford, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, United Kingdom

    Contribution

    Conceptualization, Data curation, Formal analysis, Supervision, Funding acquisition, Investigation, Writing - original draft, Project administration, Writing - review and editing

    For correspondence

    huiling.tan@ndcn.ox.ac.uk

    Competing interests

    No competing interests declared

    "This ORCID iD identifies the author of this article:" 0000-0001-8038-3029

• 747

  Page views

• 131

  Downloads

• 6

  Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.