
==== Front
Int J Clin Health Psychol
Int J Clin Health Psychol
International Journal of Clinical and Health Psychology : IJCHP
1697-2600
2174-0852
Asociacion Espanola de Psicologia Conductual

S1697-2600(24)00060-7
10.1016/j.ijchp.2024.100495
100495
Original Article
Neurophysiological markers of disease severity and cognitive dysfunction in major depressive disorder: A TMS-EEG study
Li Deyang ab
Li Xingxing abcd
Li Jiaxin ab
Liu Junyao ab
Luo Ruichenxi ab
Li Yanli ab
Wang Dongmei wangdm@psych.ac.cn
ab⁎
Zhou Dongsheng wyzhouds@sina.com
cd⁎⁎
Zhang Xiang-Yang zhangxy@psych.ac.cn
ab⁎
a CAS Key Laboratory of Mental Health, Institute of Psychology, Chinese Academy of Sciences, Beijing, China
b Department of Psychology, University of Chinese Academy of Sciences, Beijing, China
c Department of psychiatry, Affiliated Kangning Hospital of Ningbo University, Ningbo, Zhejiang, China
d Department of psychiatry, Ningbo Kangning Hospital, Ningbo, Zhejiang, China
⁎ Corresponding authors at: 16 Lincui Road, Chaoyang District, Beijing 100101, China. wangdm@psych.ac.cnzhangxy@psych.ac.cn
⁎⁎ Corresponding author at: Department of psychiatry, Affiliated Kangning Hospital of Ningbo University, Ningbo, Zhejiang 315201, China. wyzhouds@sina.com
30 8 2024
Jul-Sep 2024
30 8 2024
24 3 10049522 3 2024
12 8 2024
© 2024 The Authors. Published by Elsevier B.V.
2024

https://creativecommons.org/licenses/by-nc-nd/4.0/ This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
Background

Transcranial magnetic stimulation-electroencephalography (TMS-EEG) is a powerful technique to study the neuropathology and biomarkers of major depressive disorder (MDD). This study investigated cortical activity and its relationship with clinical symptoms and cognitive dysfunction in MDD patients by indexing TMS-EEG biomarkers in the dorsolateral prefrontal cortex (DLPFC).

Methods

133 patients with MDD and 76 healthy individuals participated in this study. Single-pulse TMS was performed on the left DLPFC to obtain TMS-evoked potential (TEP) indices. TMS-EEG waveforms and components were determined by global mean field amplitude. We used the Repeatable Battery for the Assessment of Neuropsychological Status (RBANS) to measure participants’ cognitive function.

Results

Patients with MDD had a lower excitatory P180 index compared to healthy controls, and P180 amplitude was negatively correlated with the severity of depressive and anxiety symptoms in patients with MDD. In the MDD group, P30 amplitude was negatively associated with RBANS Visuospatial/ Constructional index and total score.

Conclusions

TMS-EEG findings suggest that abnormal cortical excitation and inhibition induced by TMS on the DLPFC are associated with the severity of clinical symptoms and cognitive dysfunction in patients with MDD. P180 and P30 have the potential to serve as neurophysiological biomarkers of clinical symptoms and cognitive dysfunction in MDD patients, respectively.

Keywords

TMS-EEG
Major depressive disorder
Dorsolateral prefrontal cortex
Cognitive function
Biomarker
==== Body
pmcIntroduction

Major depressive disorder (MDD) is a prevalent mental illness defined by persistent affective, behavioral, and cognitive dysfunction (Subhas et al., 2023). The World Health Organization estimates that approximately 3.8 % of the global population suffers from depression (World Health Organization, 2023). MDD involves complex neuropathological mechanisms, including changes in neuronal structure, neurotransmission, and brain network connectivity. Notably, it has been demonstrated that the dorsolateral prefrontal cortex (DLPFC) is essential to the pathophysiology of MDD. Hypofunction of the left DLPFC in individuals with MDD has been found in different study modalities, such as neurotransmitter (Tran et al., 2023), functional imaging (Shen et al., 2015), and electroencephalography (EEG) (Kamishikiryo et al., 2022) studies. A recent meta-analysis showed that the left DLPFC is a node where placebo effect mechanisms and neuromodulatory anti-depressant mechanisms overlap, indicating that changes in the left DLPFC play a vital role in the treatment of MDD (Burke et al., 2022). In addition, MDD patients often present with anxiety symptoms, which may have a profound impact on the condition. It is also suggested that the DLPFC is involved in anxiety regulation through working memory processes (White et al., 2023). How the DLPFC contributes to the neuropathology of MDD patients warrants further exploration.

Recently, cognitive dysfunction has been recognized as one of the core features of MDD patients, involving abnormalities in the DLPFC, genomic variation, brain-derived neurotrophic factor (BDNF), and inflammatory mediators (Fossati, 2018; Morozova et al., 2022). The main domains of cognitive deficits in MDD patients include executive function, attention, and memory (Kriesche et al., 2023). Specifically, the DLPFC is considered part of the central executive network and plays a vital role in cognitive dysfunction. An imaging study found reduced functional connectivity in the left DLPFC and enhanced functional connectivity in the right DLPFC, suggesting that an imbalance between the left and right DLPFC is linked to cognitive dysfunction in MDD patients (Zhang et al., 2022). Another mechanism by which the DLPFC is involved in cognitive deficits in MDD patients is an imbalance between excitatory and inhibitory neurotransmission. Numerous studies have shown that the imbalance between excitation and inhibition in the DLPFC cerebral cortex is related to the severity of depressive symptoms (Biermann et al., 2022; Dhami et al., 2023). Notably, higher levels of activity in the left DLPFC in MDD patients may be a compensatory mechanism for the inactivation of the default mode network (DMN), which has been linked to impaired working memory and emotion regulation (Chen et al., 2023). However, further study is needed to elucidate the implications of excitatory and inhibitory modes of the DLPFC on symptoms and cognitive functioning in patients with MDD.

Transcranial magnetic stimulation (TMS) is a non-invasive brain stimulation therapy effective for MDD patients. TMS combined with EEG (TMS-EEG), enables the assessment of cortical evoked activation in the DLPFC. TMS-EEG is a powerful tool for investigating the causal role of cortical activation in mental processes. TMS evoked potentials (TEPs), including P30, N45, P60, N100, and P180, offer a window for exploring the localized neural assemblies of the DLPFC's electrical and organizational properties (Kallioniemi & Daskalakis, 2022). These components are associated with the excitatory glutamate system and the inhibitory γ-aminobutyric acid (GABA) system, reflecting the inhibitory-excitatory balance in neural circuits (Kallioniemi et al., 2022). For example, pharmacological studies have shown that P30, N45, and P180 amplitudes are associated with voltage-gated sodium channels (VGSCs) (Darmani et al., 2019; Kallioniemi & Daskalakis, 2022). P30 is thought to be mediated by GABAA receptors rather than glutamatergic N-methyl-d-aspartate (NMDA) receptor-mediated neurotransmission (Ferreri et al., 2011; Rogasch et al., 2020). P60 is thought to be linked to glutamatergic activity (Belardinelli et al., 2021). N45 and N100 are associated with inhibition involving GABAergic neurotransmission, reflecting GABAA and GABAB receptor activation, respectively (Premoli et al., 2014a; Rogasch et al., 2015). In addition, Global Mean Field Amplitude (GMFA) has been used to acquire the TEP component of the whole scalp to reduce the bias caused by limiting electrodes (Poorganji et al., 2023).

TMS-EEG has attracted increasing attention as a well-validated technique to study neuroplasticity in the DLPFC of MDD patients (Farzan, 2023). A cross-sectional study using GMFA showed that N45 amplitude holds tremendous potential as a neurophysiological biomarker of the DLPFC in detecting depression (Voineskos et al., 2019). The N100 component of the right DLPFC was also significantly predictive of MDD in adolescents (Dhami et al., 2020). Longitudinal studies have also utilized TMS-EEG biomarkers to predict the efficacy of various physical therapies for depression and to explain their neuroplasticity mechanisms, such as magnetic seizure therapy (Hadas et al., 2020) and intermittent theta pulse stimulation (Strafella et al., 2023). In addition, TMS-EEG has been used to explore neural biomarkers of cognitive function associated with the prefrontal cortex. Previous studies have found that TMS-evoked local cortical responses following stimulation of the left PFC are positively correlated with working memory and reasoning ability (Redondo-Camos et al., 2022). TMS-EEG studies have demonstrated abnormal excitability and functioning of the DLPFC in Alzheimer's disease (AD) patients, which have been associated with overall cognitive and executive functioning (Casarotto et al., 2011; Di Lazzaro et al., 2021; Joseph et al., 2021). However, there is still a lack of TMS-EEG biomarkers to predict cognitive deficits in MDD patients (Ferrarelli & Phillips, 2021). Identifying TMS-EEG biomarkers that predict symptom severity and cognitive deficits in MDD patients is vital for clinical diagnosis and treatment.

Given the critical role of DLPFC abnormalities in clinical symptoms and cognitive impairment of MDD patients, it is valuable to explore their neurophysiological mechanisms and biomarkers using TMS-EEG. To date, there is limited literature that simultaneously explores the relationship between TMS-EEG markers and symptom severity and cognitive impairment of patients with MDD. The present study aimed to (1) investigate the differences in DLPFC activity between MDD patients and healthy controls as measured by TMS-EEG, and (2) identify neurophysiological markers predicting clinical symptoms and cognitive function in MDD patients. We hypothesized that (1) TMS-EEG indexes would reflect hypofrontality in patients with MDD, and (2) TMS-EEG biomarkers would relate with the severity of clinical symptoms and cognitive deficits in MDD patients.

Methods

Participants

Between January 2022 and June 2023, 133 patients with MDD were recruited at Ningbo Kangning Hospital. All patients were right-handed. Inclusion criteria were (a) age 16 to 65 years; (b) fulfillment of DSM-5 diagnostic criteria for MDD; and (c) score ≥ 20 points on the Hamilton Depression Rating Scale-24 (HDRS-24). Exclusion criteria included (a) comorbidity with other psychiatric illnesses or neurological impairments; (b) the presence of serious medical conditions such as cardiovascular disease, immune system disorders, and infectious diseases; and (c) a history of substance dependence other than nicotine dependence. Two trained psychiatrists assessed all patients to determine whether they met the study criteria.

We also recruited 76 healthy individuals from the local community as a control group. Inclusion (16–65 years old) and exclusion (those with any history of psychiatric disorders, physical diseases, and substance or alcohol addiction) of healthy controls were done by trained researchers.

The study was authorized by the Ethics Committee of Ningbo Kangning Hospital and complied with the Declaration of Helsinki regarding informed consent and confidentiality. All participants voluntarily took part in this study and completed an informed consent form after the researchers introduced the study procedure in plain language.

Clinical assessment

The HDRS-24 was used to measure patients' current depressive episodes (Schwab et al., 1967). The Hamilton Anxiety Rating Scale-14 (HARS-14) was applied to assess the patients' anxiety symptoms (Hamilton, 1959). In addition, we used the Chinese edition of the Repeated Battery for the Assessment of Neuropsychological Status (RBANS) to measure participants' cognitive functioning. The RBANS consists of five subscales: immediate memory, delay memory, language, attention, and visuospatial/ constructional (Randolph et al., 1998). The scale was translated by our team and has shown adequate psychometric properties in Chinese populations (Zhang et al., 2008). In the present sample, the Cronbach's α for this scale was 0.81. Two psychiatrists conducted these assessments after consistent training for reliability, with the Intra-class Correlation Coefficient (ICC) exceeding 0.8.

TMS procedure

The left DLPFC region was stimulated with 100 single-pulse TMS stimuli by connecting a figure-of-8 coil and a TMS stimulator (Magstim Rapid2, UK). To locate the left DLPFC, the coil was fixed on F3. The handle of the splay lock was 45° slanted back and perpendicular to the scalp. The TMS pulse waveform used in our study conformed to the standard biphasic waveform typically applied in TMS studies (Supplementary Fig. S1). This waveform consisted of an initial phase followed by a second phase of opposite polarity, ensuring balanced stimulation and minimal residual charge. The direction of the induced current during stimulation is shown in Supplementary Fig. S2. The current generated by the coil flows in a specified direction, ensuring precise targeting of the cortical area. Before the TMS-EEG session, a resting motor threshold (RMT) test was performed on the left motor cortex to measure stimulus intensity. The RMT is the lowest intensity of stimulus that elicits a significant motor response in the right abductor pollicis brevis. While the participant was fully relaxed, we gradually adjusted the TMS stimulation intensity. Typically, the RMT elicited a motor-evoked potential (MEP) greater than 50 mV on at least 50 % of the trials (5 out of 10). We employed a stepwise approach, starting at a lower intensity and increasing by 1-2 % until the threshold was reached. The interval between stimuli was 5 s ± 10 % jitter, and the stimulus intensity was 110 % of RMT.

EEG recording and analysis

TMS-evoked potentials (TEPs) were recorded in a soundproofed, temperature-regulated, and electrically shielded room with a TMS-compatible 64-channel cap (Easycap, Germany) and a BrainVision Recorder (BrainProducts, Germany). The electrodes were grounded to AFz and referenced to FCz. The electrode impedance was kept lower than 5 kΩ during the recording, and the sampling rate was set to 25 kHz. During stimulation, participants kept their eyes open and wore earmuffs to prevent the appearance of associated auditory evoked potentials. We also placed a foam layer between the coil and the head to reduce noise (ter Braack et al., 2015).

TMS-EEG processing was conducted with EEGLAB (Delorme & Makeig, 2004), FieldTrip (Oostenveld et al., 2011), and customized MATLAB scripts (R2022b, The MathWorks, Inc.). Data containing TMS pulses (-5 to 15 ms) were first removed and recovered with linear interpolation. Later, the data were downsampled to 1 kHz, baseline corrected (-550 to-200 ms), and optimized for the TMS stimulus pulse (-2000 to 2000 ms). Extreme noise in the data was visually inspected (e.g., muscle movement, electrode damage, etc.) (Rogasch et al., 2013). We initially performed Fast Independent Component Analysis (FastICA) to automatically eliminate TMS tails and significant attenuation artifacts. In the first round of ICA, the average number of components removed per participant was 15.8. Bandpass filtering (1–100 Hz) was used to remove high-frequency noise and drift in the data. Then a band reject (48, 52 Hz) filter removed 50 Hz alternating current line artifacts. A second phase of ICA was conducted to remove remaining artefacts (blinks, eye movements, sustained muscle artefacts). The average number of components removed per participant in the second ICA round was 8.2. The cleared data were re-referenced to the mean for further investigation. EEG data analysis and artifact elimination were referred to Rogasch et al. (2017).

GMFA analysis

GMFA was computed for each participant using an equation developed by Lehmann and Skrandies (1980).GMFA(t)=∑ik(Vi(t)−Vmean(t))2/K

GMFA determines the largest amplitude of evoked electric fields and quantifies the whole brain's neurophysiological responses to TMS-EEG (Farzan et al., 2013). We calculated the amplitude of the GMFA component peaks for each participant. The time window for each component was determined based on previous studies (P30: 25 to 35 ms; N45: 40 to 50 ms; P60: 50 to 70 ms; N100: 80 to 120 ms; P180: 160 to 200 ms). The area under the GMFA curve (GMFA-AUC) was applied to probe the late N100-P180 complex. It was determined by summing the amplitudes 150–210 ms following the TMS pulse.

Statistical analysis

Each generated data was analyzed using SPSS for Windows (IBM Corporation, Armonk, NY, version 19.0). Descriptive statistics were used to show demographic and clinical assessment results. We used χ2 tests for categorical variables and the Kolmogorov-Smirnov test to check for the normality of continuous variables. Then, we logarithmically transformed non-normally distributed variables into normally distributed variables. Analysis of covariance (ANCOVA) was conducted to examine for between-group differences in demographic, GMFA components, and clinical variables, using years of education as a covariate. We also performed Pearson correlation analyses between GMFA components and clinical variables, including HDRS, HDRA, and RBANS. Further multiple linear regression analyses were performed to predict the clinical features of the MDD group. Specifically, these analyses used the HDRS, HDRA, and RBANS total or index scores as dependent variables and multiple variables as predictors, including TMS-EEG markers, age at onset, and years of education. The variance inflation factor (VIF) was introduced to evaluate multicollinearity between independent variables. All analyses applied a two-sided significance level (p < 0.05).

Results

Demographic and clinical characteristics

After EEG cleaning, 124 patients with MDD and 71 healthy controls were retained. Table 1 presents the demographic information and clinical results of these study participants. The healthy group had higher years of education than the MDD group (P < 0.001). Gender, age, and BMI were not considerably different between the two groups. For the clinical tests, before controlling for years of education, the MDD group had worse cognitive functioning than the healthy group concerning the RBANS total score and all subscale scores (P < 0.001). However, the difference of immediate memory was insignificant after controlling for years of education (p = 0.13, see Table 1). After Bonferroni correction, only the differences in RBANS language and total scores remained significant. The MDD group exhibited significantly higher HDRS and HARS scale scores than the healthy group (P < 0.001). The current medications of the MDD patients are summarized in Table 1.Table 1 Demographics and clinical information of healthy controls and MDD patients, and current medications of MDD patients.

Table 1	Healthy Controls
(n = 71)	MDD Patients
(n = 124)	F/χ2	p	
Demographics					
Age (years)	32.55 ± 9.33	30.58 ± 15.02	0.994	0.320	
Female	43(60.56 %)	82(67.21 %)	0.87	0.351	
Education (years)	14.42 ± 2.71	10.89 ± 3.19	60.98	<0.001	
Body Mass Index (BMI)	22.19 ± 2.86	23.03 ± 11.95	0.311	0.578	
Clinical					
RBANS score					
  Immediate Memory	92.45 ± 14.51	81.90 ± 16.24	2.246	a0.136	
  Visuospatial/ Constructional	103.43 ± 12.02	94.77 ± 15.31	4.037	a0.046	
  Language	101.48 ± 10.74	92.37 ± 14.06	8.758	a0.003	
  Attention	113.49 ± 12.75	101.18 ± 15.07	4.558	a0.034	
  Delayed Memory	96.80 ± 11.08	88.07 ± 14.97	4.484	a0.036	
Total score	101.48 ± 10.93	88.64 ± 13.29	13.237	a < 0.001	
HDRS score	4.38 ± 1.91	24.67 ± 4.89	768.15	a < 0.001	
HARS score	3.39 ± 1.39	15.60 ± 6.37	161.51	a < 0.001	
Episode Duration (months)	/	47.69 ± 55.80			
Age at Onset (years)	/	24.56 ± 13.32			
History of Suicide	/	42(34.7 %)			
Medications					
 Escitalopram	/	72			
  Fluoxetine	/	9			
  Sertraline	/	17			
  Venlafaxine	/	21			
  Duloxetine	/	5			
a ANCOVA to adjust for education years.

GMFA component analysis

Fig. 1 shows the butterfly plot of TEPs waveform from TMS-EEG for healthy controls and MDD patients, respectively. Fig. 2 depicts the mean GMFA curves of the healthy and MDD groups. The mean P180 amplitude was considerably larger in the healthy group than the MDD group (F = 5.412, p < 0.05, η2 = 0.027; Fig. 3A). GMFA-AUC was also considerably greater in the healthy group than in the MDD group (F = 5.349, p < 0.05, η2 = 0.027; Fig. 3B). However, these results were not significant after Bonferroni correction. The P30, N45, P60, and N100 components showed no considerable difference between the two groups.Fig. 1 Butterfly plot of TEPs waveform for healthy controls (A) and patients with MDD (B).

Fig 1

Fig. 2 Global mean field amplitude for healthy controls (top panel) and patients with MDD (bottom panel).

Fig 2

Fig. 3 Comparison of markers for TMS-EEG between healthy controls and patients with MDD. (A) Difference in P180 amplitude. (B) Difference in GMFA-AUC.

Fig 3

In correlation analysis, P180 amplitude was negatively correlated with HARS and HDRS scores in the MDD group (Fig. 4). In the MDD group, P30 amplitude was negatively associated with Visuospatial/ Constructional and total scores of the RBANS (Fig. 5). These results remained significant after Bonferroni correction. In the healthy group, no significant correlation was observed between GMFA components and cognitive function.Fig. 4 Correlation between P180 component and clinical symptoms in MDD patients. (A) Correlation between P180 amplitude and HDRS score. (B) Correlation between P180 amplitude and HARS score.

Fig 4

Fig. 5 Correlation between P30 component and cognitive function in MDD patients. (A) Correlation between P30 amplitude and RBANS Visuospatial/Constructional subscore. (B) Correlation between P30 amplitude and RBANS total score.

Fig 5

Multiple linear regression analyses

We further used linear regression analyses to explore the association between the GMFA components, cognitive function, and clinical symptoms in the MDD group. In terms of cognitive function, P30 amplitude (β = −0.21, t = −2.22, p = 0.02), age at onset (β = −0.36, t = −3.52, p = 0.001), and years of education (β = 0.33, t = 3.23, p = 0.002) were correlated with the total RBANS score. P30 amplitude (β = −0.28, t = −3.04, p = 0.003), age at onset (β = −0.35, t = −3.14, p = 0.002), and years of education (β = 0.25, t = 2.73, p = 0.007) were also correlated with the Visuospatial/ Constructional score. After Bonferroni correction, the association between P30 and RBANS total score was not significant.

For clinical symptoms, P180 was independently associated with HDRS score (β = −0.23, t = −2.58, p = 0.011) and HARS score (β = −0.22, t = −2.45, p = 0.016). The GMFA-AUC was independently associated with HDRS score (β = −0.22, t = −2.54, p = 0.012) and HARS score (β = −0.22, t = −2.44, p = 0.016). However, after Bonferroni correction, these results were not significant.

Discussion

In our previous study, we found that P60 was lower in MDD patients than in healthy controls and was negatively correlated with the severity of depression (Li et al., 2023). Based on these important results, we independently recruited a larger sample using GMFA measures, focusing on exploring the relationship between TMS-EEG markers, clinical symptoms, and cognitive deficits in MDD patients. As the previous study did not collect any data on cognitive deficits, we did not use the overlapping dataset.

To our knowledge, this study is the first to examine the association between TMS-EEG components and cognitive function in MDD patients. This study confirmed the abnormalities of cortical excitability in MDD patients and provided evidence for its correlation with clinical symptoms and cognitive function. The key findings of our study were as follows: (1) P180 amplitude was lower in MDD patients than healthy controls. (2) In the MDD group, P180 amplitude was negatively correlated with depressive and anxiety symptoms. (3) Patients with MDD showed significant neurocognitive deficits than the healthy group. In the MDD group, P30 amplitude was negatively correlated with RBANS Visuospatial/ Constructional and total scores.

Our study found that P180 was smaller in MDD patients compared to healthy controls. Although the physiological mechanisms of P180 are unknown, several studies have suggested that it may be modulated by axonal excitability (Premoli et al., 2017b). One study found that P180 late activity was especially susceptible to VGSC blockade (e.g., carbamazepine) (Darmani et al., 2019). A recent review has also shown that antiepileptic and excitability-lowering drugs significantly reduce P180 component amplitude in epileptic patients and healthy individuals (Gefferie et al., 2023). Indeed, the Na+ channel system plays a crucial role in glutamate release, and there is growing evidence that levels of glutamate metabolites are decreased in the prefrontal cortex and medial frontal cortex of MDD patients (Kantrowitz et al., 2021; Moriguchi et al., 2019). Lower levels of P180 may indicate abnormal glutamatergic neurotransmission in patients with MDD. In addition, GMFA-AUC was also smaller in the late component of the MDD group, suggesting reduced cortical excitability associated with P180. However, we failed to find significant results for N45, P60, and N100, which is inconsistent with previous studies (Dhami et al., 2020; Voineskos et al., 2019). The reason for this discrepancy could be differences in sample size, medication regimen and disease duration.

Furthermore, our findings showed that P180 amplitude was negatively correlated with depression and anxiety symptoms in MDD patients. Prior studies have shown that long-interval intracortical inhibition (LICI) mediated by GABAB receptors significantly reduces P180 amplitude (Premoli et al., 2014b). Based on the above evidence, P180 can be used to indicate cortical excitability. Our findings suggest that low excitability of the DLPFC is correlated with depression and anxiety symptoms in MDD patients, which aligns with recent studies (Pilisi et al., 2020; Yosephi et al., 2019). It was suggested that the similar clinical phenotypic of MDD and anxiety symptoms may depend on shared prefrontal alterations (Eleonora et al., 2019). These alterations may be related to glutamatergic and GABAergic-mediated excitation-inhibition balance. However, one study found that P180 did not respond to GABAergic drugs, suggesting that P180 may not be under the direct control of GABAergic neurons (Premoli et al., 2017a). Due to the lack of neuropharmacological studies on P180, we could not determine its neural mechanisms.

More importantly, our findings demonstrated that cognitive dysfunction was prevalent in patients with MDD and found that P30 amplitude was negatively associated with the level of cognitive functioning. The main types of cognitive dysfunction in MDD patients included attention, language, memory, and visuospatial/constructional dysfunction. Recent studies have shown that high P30 amplitude predicts cognitive and memory decline in AD patients, suggesting that the strength of connectivity between the left DLPFC and the right superior parietal cortex is associated with low cognitive function (Bagattini et al., 2019). Our findings on P30 and visuospatial functioning in patients with MMD appear to support this view, as the parietal cortex plays a crucial role in spatial cognitive functioning (Husain & Nachev, 2007). Enhanced prefrontal-to-parietal connections may be the result of a compensatory mechanism for the decline in parietal connectivity and function; however, this compensation is considered pathological and not effective in avoiding spatial cognitive deficits in patients (Bagattini et al., 2019; Pievani et al., 2014). Notably, recent research have found that short-range positive functional connectivity is reduced in the right superior parietal cortex in patients with MDD (Zhang et al., 2023). This evidence suggests that there may be some similarities in the pathology of visuospatial cognitive deficits in patients with MDD and AD.

In addition, previous studies have suggested that GABAA receptors may mediate the formation and regulation of P30 (Ferreri et al., 2011). As one of the major inhibitory neurotransmitters, GABAA controls most of the rapid inhibitory neurotransmission in the brain. Therefore, the negative correlation between P30 and cognitive function may support the idea that cognitive deficits in MDD patients are related to abnormal inhibitory mechanisms. Recent studies have highlighted that low GABAergic inhibition in the prefrontal cortex contributes largely to cognitive impairments in patients with MDD (Luscher et al., 2023). It has been suggested that enhanced dendritic inhibition via α5-GABAA receptor potentiation may have therapeutic effects in patients with memory impairment, age-related cognitive deficits, and depression (Jacob, 2019; Koh et al., 2020). We also found that age at onset was negatively correlated with cognitive function in patients with MDD, which can be linked to abnormalities in the GABA system in the ventromedial prefrontal cortex (Hasler et al., 2005). Thus, our findings may suggest that GABAergic deficits and hyperexcitability of the prefrontal cortex are related with cognitive dysfunction in MDD patients.

The current study has several limitations. First off, due to the cross-sectional design, we were unable to determine causal relationships between GMFA biomarkers and clinical variables. Second, patients were treated with antidepressant medication during the TMS-EEG test, so we cannot rule out the potential effect of antidepressant medication on these results. Future studies should include patients who have not received medication, which may help to address this potential confounder definitively. Third, in this study, patients with MDD and healthy controls were not matched in terms of sample size and education level. Future studies need tighter controls to eliminate the effects of demographic differences. Fourth, computing GMFA does not capture the polarity of the TEP components. Although some studies have used GMFA to determine the polarity of components (e.g., N45, N100) with significant results (Strafella et al., 2023; Voineskos et al., 2021, 2019), this is only an extrapolation based on the TEPs component shown in the butterfly plot and should be viewed with caution. Fifth, because GMFA was used and TMS targeted the DLPFC, the results may reflect a general difference in cortical activity between MDD patients and healthy controls. We were unable to draw conclusions about the specificity of DLPFC. Future studies need to limit the area of interest and include stimulation of a control area to better understand the role of the DLPFC. Finally, as a component of the N100-P200 complex in the TMS-EEG waveform, a portion of P180 is considered an auditory evoked activity evoked by TMS "clicks" (Conde et al., 2019). Based on this view, both groups underwent the same TMS-EEG procedure, performed with auditory masking to eliminate the effects of the 'click'. However, further studies are required to clarify the physiological basis of P180.

In summary, our study using TMS-EEG technology provides evidence for the relationship between abnormal TMS-EEG measurements, clinical symptom severity, and cognitive functioning in patients with MDD. P180 and P30 have the potential to serve as neurophysiological biomarkers of clinical symptoms and cognitive dysfunction, respectively, in MDD patients. This research also demonstrates that cortical excitability, associated with neurotransmission and cortical connectivity, is critical in the pathological process of MDD patients. Nevertheless, further studies are needed to investigate the neurophysiological mechanisms and clinical significance of P30 and P180 in MDD patients.

Funding

This study was funded by STI2030-Major Projects 2021ZD0202102 , the Basic Public Welfare Project of Zhejiang Province (LGF22H090055 ), the Medical Health Science and Technology Project of the Zhejiang Provincial Health Commission (2019RC079 ), and the CAS Key Lab of Mental Health.

CRediT authorship contribution statement

Deyang Li: Data curation, Methodology, Formal analysis, Writing – original draft. Xingxing Li: Data curation, Resources, Investigation. Jiaxin Li: Formal analysis, Data curation. Junyao Liu: Formal analysis, Data curation. Ruichenxi Luo: Formal analysis, Data curation. Yanli Li: Formal analysis, Data curation. Dongmei Wang: Conceptualization, Project administration, Writing – review & editing. Dongsheng Zhou: Resources, Supervision. Xiang-Yang Zhang: Project administration, Funding acquisition, Writing – review & editing.

Declaration of competing interest

The authors declare no conflicts of interest.

Appendix Supplementary materials

Image, application 1

Acknowledgments

The authors would like to thank all the participants who participated in this study. We would also like to thank the clinical psychiatrists for their significant contributions to this study.

Supplementary material associated with this article can be found, in the online version, at doi:10.1016/j.ijchp.2024.100495.
==== Refs
References

Bagattini C. Mutanen T.P. Fracassi C. Manenti R. Cotelli M. Ilmoniemi R.J. Predicting Alzheimer's disease severity by means of TMS-EEG coregistration Neurobiology of Aging 80 2019 38 45 10.1016/j.neurobiolaging.2019.04.008 31077959
Belardinelli P. König F. Liang C. Premoli I. Desideri D. Müller-Dahlhaus F. TMS-EEG signatures of glutamatergic neurotransmission in human cortex Scientific Reports 11 1 2021 8159 10.1038/s41598-021-87533-z 33854132
Biermann L. Wunram H.L. Pokorny L. Breitinger E. Grossheinrich N. Jarczok T.A. Changes in the TMS-evoked potential N100 in the dorsolateral prefrontal cortex as a function of depression severity in adolescents Journal of Neural Transmission 129 11 2022 1339 1352 10.1007/s00702-022-02539-9 36029418
Burke M.J. Romanella S.M. Mencarelli L. Greben R. Fox M.D. Kaptchuk T.J. Placebo effects and neuromodulation for depression: A meta-analysis and evaluation of shared mechanisms Molecular Psychiatry 27 3 2022 1658 1666 10.1038/s41380-021-01397-3 34903861
Casarotto S. Määttä S. Herukka S.K. Pigorini A. Napolitani M. Gosseries O. Transcranial magnetic stimulation-evoked EEG/cortical potentials in physiological and pathological aging Neuroreport 22 12 2011 592 597 10.1097/WNR.0b013e328349433a 21753711
Chen L.J. Wang Q. Xu T.C. Working memory function in patients with major depression disorder: A narrative review Clinical Psychology & Psychotherapy 30 2 2023 281 293 10.1002/cpp.2811 36510396
Conde V. Tomasevic L. Akopian I. Stanek K. Saturnino G.B. Thielscher A. The non-transcranial TMS-evoked potential is an inherent source of ambiguity in TMS-EEG studies NeuroImage 185 2019 300 312 10.1016/j.neuroimage.2018.10.052 30347282
Darmani G. Bergmann T.O. Zipser C. Baur D. Müller-Dahlhaus F. Ziemann U. Effects of antiepileptic drugs on cortical excitability in humans: A TMS-EMG and TMS-EEG study Human Brain Mapping 40 4 2019 1276 1289 10.1002/hbm.24448 30549127
Delorme A. Makeig S. EEGLAB: An open source toolbox for analysis of single-trial EEG dynamics including independent component analysis Journal of Neuroscience Methods 134 1 2004 9 21 10.1016/j.jneumeth.2003.10.009 15102499
Dhami P. Atluri S. Lee J.C. Knyahnytska Y. Croarkin P.E. Blumberger D.M. Prefrontal cortical reactivity and connectivity markers distinguish youth depression from healthy youth Cerebral Cortex 30 7 2020 3884 3894 10.1093/cercor/bhaa004 32118262
Dhami P. Moreno S. Croarkin P.E. Blumberger D.M. Daskalakis Z.J. Farzan F. Baseline markers of cortical excitation and inhibition predict response to theta burst stimulation treatment for youth depression Scientific Reports 13 1 2023 19115 10.1038/s41598-023-45107-1 37925557
Di Lazzaro V. Bella R. Benussi A. Bologna M. Borroni B. Capone F. Diagnostic contribution and therapeutic perspectives of transcranial magnetic stimulation in dementia Clinical Neurophysiology 132 10 2021 2568 2607 10.1016/j.clinph.2021.05.035 34482205
Eleonora M. Giuseppe D. Marika G. Marco G. Sara P. Carolina B. Common and different neural markers in major depression and anxiety disorders: A pilot structural magnetic resonance imaging study Psychiatry Research-Neuroimaging 290 2019 42 50 10.1016/j.pscychresns.2019.06.006 31279954
Farzan F. Transcranial magnetic stimulation-electroencephalography for biomarker discovery in psychiatry Biological Psychiatry 95 6 2023 564 580 10.1016/j.biopsych.2023.12.018 38142721
Farzan F. Barr M.S. Hoppenbrouwers S.S. Fitzgerald P.B. Chen R. Pascual-Leone A. The EEG correlates of the TMS-induced EMG silent period in humans NeuroImage 83 2013 120 134 10.1016/j.neuroimage.2013.06.059 23800790
Ferrarelli F. Phillips M.L. Examining and modulating neural circuits in psychiatric disorders with transcranial magnetic stimulation and electroencephalography: Present practices and future developments American Journal of Psychiatry 178 5 2021 400 413 10.1176/appi.ajp.2020.20071050 33653120
Ferreri F. Pasqualetti P. Määttä S. Ponzo D. Ferrarelli F. Tononi G. Human brain connectivity during single and paired pulse transcranial magnetic stimulation NeuroImage 54 1 2011 90 102 10.1016/j.neuroimage.2010.07.056 20682352
Fossati P. Is major depression a cognitive disorder? Revue Neurologique 174 4 2018 212 215 10.1016/j.neurol.2018.01.365 29618408
Gefferie S.R. Jiménez-Jiménez D. Visser G.H. Helling R.M. Sander J.W. Balestrini S. Transcranial magnetic stimulation-evoked electroencephalography responses as biomarkers for epilepsy: A review of study design and outcomes Human Brain Mapping 44 8 2023 3446 3460 10.1002/hbm.26260 36896753
Hadas I. Zomorrodi R. Hill A.T. Sun Y.M. Fitzgerald P.B. Blumberger D.M. Subgenual cingulate connectivity and hippocampal activation are related to MST therapeutic and adverse effects Translational Psychiatry 10 1 2020 392 10.1038/s41398-020-01042-7 33173028
Hamilton M. The assessment of anxiety states by rating The British Journal of Medical Psychology 32 1 1959 50 55 10.1111/j.2044-8341.1959.tb00467.x 13638508
Hasler G. Neumeister A. van der Veen J.W. Tumonis T. Bain E.E. Shen J. Normal prefrontal gamma-aminobutyric acid levels in remitted depressed subjects determined by proton magnetic resonance spectroscopy Biological Psychiatry 58 12 2005 969 973 10.1016/j.biopsych.2005.05.017 16043137
Husain M. Nachev P. Space and the parietal cortex Trends in Cognitive Sciences 11 1 2007 30 36 10.1016/j.tics.2006.10.011 17134935
Jacob T.C. Neurobiology and therapeutic potential of α5-GABA type A receptors Frontiers in Molecular Neuroscience 12 2019 179 10.3389/fnmol.2019.00179 31396049
Joseph S. Knezevic D. Zomorrodi R. Blumberger D.M. Daskalakis Z.J. Mulsant B.H. Dorsolateral prefrontal cortex excitability abnormalities in Alzheimer's Dementia: Findings from transcranial magnetic stimulation and electroencephalography study International Journal of Psychophysiology 169 2021 55 62 10.1016/j.ijpsycho.2021.08.008 34499960
Kallioniemi E. Daskalakis Z.J. Identifying novel biomarkers with TMS-EEG - Methodological possibilities and challenges Journal of Neuroscience Methods 377 2022 109631 10.1016/j.jneumeth.2022.109631
Kallioniemi E. Saari J. Ferreri F. Määttä S. TMS-EEG responses across the lifespan: Measurement, methods for characterisation and identified responses Journal of Neuroscience Methods 366 2022 109430 10.1016/j.jneumeth.2021.109430
Kamishikiryo T. Okada G. Itai E. Masuda Y. Yokoyama S. Takamura M. Left DLPFC activity is associated with plasma kynurenine levels and can predict treatment response to escitalopram in major depressive disorder Psychiatry and Clinical Neurosciences 76 8 2022 367 376 10.1111/pcn.13373 35543406
Kantrowitz J.T. Dong Z.C. Milak M.S. Rashid R. Kegeles L.S. Javitt D.C. Ventromedial prefrontal cortex/anterior cingulate cortex Glx, glutamate, and GABA levels in medication-free major depressive disorder Translational Psychiatry 11 1 2021 419 10.1038/s41398-021-01541-1 34354048
Koh M.T. Branch A. Haberman R. Gallagher M. Significance of inhibitory recruitment in aging with preserved cognition: Limiting gamma-aminobutyric acid type A α5 function produces memory impairment Neurobiology of Aging 91 2020 1 4 10.1016/j.neurobiolaging.2020.02.019 32240868
Kriesche D. Woll C.F.J. Tschentscher N. Engel R.R. Karch S. Neurocognitive deficits in depression: A systematic review of cognitive impairment in the acute and remitted state European Archives of Psychiatry and Clinical Neuroscience 273 5 2023 1105 1128 10.1007/s00406-022-01479-5 36048295
Lehmann D. Skrandies W. Reference-free identification of components of checkerboard-evoked multichannel potential fields Electroencephalography and Clinical Neurophysiology 48 6 1980 609 621 10.1016/0013-4694(80)90419-8 6155251
Li X.X. Chen M. Liu Q.Q. Zheng C. Yu C. Hou G.W. TMS-evoked potential in the dorsolateral prefrontal cortex to assess the severity of depression disease: A TMS-EEG study Frontiers in Pharmacology 14 2023 10.3389/fphar.2023.1207020
Luscher B. Maguire J.L. Rudolph U. Sibille E. GABAA receptors as targets for treating affective and cognitive symptoms of depression Trends in Pharmacological Sciences 44 9 2023 586 600 10.1016/j.tips.2023.06.009 37543478
Moriguchi S. Takamiya A. Noda Y. Horita N. Wada M. Tsugawa S. Glutamatergic neurometabolite levels in major depressive disorder: A systematic review and meta-analysis of proton magnetic resonance spectroscopy studies Molecular Psychiatry 24 7 2019 952 964 10.1038/s41380-018-0252-9 30315224
Morozova A. Zorkina Y. Abramova O. Pavlova O. Pavlov K. Soloveva K. Neurobiological highlights of cognitive impairment in psychiatric disorders International Journal of Molecular Sciences 23 3 2022 1217 10.3390/ijms23031217 35163141
Oostenveld R. Fries P. Maris E. Schoffelen J.M. FieldTrip: Open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data Computational Intelligence and Neuroscience 2011 10.1155/2011/156869 2011
Pievani M. Filippini N. van den Heuvel M.P. Cappa S.F. Frisoni G.B. Brain connectivity in neurodegenerative diseases-from phenotype to proteinopathy Nature Reviews Neurology 10 11 2014 620 633 10.1038/nrneurol.2014.178 25287597
Pilisi R. Tényi T. Büki A. Kovács N. Zemplényi A. Sebestyén G. The role of repetitive transcranial magnetic stimulation in the treatment of mental disorders, especially in treatment-resistant major depressive disorder Orvosi Hetilap 161 1 2020 3 10 10.1556/650.2020.31611 31884813
Poorganji M. Zomorrodi R. Hawco C. Hill A.T. Hadas I. Zrenner C. Isolating sensory artifacts in the suprathreshold TMS-EEG signal over DLPFC Scientific Reports 13 1 2023 6796 10.1038/s41598-023-29920-2 37100795
Premoli I. Bergmann T.O. Fecchio M. Rosanova M. Biondi A. Belardinelli P. The impact of GABAergic drugs on TMS-induced brain oscillations in human motor cortex NeuroImage 163 2017 1 12 10.1016/j.neuroimage.2017.09.023 28917695
Premoli I. Castellanos N. Rivolta D. Belardinelli P. Bajo R. Zipser C. TMS- EEG signatures of GABAergic neurotransmission in the human cortex Journal of Neuroscience 34 16 2014 5603 5612 10.1523/JNEUROSCI.5089-13.2014 24741050
Premoli I. Costantini A. Rivolta D. Biondi A. Richardson M.P. The effect of lamotrigine and levetiracetam on TMS-evoked EEG responses depends on stimulation intensity Frontiers in Neuroscience 11 2017 585 10.3389/fnins.2017.00585 29104528
Premoli I. Rivolta D. Espenhahn S. Castellanos N. Belardinelli P. Ziemann U. Characterization of GABAB-receptor mediated neurotransmission in the human cortex by paired-pulse TMS-EEG NeuroImage 103 2014 152 162 10.1016/j.neuroimage.2014.09.028 25245814
Randolph C. Tierney M.C. Mohr E. Chase T.N. The repeatable battery for the assessment of neuropsychological status (RBANS): Preliminary clinical validity Journal of Clinical and Experimental Neuropsychology 20 3 1998 310 319 10.1076/jcen.20.3.310.823 9845158
Redondo-Camos M. Cattaneo G. Perellon-Alfonso R. Alviarez-Schulze V. Morris T.P. Solana-Sanchez J. Local prefrontal cortex TMS-induced reactivity is related to working memory and reasoning in middle-aged adults Frontiers in Psychology 13 2022 813444 10.3389/fpsyg.2022.813444
Rogasch N.C. Daskalakis Z.J. Fitzgerald P.B. Cortical inhibition of distinct mechanisms in the dorsolateral prefrontal cortex is related to working memory performance: A TMS-EEG study Cortex 64 2015 68 77 10.1016/j.cortex.2014.10.003 25461708
Rogasch N.C. Sullivan C. Thomson R.H. Rose N.S. Bailey N.W. Fitzgerald P.B. Analysing concurrent transcranial magnetic stimulation and electroencephalographic data: A review and introduction to the open-source TESA software NeuroImage 147 2017 934 951 10.1016/j.neuroimage.2016.10.031 27771347
Rogasch N.C. Thomson R.H. Daskalakis Z.J. Fitzgerald P.B. Short-latency artifacts associated with concurrent TMS-EEG Brain Stimulation 6 6 2013 868 876 10.1016/j.brs.2013.04.004 23651674
Rogasch N.C. Zipser C. Darmani G. Mutanen T.P. Biabani M. Zrenner C. The effects of NMDA receptor blockade on TMS-evoked EEG potentials from prefrontal and parietal cortex Scientific Reports 10 1 2020 3168 10.1038/s41598-020-59911-6 32081901
Schwab J.J. Bialow M.R. Clemmons R.S. Holzer C.E. Hamilton rating scale for depression with medical in-patients British Journal of Psychiatry 113 494 1967 83 88 10.1192/bjp.113.494.83
Shen T. Li C. Wang B. Yang W.M. Zhang C. Wu Z.G. Increased cognition connectivity network in major depression disorder: A fMRI study Psychiatry Investigation 12 2 2015 227 234 10.4306/pi.2015.12.2.227 25866524
Strafella R. Momi D. Zomorrodi R. Lissemore J. Noda Y. Chen R. Identifying neurophysiological markers of intermittent theta burst stimulation in treatment-resistant depression using transcranial magnetic stimulation- electroencephalography Biological Psychiatry 94 6 2023 454 465 10.1016/j.biopsych.2023.04.011 37084864
Subhas N. Ang J.K. Tan K.A. Ahmad S.N.A. Relations between clinical characteristics and cognitive deficits among adult patients diagnosed with major depressive disorder International Journal of Psychiatry in Clinical Practice 27 3 2023 219 231 10.1080/13651501.2022.2149415 36448673
ter Braack E.M. de Vos C.C. van Putten M. Masking the auditory evoked potential in TMS-EEG: A comparison of various methods Brain Topography 28 3 2015 520 528 10.1007/s10548-013-0312-z 23996091
Tran K.H. Luki J. Hanstock S. Hanstock C.C. Seres P. Aitchison K. Decreased GABA plus ratios referenced to creatine and phosphocreatine in the left dorsolateral prefrontal cortex of females of reproductive age with major depression Journal of Psychiatry & Neuroscience 48 4 2023 E285 E294 10.1503/jpn.230016 37607825
Voineskos D. Blumberger D.M. Rogasch N.C. Zomorrodi R. Farzan F. Foussias G. Neurophysiological effects of repetitive transcranial magnetic stimulation (rTMS) in treatment resistant depression Clinical Neurophysiology 132 9 2021 2306 2316 10.1016/j.clinph.2021.05.008 34167891
Voineskos D. Blumberger D.M. Zomorrodi R. Rogasch N.C. Farzan F. Foussias G. Altered transcranial magnetic stimulation-electroencephalographic markers of inhibition and excitation in the dorsolateral prefrontal cortex in major depressive disorder Biological Psychiatry 85 6 2019 477 486 10.1016/j.biopsych.2018.09.032 30503506
White L.K. Makhoul W. Teferi M. Sheline Y.I. Balderston N.L. The role of dlPFC laterality in the expression and regulation of anxiety Neuropharmacology 224 2023 109355 10.1016/j.neuropharm.2022.109355
World Health Organization. (2023). World Health Organization fact sheet of depressive disorder. Retrieved December 11 2023 from https://www.who.int/news-room/fact-sheets/detail/depression.
Yosephi M.H. Ehsani F. Daghiani M. Zoghi M. Jaberzadeh S. The effects of trans-cranial direct current stimulation intervention on fear: A systematic review of literature Journal of Clinical Neuroscience 62 2019 7 13 10.1016/j.jocn.2019.01.011 30660476
Zhang B.H. Tan Y.L. Zhang W.F. Wang Z. Yang R. G G. Repeatable battery for the assessment of neuropsychologicar status as a screening test in Chinese: Reliability and validity Chinese Mental Health Journal 22 12 2008 865 869
Zhang L.L. Cui X.L. Ou Y.P. Liu F. Li H.B. Xie G.J. Abnormal long- and short-range functional connectivity in patients with first-episode drug-naive melancholic and non-melancholic major depressive disorder Journal of Affective Disorders 320 2023 360 369 10.1016/j.jad.2022.09.161 36206876
Zhang X.M. Zhang R.R. Lv L.L. Qi X.Y. Shi J.P. Xie S.P. Correlation between cognitive deficits and dorsolateral prefrontal cortex functional connectivity in first-episode depression Journal of Affective Disorders 312 2022 152 158 10.1016/j.jad.2022.06.024 35752217
