Next Article in Journal
Acute Acalculous Cholecystitis Associated with Abscesses—An Unknown Dual Pathology
Next Article in Special Issue
Characterization of Vagus Nerve Stimulation (VNS) Dose-Dependent Effects on EEG Power Spectrum and Synchronization
Previous Article in Journal
A Preliminary Study on Microbiota Characteristics of Bronchoalveolar Lavage Fluid in Patients with Pulmonary Nodules Based on Metagenomic Next-Generation Sequencing
Previous Article in Special Issue
Transcutaneous Cervical Spinal Cord Stimulation Combined with Robotic Exoskeleton Rehabilitation for the Upper Limbs in Subjects with Cervical SCI: Clinical Trial
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Online Left-Hemispheric In-Phase Frontoparietal Theta tACS Modulates Theta-Band EEG Source-Based Large-Scale Functional Network Connectivity in Patients with Schizophrenia: A Randomized, Double-Blind, Sham-Controlled Clinical Trial

1
Department of Psychiatry, Tri-Service General Hospital, National Defense Medical Center, Taipei 114202, Taiwan
2
Department of Life Sciences, National Central University, Taoyuan 320317, Taiwan
3
Department of Nuclear Medicine, College of Medicine, The Catholic University of Korea, Seoul 07345, Republic of Korea
4
Department of Psychology, Seoul National University, Seoul 08826, Republic of Korea
5
Department of Emergency Medicine, Tri-Service General Hospital, National Defense Medical Center, Taipei 114202, Taiwan
6
Department of Emergency Medicine, Taoyuan Armed Forces General Hospital, Taoyuan 325208, Taiwan
7
Department of Psychiatry, Tri-Service General Hospital Beitou Branch, National Defense Medical Center, Taipei 112003, Taiwan
*
Author to whom correspondence should be addressed.
Biomedicines 2023, 11(2), 630; https://doi.org/10.3390/biomedicines11020630
Submission received: 3 January 2023 / Revised: 13 February 2023 / Accepted: 16 February 2023 / Published: 20 February 2023
(This article belongs to the Special Issue Neuromodulation: From Theories to Therapies)

Abstract

:
EEG studies indicated that schizophrenia patients had increased resting-state theta-band functional connectivity, which was associated with negative symptoms. We recently published the first study showing that theta (6 Hz) transcranial alternating current stimulation (tACS) over left prefrontal and parietal cortices during a working memory task for accentuating frontoparietal theta-band synchronization (in-phase theta-tACS) reduced negative symptoms in schizophrenia patients. Here, we hypothesized that in-phase theta-tACS can modulate theta-band large-scale networks connectivity in schizophrenia patients. In this randomized, double-blind, sham-controlled trial, patients received twice-daily, 2 mA, 20-min sessions of in-phase theta-tACS for 5 consecutive weekdays (n = 18) or a sham stimulation (n = 18). Resting-state electroencephalography data were collected at baseline, end of stimulation, and at one-week follow-up. Exact low resolution electromagnetic tomography (eLORETA) was used to compute intra-cortical activity. Lagged phase synchronization (LPS) was used to measure whole-brain source-based functional connectivity across 84 cortical regions at theta frequency (5–7 Hz). EEG data from 35 patients were analyzed. We found that in-phase theta-tACS significantly reduced the LPS between the posterior cingulate (PC) and the parahippocampal gyrus (PHG) in the right hemisphere only at the end of stimulation relative to sham (p = 0.0009, corrected). The reduction in right hemispheric PC-PHG LPS was significantly correlated with negative symptom improvement at the end of the stimulation (r = 0.503, p = 0.039). Our findings suggest that in-phase theta-tACS can modulate theta-band large-scale functional connectivity pertaining to negative symptoms. Considering the failure of right hemispheric PC-PHG functional connectivity to predict improvement in negative symptoms at one-week follow-up, future studies should investigate whether it can serve as a surrogate of treatment response to theta-tACS.

1. Introduction

Schizophrenia is a chronic, debilitating mental illness. Its core components include positive symptoms, negative symptoms, and cognitive deficits. Primary negative symptoms (alogia, asociality, amotivation, anhedonia, and affective flattening) account for a large part of poor functional outcomes and long-term disability in patients with schizophrenia [1]. Efficacious treatments for negative symptoms remain an important unmet need in schizophrenia.
Recent research indicates that transcranial electrical stimulation (tES) is a representative non-invasive brain stimulation method that may serve as an effective therapy in counteracting negative symptoms [2]. Common tES approaches are transcranial direct current stimulation (tDCS), transcranial alternating current stimulation (tACS), and transcranial random noise stimulation (tRNS). Our earlier research demonstrated the efficacy of adjunct tDCS and tRNS over the prefrontal cortex in improving negative symptoms of schizophrenia [3,4]. Results of a small case series indicated that tACS applied at theta frequency over the prefrontal cortex and delivered during a resting state can improve negative symptoms of schizophrenia [5]. tACS applies a low-intensity (0.5 to 4 mA) sinusoidal (frequency: 0.1 to 80 Hz; phase: 0–360 degrees; DC offset: with/without) electrical currents to the targeted brain regions through the specific electrode montages on the scalp. tACS can enhance local oscillatory activity (i.e., an increase in the power of endogenous brain oscillations in the range of the stimulated frequency) and entrain (i.e., an alignment of local brain endogenous oscillations and the tACS phase) the brain rhythms of the targeted brain region [6]. The sustained aftereffects of a single session of tACS can last up to 70 min after stimulation [7].
Functional magnetic resonance imaging (fMRI) studies identified working memory (WM) related activities in a widely distributed neural network involving the dorsolateral prefrontal cortex (DLPFC), medial prefrontal cortex (mPFC), insula, hippocampus, posterior parietal cortex (PPC), and occipito-temporal cortex [8]. Data from simultaneous fMRI-EEG scans during the retention period of a WM task demonstrated negative correlations between the EEG theta (5–7 Hz) frequency band and the blood oxygen level dependency (BOLD) signal, as well as a strong anatomical overlap in the default mode network (DMN) for the BOLD activation pattern and theta EEG-BOLD signal correlations, which were not present in other frequency domains [9]. A growing body of evidence suggests theta (~6 Hz) phase synchronization (i.e., theta coupling at ~0° relative phase) of distant cortical areas (e.g., frontoparietal regions) as vital neural mechanisms in WM performance, requiring the large-scale integration of functional cortical circuits [10]. For example, research indicated a significant increase of theta (~6 ± 1 Hz) phase synchronization between left DLPFC and PPC occurring ~200–500 ms after memory probe onset while the participants performed a WM task [10]. Polania et al. used tACS simultaneously applied at 6 Hz over the left prefrontal and parietal cortices with a relative 0° (‘‘synchronized’’ or “in-phase” condition) or 180° (‘‘desynchronized’’ or “anti-phase” condition) phase difference versus the placebo stimulation, while healthy individuals performed a working memory task (i.e., online stimulation) [10]. The results indicated that as compared to the placebo stimulation, exogenously induced frontoparietal theta synchronization by in-phase theta-tACS significantly improved cognitive performance, and exogenously induced frontoparietal theta desynchronization by anti-phase theta-tACS deteriorated performance. In our recent double-blind randomized sham control trial, we applied in-phase tACS tuned at theta frequency (6 Hz) and delivered this while the participants performed a working memory task (online theta-tACS), aiming to artificially induce coupling behaviorally relevant brain rhythms between left frontal and parietal regions. The results indicated that 10 sessions of online theta-tACS significantly improved negative symptoms (primary outcome) as well as working memory performance (secondary outcome) in stabilized patients with schizophrenia [11]. It remains unknown whether the repeated application of online theta-tACS can alter or modulate resting-state cortical networks in patients with schizophrenia.
Research indicates that a fundamental pathophysiological feature of schizophrenia consists of a functional disconnection syndrome, in which neural oscillations at the theta frequency are robustly implicated [12]. For example, evidence from studies of patients with schizophrenia indicates increased theta activity during resting conditions and reduced activity for targets during active paradigms [13]. Excessive theta activity (4–8 Hz) during resting conditions is involved in impaired cognitive processes (e.g., working memory, processing speed, verbal learning, reasoning, and problem-solving) [14,15] and negative symptom burden [16] in patients with schizophrenia. Research that performed EEG-based source localization analyses to investigate the resting-state neural connectivity in patients with schizophrenia as compared to healthy controls reported increased functional connectivity between brain networks in the theta frequency [17,18] and an increase in the synchronization between the posterior cingulate cortex, cuneus, and precuneus [19]. Moreover, the increased theta-band functional connectivity between brain networks was associated with a delay in the initiation of information processing [19] and a worse processing speed and verbal memory [17,18]. A recent study of drug naïve first-episode psychosis patients indicated that theta-band connectivity between right frontal and parietal regions underlying the active detection of auditory novelty stimuli during an auditory involuntary attention task was inversely correlated with negative symptoms [20].
Motivated by the above-mentioned findings, suggesting long-range EEG phase synchronization in a narrow theta band of 5–7 Hz as an underlying neural mechanism of WM related activities [10], this study aimed to investigate whether applying 10 sessions of in-phase theta (6 Hz)-tACS for externally inducing left frontoparietal synchronization during a WM task would modulate the resting-state theta-band (5–7 Hz) EEG source-based large-scale functional network connectivity in patients with schizophrenia and determine if the characteristics of network connectivity serve as predictors or surrogates of treatment response. We did not construct specific hypotheses regarding the direction of changes in theta-band network connectivity due to the paucity of relevant previous studies.

2. Materials and Methods

2.1. Participants

The present study is a double-blind, randomized, sham-controlled trial, approved by the ethics committee of Tri-Service General Hospital, Taipei, Taiwan and registered (ClinicalTrials.gov (accessed on 15 February 2023), ID: NCT04545294; Table S1). The inclusion criteria were: (1) patients aged 20–65 with DSM-5-defined schizophrenia or schizoaffective disorder; (2) duration of illness >2 years; (3) being clinically stable and on an adequate therapeutic dose of antipsychotics for at least 8 weeks before enrolment; (4) agreement to participate in the study and provide written informed consent. The exclusion criteria were (1) unstable medical conditions, current psychiatric comorbidity, prominent mood symptoms, or active substance use disorders (with the exception of caffeine and/or tobacco); (2) a history of seizures, meningitis, or encephalitis; (3) contraindications for transcranial electrical and magnetic stimulation; (4) a history of intracranial neoplasms or surgery, or a history of severe head injuries or cerebrovascular diseases; (5) pregnancy or breastfeeding at enrollment; (6) scalp skin lesions at the area of electrode application.

2.2. Study Design

This randomized, double-blind, and sham-controlled clinical trial was performed between August 2019 and April 2020. Each participant was randomly assigned to one of two groups (i.e., tACS: sham = 1:1). A study coordinator not involved in the trial created 5-digit random numbers and assigned the randomization numbers to the participants to maintain blinding integrity (Supplementary Materials). The participants and researchers were unaware of group assignment until the unblinding of the trial. The primary outcome, the change in the negative symptom severity over time, was measured by the negative symptom subscale score of the Positive and Negative Syndrome Scale (PANSS) at baseline, end of stimulation, and at 1-week and 1-month follow-up visits (Figure S1 in Supplementary Materials). Detailed analyses of clinical data have been reported in a separate paper [11].

2.3. Brain Stimulation

Online in-phase theta-tACS was applied during the performance of a dual n-back task involving both visuospatial and auditory-verbal WM paradigms (see Supplementary Materials). Sinusoidal tACS was delivered by two battery-operated devices (Eldith DC stimulator Plus, NeuroConn, Ilmenau, Germany) connected with two 4 × 1 wire adaptors (Equalizer Box, NeuroConn), via 10 carbon rubber electrodes (1 cm radius, high-density 4 × 1 configuration with a gel layer of 2.0 mm), at 6 Hz frequency, 2 mA current intensity without DC offset, with 100 cycles ramp-up/ramp-down and a 0° relative phase, for 20 min, twice-daily on 5 consecutive weekdays. In the sham group, short continuous currents without neuromodulatory effects (i.e., 2 mA normal-like stimulation for the initial 30 s followed by a tiny current pulse of 110 μA over 15 ms for impedance control taking place every 550 ms during the remaining time) were applied to mimic real-stimulation sensations. Stimulation electrodes of one DC stimulator were placed at the international 10-10 electrode position F1, F5, AF3, and FC3 while the return electrode was placed at CPz. Stimulation electrodes of the other DC stimulator were placed at P1, P5, CP3, and PO3 while the return electrode was placed at FCz (Figure S2 in Supplementary Materials). These electrodes were held in place with conductive paste to keep the combined impedance of all electrodes <15 kΩ. A custom-made pulse generator controlled the two stimulators to achieve a synchronous (in-phase) setup (0° relative phase difference between the output signals of the two stimulators). The chosen montage was adapted from Polania et al. [10] for theta-tACS to externally induce left frontoparietal synchronization. HD-Explore® (Soterix Medical, New York, NY, USA) was used to simulate the numerical computation of the electric field of the current montage (Figure 1).

2.4. Electrical Source Estimation of Resting-State EEG

Resting-state EEG (rsEEG) was collected at baseline, end of stimulation, and one-week follow-up, using a 32-channel EEG cap (NP32, GmbH, Ilmenau, Germany) with Ag/AgCl sintered ring electrodes placed according to the international 10–20 system and referenced to the tip of the nose, together with Neuro Prax® TMS/tES compatible full band DC-EEG system (NeuroConn GmbH, Ilmenau, Germany) with a sampling frequency of 4000 Hz, an analogue-digital precision of 24 bits, and an analogous bandpass filter (0–1200 Hz). Patients were seated comfortably in a recliner in a light and sound-attenuated room. They were instructed not to drink caffeinated beverages one hour prior to, and alcohol 24 h prior to EEG recording to avoid caffeine- or alcohol-induced changes in the EEG stream. The ground electrode was placed at Fpz. The horizontal electrooculogram (HEOG) was recorded by two electrodes placed at 1 cm from the outer canthi of both eyes. Two electrodes were placed above and below the left eye, respectively, to record the blinks and vertical electrooculogram (VEOG). The impedance of each electrode was checked to remain below 5 kΩ. Before starting the EEG recording, the patients performed 3-min calibration tasks to estimate the influence of horizontal/vertical movements and blink artifacts on EEG, which was processed and stored in the Neuro Prax® EEG system built-in software providing fully automatic correction of real-time EEG artifacts caused by blinking and eye or body movement during the subsequent EEG recording (Supplementary Materials). rsEEG with eyes open (5 min) and eyes closed (5 min) were recorded for a total of 10 min and the sequence was randomized and counterbalanced across the patients. Patients were instructed to visually fixate on a crosshair in front of them during the eyes-open condition or stay relaxed in a state of mind wandering (i.e., without goal-oriented mental activity) with their eyes closed during the eyes-close condition. EEGLAB v2020.0 [21], an interactive Matlab (MathWorks, Natick, MA, USA) toolbox was used for preprocessing the offline EEG data, e.g., data being downsampled to 500 Hz, band-pass filtered to 1–100 Hz with the finite impulse response method and analog 60 Hz-notch filtered. Bad channels were automatically detected and removed based on artifact subspace reconstruction (ASR) [22]. Independent component analysis (ICA) followed by ICLabel [23] was used to automatically remove artifacts caused by muscle activity, heartbeats, eye movements, and eye blinks (see Supplementary Materials). Considering that theta-tACS was applied in an eyes-open state, only accepted epochs of eyes-open EEG data collected in a resting state were selected for electrical source estimation. All participants had one artifact-free epoch with a length of a minimum of 60 s, respecting the guidelines from previous research [24].
All source imaging was performed with the exact low-resolution brain electromagnetic tomography (eLORETA), a linear inverse solution used to compute the standardized current density or distribution of current density across voxels in the brain by localizing and reconstructing the intracerebral electrical sources underlying the scalp-recorded activity [24] in a realistic head model [25], using the Montreal Neurological Institute (MNI; Montreal, Quebec, Canada) MNI152 template [26], with the three-dimensional eLORETA inverse solution space (i.e., intracerebral volume) restricted to cortical gray matter and hippocampi, as determined by the probabilistic Talairach atlas [27], and partitioned in 6239 voxels (voxel size 5 × 5 × 5 mm).

2.5. Whole-Brain Electrical Source-Based Functional Connectivity

In this study, whole-brain EEG source-based functional connectivity analyses were performed by the ‘‘whole-brain Brodmann areas (BAs)’’ approach, using the anatomical definitions of 84 BAs provided by eLORETA software package and based on the Talairach Daemon “http://www.talairach.org/ (accessed on 15 February 2023)” [28]. Standard electrode positions on the MNI152 scalp were used and the intracerebral volume’s partitioning was restricted to 6239 voxels of 5 × 5 × 5 mm spatial resolution. The electric activity at each voxel in the neuroanatomic MNI space as the exact magnitude of the estimated current density is represented by the eLORETA images. We used Brodmann areas as anatomical labels using MNI space with the correction to Talairach space. We chose a ROI-maker 2 method (available in eLORETA) for the construction of the regions of interest (ROIs). Brain activity was computed with eLORETA in 42 regions, positioning the center in Brodmann Areas (BAs: 1–11, 13, 17–25, 27–47) for each hemisphere (see Table S2 for MNI coordinates of each seed), to obtain a topographic view. Pairs of BAs were analyzed using the values of single voxel with the highest F-ratio value at the centroid of each BA. Analyses of whole-brain functional connectivity were based on lagged phase synchronization (LPS) and were performed by using all 42 BAs in each hemisphere as regions of interest (ROIs) to test interregional functional correlations between any pair of BAs LPS. Thus, connectivity analyses were based on (84 × 83/2=) 3486 pairs of sources distributed throughout the cortex, sufficient for obtaining detailed estimates of connectivity within the constraints posed by the limited spatial resolution of EEG. LPS is a method associated with nonlinear functional connectivity and represents the connectivity of two signals after excluding the instantaneous zero-lag component of phase synchronization caused by intrinsic artifacts or non-physiological effects [24]. A value of 1 indicates perfect synchronization and a value of 0 indicates no synchronization. LPS between 84 ROIs was computed for each artifact-free EEG segment in the frequency domain using normalized Fourier transforms. The data in the theta frequency range (5–7 Hz) were selected for statistical analyses. For additional details see Supplement Materials.

2.6. Statistical Analyses

Statistical analyses were performed using either IBM SPSS Statistics 21.0 software (IBM SPSS Inc., Chicago, IL, USA) or the implemented statistical eLORETA nonparametric mapping (SnPM) tool [24]. Repeated-measures analyses of variance were used to assess the effects of theta-tACS on the PANSS negative symptom subscale score over time, with “time” as the within-group factor and “treatment group” as the between-group factor. When significant “time” × “treatment group” interaction was found, the post-hoc Student’s t-tests were used to compare the between-group differences at post-baseline visits. The SnPM analysis tool includes a correction for multiple comparisons. The statistical analyses of between-group changes in theta-band (5–7 Hz) electrical source estimation and source functional connectivity from baseline to post-baseline visits were conducted using t-tests that were corrected for multiple comparisons using a non-parametric permutation procedure (5000 randomizations). Spearman rank correlations were used to analyze the relationships between the significant changes in EEG-based source functional connectivity from baseline to post-baseline visits and treatment response to in-phase theta-tACS. Statistical significance for the results was set at p < 0.05 (two-tailed) and the false discovery rate (FDR) was used for multiple comparisons corrections.

3. Results

Thirty-six patients were randomly allocated to either receive in-phase theta-tACS (n = 18) or a sham stimulation (n = 18) (Figure S3: CONSORT Flowchart) and all of them completed all 10 sessions of the respective arm. The current study analyzed 35 patients with complete EEG data. Table 1 shows the concise demographics and clinical assessment of the participants. There were no significant differences in the sociodemographic and clinical characteristics at baseline between the theta-tACS group and sham group. A significant group-by-time interaction for the PANSS negative subscale score was found (F = 3.05, p = 0.0027). Post-hoc analyses showed significant between-group differences at all post-baseline assessments (Figure 2). The negative symptom severity significantly improved at the end of theta-tACS relative to sham and the effect was maintained at the follow-up visits.

3.1. Effects of Theta-tACS on Whole-Brain EEG Source-Based Theta-Band Functional Connectivity

SnPM showed no significant difference in eLORETA seed-based whole-brain analyses of functional connectivity in the theta-band frequency between theta-tACS and sham group at baseline. As can be seen in Figure 3, rows and columns in these matrices represent the 84 ROIs. The element of the matrix was the t value of the independent t test that determines whether the mean change for LPS of any pair of 84 ROIs is significantly different from baseline to postbaseline visits (i.e., end of stimulation and 1-week follow-up) among participants treated with theta-tACS relative to the sham stimulation (i.e., left panels of Figure 3A,B). To avoid ambiguity, the elements of the matrix (i.e., t values) greater than the t-value threshold for statistical significance returned 1 and those less than the threshold returned 0 in the map (i.e., right panels of Figure 3A,B). There was a significant between-group change from baseline to end of stimulation in theta-band functional connectivity (only between one pair of ROI 66 and ROI 71, see Figure 3A, right panel). ROI 66 and ROI 71 are posterior cingulate and parahippocampal gyrus, respectively (Table S2). The between-group change from baseline to one-week follow-up in theta-band functional connectivity was non-significant (Figure 3B, right panel).
To be more specific, the LPS between a region in the posterior cingulate (PC, xyz = 5, −50, 5; BA 29; right hemisphere; Figure 4A) and a region in the parahippocampal gyrus (PHG, xyz = 15, 0, −20; BA 34; right hemisphere) was significantly reduced in theta-tACS group compared with the sham (p value = 0.0009, corrected; Figure 4B). The inclusion of the antipsychotic medication dose (in olanzapine equivalents) did not substantially alter the results. In addition, the effects of theta-tACS on functional connectivity were restricted only to the specific theta-frequency (5–7 Hz) band because the results were non-significant for the broader theta band 4.5–7.5 Hz.

3.2. Correlation Analyses

The reduction in the LPS between the PC and the PHG from baseline to the end of stimulation was significantly correlated with the improvement of negative symptoms at the end of stimulation (r = 0.503, p = 0.039, Figure 5) in theta-tACS group. The inclusion of the antipsychotic medication dose (in olanzapine equivalents) did not substantially alter the results.

4. Discussion

To our knowledge, this is the first randomized, double-blind, sham-controlled clinical trial providing evidence that adjunct online in-phase theta (6 Hz)-tACS improves negative symptoms of schizophrenia through modulating theta-band (5–7 Hz) EEG source-based large-scale functional network connectivity. Recent research showed that a single session of online in-phase theta (6 Hz)-tACS with a low current intensity (1 mA) did not result in any significant change in resting-state EEG scalp-level theta-band connectivity between prefrontal and parietal areas in healthy individuals [29]. There was a lack of evidence for the acute and longer-lasting effects of 10 sessions of online in-phase theta (6 Hz)-tACS using a protocol comprising of 2 mA for 20 min (i.e., a total stimulation duration of 200 min) on resting-state EEG source-level theta-band long-range functional connectivity in patients with schizophrenia. In the present study, online in-phase theta (6 Hz)-tACS was applied to externally induce left frontoparietal synchronization during a working memory task. The results of clinical data showed that online in-phase theta (6 Hz)-tACS effectively improved negative symptoms in stable patients with schizophrenia [11] while the EEG results showed that resting-state theta-band (5–7 Hz) functional connectivity between the posterior cingulate and the parahippocampal gyrus in the right hemisphere (Figure 4) was significantly reduced at the end of stimulation among participants treated with in-phase theta (6 Hz)-tACS relative to the sham condition.
Evidence indicates that patients with schizophrenia showed increased EEG-based resting-state functional connectivity at lower frequencies (i.e., hyper-synchronizations in the theta and delta frequencies) compared to the healthy controls, especially higher right fronto-parietal connectivity and higher connectivity between central and parietal areas in the theta band [17], the increased theta-band resting-state connectivity across the midline, sensorimotor, orbitofrontal regions, and the left temporoparietal junction [18], and abnormally strengthened coupling between sources located within the posterior hub of the DMN, containing PC, precuneus, and cuneus [19]. However, these cross-sectional case-control studies cannot determine whether increased theta-band long-range functional connectivity is causal to clinical manifestations of the illness, functionally compensatory processes, or purely epiphenomena.
Our study demonstrated that adjunct 10 sessions of in-phase theta (6 Hz)-tACS during a working memory task resulted in theta-band (5–7 Hz) interregional de-synchronization between two nodes in the posterior part of DMN (i.e., the PC and the PHG) at the end of stimulation along with an improvement in negative symptoms. Furthermore, the reduction in resting-state theta-band connectivity within the above network over time was positively correlated with negative symptom improvement at the end of stimulation, thus supporting the value of right hemispheric PC-PHG LPS as a surrogate biomarker for the clinical efficacy of online in-phase theta-tACS in treating negative symptoms of schizophrenia.
A large number of fMRI/DTI studies indicate that functional hyperconnectivity within the DMN may be the most common finding in comparisons of schizophrenia patients with healthy individuals [30,31] and suggest that increased functional connectivity of the intrinsic networks may be an endophenotype or part of the core pathophysiology of schizophrenia and underlie risks for this serious mental illness [30]. A growing body of evidence also suggests that functional hyperconnectivity within the DMN may play a vital role in the manifestation of negative symptoms of schizophrenia. For example, a previous study indicated that DMN areas with higher activity or connectivity (e.g., frontal polar cortex DMN hyperactivity) contributed to regional functional pathology in schizophrenia and bore significance for negative symptoms [32]. A recent study of neuroimaging coupled with machine learning in patients with schizophrenia-spectrum disorder indicated that resting-state connectivity within the DMN measured at baseline accurately predicted the change in negative symptom severity (accuracy: 83%) at 1-year follow-up, providing its utility in guiding long-term prognostication in schizophrenia [33]. Specifically, the increased connectivity of a left-lateralized sub-network of DMN with direct involvement of hub regions, such as between the posterior cingulate and anterior medial prefrontal cortex (amPFC) and between ventromedial prefrontal cortex (vmPFC) and parahippocampal cortex, predicted the worsening of negative symptoms at 1-year follow-up. A proposed explanation for these findings was based on functional segregation within DMN regions and the putative introspective and extrospective modes of the DMN (i.e., the DMN toggling between an introspective self-referential mode and an extrospective mode that remains alert to changes in the external environment, while normally the internalization and self-reflective thinking under normal conditions are frequently interrupted and shifted toward an externalized state vigilant to external stimuli).
The DMN can be separated into two subsystems (i.e., dorsal medial and medial temporal subsystems), tightly linked by functional hubs also known as a midline core (i.e., amPFC, PC, and angular gyrus, which activate preferentially for self-relevant conditions) [34]. The dorsal medial subsystem comprises the dorsal medial prefrontal cortex (dmPFC), lateral temporal cortex, temporoparietal junction (TPJ) and temporal pole, which activate when subjects reflect on their present mental states. The medial temporal subsystem comprises the vmPFC, hippocampal formation, PHC, retrosplenial cortex, and posterior inferior parietal lobule, which are related to the subjects’ autobiographical memory and future simulations. Hyper-connectivity within the DMN functional hubs or subsystems may contribute to more internalization and fewer shifts to the externalized mode, which just characterizes the internalizing negative symptoms of schizophrenia [33].
Our results suggest a potential relationship between theta-tACS inducing a decrease in right hemispheric PC-PHG functional connectivity in the theta frequency (5–7 Hz) and improvements in negative symptoms, suggesting that theta-band long-range functional hyper-connectivity may be a neural signature of negative symptoms of schizophrenia, and theta-tACS controlled negative symptoms by effectively regulating the PC-PHG connectivity within DMN. It should be noted that a recent fMRI study indicated that hyperconnectivity within the DMN nodes, i.e., the functional connectivity between angular gyrus (AG) and middle temporal gyrus (MTG), was further enhanced in medication-resistant schizophrenia patients receiving a regular course of bitemporal electroconvulsive therapy (ECT) combined with antipsychotics for 4 weeks, as compared to those receiving only antipsychotics, suggesting that ECT did not reverse the original state of hyperconnectivity within the DMN, but provoked a compensatory increase in DMN connectivity [35]. In addition, the increases in functional connectivity between left AG and right MTG and between right AG and right MTG were positively correlated with the reduction in negative symptoms. This study implied that the hyperconnectivity of the DMN in medication-resistant schizophrenia patients reflected a compensatory mechanism to counteract the psychopathological symptoms and also highlighted the need for a separate consideration of stages of illness and various phenotypes. Specifically, evidence indicates that increased functional connectivity or an early compensatory activity is more prominent in early-stage schizophrenia, whereas decreased functional connectivity or decompensation is typical of late-stage schizophrenia [36]. Since our study recruited patients with stabilized schizophrenia and used EEG as a neuroimaging modality that allows for studying neuronal activity in different frequencies. The interpretation of our EEG-based results regarding the theta frequency-specific change of functional connectivity between nodes within DMN after theta-tACS warrants further discussion.
A previous study used magnetoencephalography (MEG) to investigate resting-state brain oscillation and the DMN based on a source space in different frequency bands (i.e., theta, alpha, beta, and gamma) and found that the spatial distribution of DMN activity in the alpha frequency was similar to that reported in fMRI studies, and schizophrenia patients had an increased resting activity of the posterior hub of DMN (i.e., PC), predominantly in the theta band [37]. It is known that theta-band oscillations are generated by the septum pellucidum, medial septal nucleus, and especially the nucleus of the vertical limb of the diagonal band of Broca, which has reciprocal connections with the hippocampal formation (e.g., PHG and surrounding structures such as PC and entorhinal cortex) and is responsible for the generation of theta waves in the brain regions [38]. The PC activity is involved in the processing of information regarding the self (e.g., during any task related to the self and others, remembering the past, and thinking about the future). The activity of the hippocampal formation including PHG is involved in spatial and scene recognition and future simulations. Both PHG and PC play an important role in internalizing negative symptoms of schizophrenia [33] and are also key theta rhythm generators in the brain. The over-synchronization between the two DMN structures may reflect the aberrant and dysfunctional reciprocal connections and represent a primary abnormality underlying the negative symptoms. Taken together, theta-tACS intervention in schizophrenia could be working by desynchronizing the neuronal networks whose over-synchronization accounts for the manifestation of negative symptoms. Future studies are needed to investigate the feasibility of using novel machine learning classifiers (e.g., deep convolutional neural network or fuzzy classifier) for pre- and post-treatment resting-state EEG connectivity to identify schizophrenia patients responding to theta-tACS [39,40].
Our study has a few limitations. First, EEG directly measures fast neural network dynamics with high temporal resolution (ms) for investigating the high temporal dynamics of the whole-brain network functional connectivity but it is limited by its ability to detect sources of electrical activity in deep structures (e.g., thalamus and cerebellum) that were reported to have aberrant connectivity in schizophrenia patients [41], as well as by its spatial resolution recorded using the 32-channel array of scalp electrodes for EEG source localization [42]. Our findings await further investigation using simultaneous EEG-fMRI for the analysis of the brain networks’ functional connectivity [43]. Second, the efficacy of theta-tACS in reducing negative symptom severity was maintained at 1-week follow-up visit, but the between-group difference in the change of right hemispheric PC-PHG functional connectivity disappeared at that timepoint. Therefore, we cannot exclude the possibility that the effect of theta-tACS on right hemispheric PC-PHG functional connectivity was a random finding. We also consider it necessary to exclude the possibility that the effect was purely an epiphenomenon from the impact of theta-tACS on the connectivity of other important large-scale networks implicated in schizophrenia due to the inherent limitations of the eLORETA seed-based approach for functional connectivity analysis with pre-defined seeds or ROIs. Third, it is worth noting that atypical antipsychotics may have a positive impact on the modulation of DMN connectivity [44,45] and thus we cannot exclude the possibility that the changes observed in PC-PHG functional connectivity were mediated by the effects of the interaction between theta-tACS and the antipsychotic medications that the participants were currently receiving. However, the direct effects of antipsychotics seem unlikely in our study since olanzapine equivalent dose had no significant effect on the results. Finally, our study did not include healthy individuals as a control group and thus the interpretation of the statement of theta-band hyper-connectivity in our patients with schizophrenia should be treated with caution. Further studies including both patients and healthy subjects are required to confirm our results.

5. Conclusions

Our study indicates that adjunct theta-tACS might contribute to the improvement of negative symptoms by modulating theta-band functional connectivity between spatially separated brain structures within the DMN in patients with schizophrenia. Furthermore, the change in right hemispheric PC-PHG functional connectivity in response to theta-tACS may serve as a biomarker to track the therapeutic effect of this novel non-invasive brain stimulation. Our findings may be a new starting point for the development of individualized and effective brain stimulation therapies in treating negative symptoms of schizophrenia.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biomedicines11020630/s1, Figure S1: Timeline for treatment and assessments; Figure S2: The 2D head model of left-hemispheric in-phase theta-rhythm frontoparietal transcranial alternating current stimulation (theta-tACS); Figure S3: A CONSORT Flow Diagram of this clinical trial; Table S1: A CONSORT checklist of this clinical trial; Table S2: Coordinates of the 84 cortical points (Regions of interest, ROIs) used for connectivity analyses [21,23,24,46,47,48,49,50,51,52,53].

Author Contributions

Conceptualization, C.-C.M. and N.-S.T.; methodology, H.-A.C. and C.C.-Y.H.; software, Y.-Y.L.; validation, Y.-A.C., S.Y.P., J.J.I. and Y.-Y.L.; formal analysis, H.-A.C. and T.-C.Y.; data curation, H.-A.C. and T.-C.Y.; writing: original draft preparation, T.-C.Y.; manuscript editing, Y.-A.C., S.Y.P., and J.J.I.; visualization, Y.-Y.L.; supervision, Y.-A.C., S.Y.P. and J.J.I.; project administration, H.-A.C.; funding acquisition, H.-A.C. and C.C.-Y.H. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported in part by grants from the Medical Affairs Bureau, Ministry of National Defense (MND-MAB-D-111064 and MND-MAB-D-112065) and the National Science and Technology Council of Taiwanese Government (NSTC 111-2628-B-016-004) and the Tri-Service General Hospital (TSGH-A-111002 and TSGH-A-112011).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board of Tri-Service General Hospital, Taipei, Taiwan (protocol code: 2-106-05-123 and date of approval: 9 November 2017).

Informed Consent Statement

Written informed consent has been obtained from the patients to publish this paper.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We appreciate all the patients who have participated in this study.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

AGangular gyrus
amPFCanterior medial prefrontal cortex
ASRartifact subspace reconstruction
BAsBrodmann areas
BOLDblood oxygen level dependency
DLPFCdorsolateral prefrontal cortex
dmPFCdorsal medial prefrontal cortex
DMNdefault mode network
DSM-5Diagnostic and Statistical Manual of Mental Disorders, 5th Edition
DTIdiffusion tensor imaging
eLORETAexact low-resolution brain electromagnetic tomography
ECTelectroconvulsive therapy
EEGelectroencephalography
FDRfalse discovery rate
fMRIfunctional magnetic resonance imaging
HEOGhorizontal electrooculogram
ICAindependent component analysis
LPSlagged phase synchronization
MEGmagnetoencephalography
MNIMontreal Neurological Institute
mPFCmedial prefrontal cortex
MTGmiddle temporal gyrus
PANSSPositive and Negative Syndrome Scale
PCposterior cingulate
PHGparahippocampal gyrus
PPCposterior parietal cortex
ROIsregions of interest
rsEEGresting-state EEG
SANSScale for the Assessment of Negative Symptoms
SnPMstatistical nonparametric mapping
tACStranscranial alternating current stimulation
tDCStranscranial direct current stimulation
tEStranscranial electrical stimulation
TMStranscranial magnetic stimulation
TPJtemporoparietal junction
tRNStranscranial random noise stimulation
VEOGvertical electrooculogram
vmPFCventromedial prefrontal cortex
WMworking memory

References

  1. Correll, C.U.; Schooler, N.R. Negative Symptoms in Schizophrenia: A Review and Clinical Guide for Recognition, Assessment, and Treatment. Neuropsychiatr. Dis. Treat. 2020, 16, 519–534. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Osoegawa, C.; Gomes, J.S.; Grigolon, R.B.; Brietzke, E.; Gadelha, A.; Lacerda, A.L.T.; Dias, A.M.; Cordeiro, Q.; Laranjeira, R.; de Jesus, D.; et al. Non-invasive brain stimulation for negative symptoms in schizophrenia: An updated systematic review and meta-analysis. Schizophr. Res. 2018, 197, 34–44. [Google Scholar] [CrossRef] [PubMed]
  3. Chang, C.C.; Kao, Y.C.; Chao, C.Y.; Tzeng, N.S.; Chang, H.A. Examining bi-anodal transcranial direct current stimulation (tDCS) over bilateral dorsolateral prefrontal cortex coupled with bilateral extracephalic references as a treatment for negative symptoms in non-acute schizophrenia patients: A randomized, double-blind, sham-controlled trial. Prog. Neuropsychopharmacol. Biol. Psychiatry 2020, 96, 109715. [Google Scholar] [CrossRef] [PubMed]
  4. Chang, C.C.; Lin, Y.Y.; Tzeng, N.S.; Kao, Y.C.; Chang, H.A. Adjunct high-frequency transcranial random noise stimulation over the lateral prefrontal cortex improves negative symptoms of schizophrenia: A randomized, double-blind, sham-controlled pilot study. J. Psychiatr. Res. 2021, 132, 151–160. [Google Scholar] [CrossRef] [PubMed]
  5. Kallel, L.; Mondino, M.; Brunelin, J. Effects of theta-rhythm transcranial alternating current stimulation (4.5 Hz-tACS) in patients with clozapine-resistant negative symptoms of schizophrenia: A case series. J. Neural. Transm. 2016, 123, 1213–1217. [Google Scholar] [CrossRef] [PubMed]
  6. Elyamany, O.; Leicht, G.; Herrmann, C.S.; Mulert, C. Transcranial alternating current stimulation (tACS): From basic mechanisms towards first applications in psychiatry. Eur. Arch. Psychiatry Clin. Neurosci. 2021, 271, 135–156. [Google Scholar] [CrossRef]
  7. Kasten, F.H.; Dowsett, J.; Herrmann, C.S. Sustained Aftereffect of alpha-tACS Lasts Up to 70 min after Stimulation. Front. Hum. Neurosci. 2016, 10, 245. [Google Scholar] [CrossRef] [Green Version]
  8. Wang, X.; Cheng, B.; Roberts, N.; Wang, S.; Luo, Y.; Tian, F.; Yue, S. Shared and distinct brain fMRI response during performance of working memory tasks in adult patients with schizophrenia and major depressive disorder. Hum. Brain Mapp. 2021, 42, 5458–5476. [Google Scholar] [CrossRef]
  9. Michels, L.; Bucher, K.; Luchinger, R.; Klaver, P.; Martin, E.; Jeanmonod, D.; Brandeis, D. Simultaneous EEG-fMRI during a working memory task: Modulations in low and high frequency bands. PLoS ONE 2010, 5, e10298. [Google Scholar] [CrossRef]
  10. Polania, R.; Nitsche, M.A.; Korman, C.; Batsikadze, G.; Paulus, W. The importance of timing in segregated theta phase-coupling for cognitive performance. Curr. Biol. 2012, 22, 1314–1318. [Google Scholar] [CrossRef] [Green Version]
  11. Chang, C.C.; Huang, C.C.; Chung, Y.A.; Im, J.J.; Lin, Y.Y.; Ma, C.C.; Tzeng, N.S.; Chang, H.A. Online Left-Hemispheric In-Phase Frontoparietal Theta tACS for the Treatment of Negative Symptoms of Schizophrenia. J. Pers. Med. 2021, 11, 1114. [Google Scholar] [CrossRef]
  12. Lisman, J.; Buzsaki, G. A neural coding scheme formed by the combined function of gamma and theta oscillations. Schizophr. Bull. 2008, 34, 974–980. [Google Scholar] [CrossRef] [Green Version]
  13. Perrottelli, A.; Giordano, G.M.; Brando, F.; Giuliani, L.; Mucci, A. EEG-Based Measures in At-Risk Mental State and Early Stages of Schizophrenia: A Systematic Review. Front. Psychiatry 2021, 12, 653642. [Google Scholar] [CrossRef]
  14. Liemburg, E.; Castelein, S.; Stewart, R.; van der Gaag, M.; Aleman, A.; Knegtering, H.; Genetic, R.; Outcome of Psychosis, I. Two subdomains of negative symptoms in psychotic disorders: Established and confirmed in two large cohorts. J. Psychiatr. Res. 2013, 47, 718–725. [Google Scholar] [CrossRef] [Green Version]
  15. Cao, Y.; Han, C.; Peng, X.; Su, Z.; Liu, G.; Xie, Y.; Zhang, Y.; Liu, J.; Zhang, P.; Dong, W.; et al. Correlation Between Resting Theta Power and Cognitive Performance in Patients with Schizophrenia. Front. Hum. Neurosci. 2022, 16, 853994. [Google Scholar] [CrossRef]
  16. Garakh, Z.; Zaytseva, Y.; Kapranova, A.; Fiala, O.; Horacek, J.; Shmukler, A.; Gurovich, I.Y.; Strelets, V.B. EEG correlates of a mental arithmetic task in patients with first episode schizophrenia and schizoaffective disorder. Clin. Neurophysiol. 2015, 126, 2090–2098. [Google Scholar] [CrossRef]
  17. Krukow, P.; Jonak, K.; Karakula-Juchnowicz, H.; Podkowinski, A.; Jonak, K.; Borys, M.; Harciarek, M. Disturbed functional connectivity within the left prefrontal cortex and sensorimotor areas predicts impaired cognitive speed in patients with first-episode schizophrenia. Psychiatry Res. Neuroimaging 2018, 275, 28–35. [Google Scholar] [CrossRef] [PubMed]
  18. Andreou, C.; Leicht, G.; Nolte, G.; Polomac, N.; Moritz, S.; Karow, A.; Hanganu-Opatz, I.L.; Engel, A.K.; Mulert, C. Resting-state theta-band connectivity and verbal memory in schizophrenia and in the high-risk state. Schizophr. Res. 2015, 161, 299–307. [Google Scholar] [CrossRef] [PubMed]
  19. Krukow, P.; Jonak, K.; Grochowski, C.; Plechawska-Wojcik, M.; Karakula-Juchnowicz, H. Resting-state hyperconnectivity within the default mode network impedes the ability to initiate cognitive performance in first-episode schizophrenia patients. Prog. Neuropsychopharmacol. Biol. Psychiatry 2020, 102, 109959. [Google Scholar] [CrossRef] [PubMed]
  20. Solis-Vivanco, R.; Mondragon-Maya, A.; Reyes-Madrigal, F.; de la Fuente-Sandoval, C. Impairment of novelty-related theta oscillations and P3a in never medicated first-episode psychosis patients. NPJ Schizophr. 2021, 7, 15. [Google Scholar] [CrossRef] [PubMed]
  21. Delorme, A.; Makeig, S. EEGLAB: An open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. J. Neurosci. Methods 2004, 134, 9–21. [Google Scholar] [CrossRef] [Green Version]
  22. Chang, C.Y.; Hsu, S.H.; Pion-Tonachini, L.; Jung, T.P. Evaluation of Artifact Subspace Reconstruction for Automatic Artifact Components Removal in Multi-Channel EEG Recordings. IEEE Trans. Biomed. Eng. 2020, 67, 1114–1121. [Google Scholar] [CrossRef]
  23. Pion-Tonachini, L.; Kreutz-Delgado, K.; Makeig, S. ICLabel: An automated electroencephalographic independent component classifier, dataset, and website. Neuroimage 2019, 198, 181–197. [Google Scholar] [CrossRef] [Green Version]
  24. Pascual-Marqui, R.D.; Lehmann, D.; Koukkou, M.; Kochi, K.; Anderer, P.; Saletu, B.; Tanaka, H.; Hirata, K.; John, E.R.; Prichep, L.; et al. Assessing interactions in the brain with exact low-resolution electromagnetic tomography. Philos. Trans. A Math. Phys. Eng. Sci. 2011, 369, 3768–3784. [Google Scholar] [CrossRef]
  25. Fuchs, M.; Kastner, J.; Wagner, M.; Hawes, S.; Ebersole, J.S. A standardized boundary element method volume conductor model. Clin. Neurophysiol. 2002, 113, 702–712. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Mazziotta, J.; Toga, A.; Evans, A.; Fox, P.; Lancaster, J.; Zilles, K.; Woods, R.; Paus, T.; Simpson, G.; Pike, B.; et al. A probabilistic atlas and reference system for the human brain: International Consortium for Brain Mapping (ICBM). Philos. Trans. R. Soc. Lond. B Biol. Sci. 2001, 356, 1293–1322. [Google Scholar] [CrossRef] [PubMed]
  27. Lancaster, J.L.; Woldorff, M.G.; Parsons, L.M.; Liotti, M.; Freitas, C.S.; Rainey, L.; Kochunov, P.V.; Nickerson, D.; Mikiten, S.A.; Fox, P.T. Automated Talairach atlas labels for functional brain mapping. Hum. Brain Mapp. 2000, 10, 120–131. [Google Scholar] [CrossRef] [PubMed]
  28. Olbrich, S.; Mulert, C.; Karch, S.; Trenner, M.; Leicht, G.; Pogarell, O.; Hegerl, U. EEG-vigilance and BOLD effect during simultaneous EEG/fMRI measurement. Neuroimage 2009, 45, 319–332. [Google Scholar] [CrossRef]
  29. Alekseichuk, I.; Pabel, S.C.; Antal, A.; Paulus, W. Intrahemispheric theta rhythm desynchronization impairs working memory. Restor. Neurol. Neurosci. 2017, 35, 147–158. [Google Scholar] [CrossRef]
  30. Hu, M.L.; Zong, X.F.; Mann, J.J.; Zheng, J.J.; Liao, Y.H.; Li, Z.C.; He, Y.; Chen, X.G.; Tang, J.S. A Review of the Functional and Anatomical Default Mode Network in Schizophrenia. Neurosci. Bull. 2017, 33, 73–84. [Google Scholar] [CrossRef] [Green Version]
  31. Galindo, L.; Berge, D.; Murray, G.K.; Mane, A.; Bulbena, A.; Perez, V.; Vilarroya, O. Default Mode Network Aberrant Connectivity Associated with Neurological Soft Signs in Schizophrenia Patients and Unaffected Relatives. Front. Psychiatry 2017, 8, 298. [Google Scholar] [CrossRef] [Green Version]
  32. Mingoia, G.; Wagner, G.; Langbein, K.; Maitra, R.; Smesny, S.; Dietzek, M.; Burmeister, H.P.; Reichenbach, J.R.; Schlosser, R.G.; Gaser, C.; et al. Default mode network activity in schizophrenia studied at resting state using probabilistic ICA. Schizophr. Res. 2012, 138, 143–149. [Google Scholar] [CrossRef]
  33. Kottaram, A.; Johnston, L.A.; Tian, Y.; Ganella, E.P.; Laskaris, L.; Cocchi, L.; McGorry, P.; Pantelis, C.; Kotagiri, R.; Cropley, V.; et al. Predicting individual improvement in schizophrenia symptom severity at 1-year follow-up: Comparison of connectomic, structural, and clinical predictors. Hum. Brain Mapp. 2020, 41, 3342–3357. [Google Scholar] [CrossRef] [PubMed]
  34. Andrews-Hanna, J.R.; Smallwood, J.; Spreng, R.N. The default network and self-generated thought: Component processes, dynamic control, and clinical relevance. Ann. N. Y. Acad. Sci. 2014, 1316, 29–52. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Hu, Q.; Huang, H.; Jiang, Y.; Jiao, X.; Zhou, J.; Tang, Y.; Zhang, T.; Sun, J.; Yao, D.; Luo, C.; et al. Temporoparietal Connectivity Within Default Mode Network Associates with Clinical Improvements in Schizophrenia Following Modified Electroconvulsive Therapy. Front. Psychiatry 2021, 12, 768279. [Google Scholar] [CrossRef]
  36. Fornito, A.; Zalesky, A.; Pantelis, C.; Bullmore, E.T. Schizophrenia, neuroimaging and connectomics. Neuroimage 2012, 62, 2296–2314. [Google Scholar] [CrossRef] [PubMed]
  37. Kim, J.S.; Shin, K.S.; Jung, W.H.; Kim, S.N.; Kwon, J.S.; Chung, C.K. Power spectral aspects of the default mode network in schizophrenia: An MEG study. BMC Neurosci. 2014, 15, 104. [Google Scholar] [CrossRef] [Green Version]
  38. Vertes, R.P.; Hoover, W.B.; Viana Di Prisco, G. Theta rhythm of the hippocampus: Subcortical control and functional significance. Behav. Cogn. Neurosci. Rev. 2004, 3, 173–200. [Google Scholar] [CrossRef] [Green Version]
  39. Lai, C.Q.; Ibrahim, H.; Suandi, S.A.; Abdullah, M.Z. Convolutional Neural Network for Closed-Set Identification from Resting State Electroencephalography. Mathematics 2022, 10, 3442. [Google Scholar] [CrossRef]
  40. Rabcan, J.; Levashenko, V.; Zaitseva, E.; Kvassay, M. Review of Methods for EEG Signal Classification and Development of New Fuzzy Classification-Based Approach. IEEE Access 2020, 8, 189720–189734. [Google Scholar] [CrossRef]
  41. Brady, R.O., Jr.; Gonsalvez, I.; Lee, I.; Ongur, D.; Seidman, L.J.; Schmahmann, J.D.; Eack, S.M.; Keshavan, M.S.; Pascual-Leone, A.; Halko, M.A. Cerebellar-Prefrontal Network Connectivity and Negative Symptoms in Schizophrenia. Am. J. Psychiatry 2019, 176, 512–520. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Mahmud, M.S.; Ahmed, F.; Al-Fahad, R.; Moinuddin, K.A.; Yeasin, M.; Alain, C.; Bidelman, G.M. Decoding Hearing-Related Changes in Older Adults’ Spatiotemporal Neural Processing of Speech Using Machine Learning. Front. Neurosci. 2020, 14, 748. [Google Scholar] [CrossRef]
  43. Violante, I.R.; Li, L.M.; Carmichael, D.W.; Lorenz, R.; Leech, R.; Hampshire, A.; Rothwell, J.C.; Sharp, D.J. Externally induced frontoparietal synchronization modulates network dynamics and enhances working memory performance. Elife 2017, 6, e22001. [Google Scholar] [CrossRef] [PubMed]
  44. Sambataro, F.; Blasi, G.; Fazio, L.; Caforio, G.; Taurisano, P.; Romano, R.; Di Giorgio, A.; Gelao, B.; Lo Bianco, L.; Papazacharias, A.; et al. Treatment with olanzapine is associated with modulation of the default mode network in patients with Schizophrenia. Neuropsychopharmacology 2010, 35, 904–912. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Surguladze, S.A.; Chu, E.M.; Marshall, N.; Evans, A.; Anilkumar, A.P.; Timehin, C.; McDonald, C.; Ecker, C.; Phillips, M.L.; David, A.S. Emotion processing in schizophrenia: fMRI study of patients treated with risperidone long-acting injections or conventional depot medication. J. Psychopharmacol. 2011, 25, 722–733. [Google Scholar] [CrossRef] [PubMed]
  46. Palm, U.; Reisinger, E.; Keeser, D.; Kuo, M.F.; Pogarell, O.; Leicht, G.; Mulert, C.; Nitsche, M.A.; Padberg, F. Evaluation of sham transcranial direct current stimulation for randomized, placebo-controlled clinical trials. Brain Stimul. 2013, 6, 690–695. [Google Scholar] [CrossRef]
  47. Jaeggi, S.M.; Buschkuehl, M.; Jonides, J.; Perrig, W.J. Improving fluid intelligence with training on working memory. Proc. Natl. Acad. Sci. USA 2008, 105, 6829–6833. [Google Scholar] [CrossRef] [Green Version]
  48. Li, X.; Chu, M.Y.; Lv, Q.Y.; Hu, H.X.; Li, Z.; Yi, Z.H.; Wang, J.H.; Zhang, J.Y.; Lui, S.S.Y.; Cheung, E.F.C.; et al. The remediation effects of working memory training in schizophrenia patients with prominent negative symptoms. Cogn. Neuropsychiatry 2019, 24, 434–453. [Google Scholar] [CrossRef]
  49. Schlegelmilch, F.; Markert, S.; Berkes, S.; Schellhorn, K. Online ocular artifact removal for dc-EEG-signals: Estimation of dc-level. Biomed. Tech. 2004, 2, 340–341. [Google Scholar]
  50. Pion-Tonachini, L.; Kreutz-Delgado, K.; Makeig, S. The ICLabel dataset of electroencephalographic (EEG) independent component (IC) features. Data Brief 2019, 25, 104101. [Google Scholar] [CrossRef]
  51. Delorme, A.; Mullen, T.; Kothe, C.; Akalin Acar, Z.; Bigdely-Shamlo, N.; Vankov, A.; Makeig, S. EEGLAB, SIFT, NFT, BCILAB, and ERICA: New tools for advanced EEG processing. Comput. Intell. Neurosci. 2011, 2011, 130714. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Whitton, A.E.; Deccy, S.; Ironside, M.L.; Kumar, P.; Beltzer, M.; Pizzagalli, D.A. Electroencephalography Source Functional Connectivity Reveals Abnormal High-Frequency Communication Among Large-Scale Functional Networks in Depression. Biol. Psychiatry Cogn. Neurosci. Neuroimaging 2018, 3, 50–58. [Google Scholar] [CrossRef] [PubMed]
  53. Nichols, T.E.; Holmes, A.P. Nonparametric permutation tests for functional neuroimaging: A primer with examples. Hum. Brain Mapp. 2002, 15, 1–25. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. The 2D (upper panel) and 3D (lower panel) representation of electric field simulation of left-hemispheric frontoparietal theta in-phase transcranial alternating current stimulation (theta-tACS) by HD-Explore® (Soterix Medical, New York, NY, USA), which utilizes a finite element model of brain current flow based on an MRI derived MNI-152 standard brain template. The colorbar indicates the intensity of electrical field (V/m). In the upper panel, the center of the white circle represents the stereotaxic coordinate (x − 11, y − 18, z + 40) where the field intensity is 0.28 V/m.
Figure 1. The 2D (upper panel) and 3D (lower panel) representation of electric field simulation of left-hemispheric frontoparietal theta in-phase transcranial alternating current stimulation (theta-tACS) by HD-Explore® (Soterix Medical, New York, NY, USA), which utilizes a finite element model of brain current flow based on an MRI derived MNI-152 standard brain template. The colorbar indicates the intensity of electrical field (V/m). In the upper panel, the center of the white circle represents the stereotaxic coordinate (x − 11, y − 18, z + 40) where the field intensity is 0.28 V/m.
Biomedicines 11 00630 g001
Figure 2. The difference between active versus sham group in the change from baseline for the Positive and Negative Syndrome Scale (PANSS) negative subscale score at each post-baseline visit. Error bars indicate standard errors of the means.
Figure 2. The difference between active versus sham group in the change from baseline for the Positive and Negative Syndrome Scale (PANSS) negative subscale score at each post-baseline visit. Error bars indicate standard errors of the means.
Biomedicines 11 00630 g002
Figure 3. The left panels indicate the matrix of the t value for the mean change for lagged phase synchronization (LPS) of any pair of 84 regions of interest (ROIs) from baseline to the end of stimulation (A) and from baseline to the one-week follow-up (B). The right panels indicate the maps where t values in the matrix of the left panels that are greater than t-value threshold for statistical significance return 1 and those less than the threshold return 0.
Figure 3. The left panels indicate the matrix of the t value for the mean change for lagged phase synchronization (LPS) of any pair of 84 regions of interest (ROIs) from baseline to the end of stimulation (A) and from baseline to the one-week follow-up (B). The right panels indicate the maps where t values in the matrix of the left panels that are greater than t-value threshold for statistical significance return 1 and those less than the threshold return 0.
Biomedicines 11 00630 g003
Figure 4. (A) Relative to the sham group, the theta-tACS group had a significant decrease in theta-band (5–7 Hz) lagged phase synchronization (LPS) between the posterior cingulate (PC) and the parahippocampal gyrus (PHG) only (B) from baseline to the end of stimulation (p value = 0.0009, corrected). Color coding: PC (right hemisphere), red; PHG (right hemisphere), blue. The figure was created using eLORETA and BrainNet Viewer. Regions of interest (ROIs) shown here are displayed on a 5 × 5 × 5 MNI template brain in eLORETA for analyses (5 mm resolution is used). Error bars indicated standard errors. * p < 0.05 (corrected).
Figure 4. (A) Relative to the sham group, the theta-tACS group had a significant decrease in theta-band (5–7 Hz) lagged phase synchronization (LPS) between the posterior cingulate (PC) and the parahippocampal gyrus (PHG) only (B) from baseline to the end of stimulation (p value = 0.0009, corrected). Color coding: PC (right hemisphere), red; PHG (right hemisphere), blue. The figure was created using eLORETA and BrainNet Viewer. Regions of interest (ROIs) shown here are displayed on a 5 × 5 × 5 MNI template brain in eLORETA for analyses (5 mm resolution is used). Error bars indicated standard errors. * p < 0.05 (corrected).
Biomedicines 11 00630 g004
Figure 5. Correlation analyses in patients treated with active online theta-tACS showed that a greater reduction in the LPS between the posterior cingulate and the parahippocampal gyrus from baseline to the end of stimulation was associated with a greater reduction in PANSS negative subscale score from baseline to the end of stimulation. The regression line for the linear regression slope was shown.
Figure 5. Correlation analyses in patients treated with active online theta-tACS showed that a greater reduction in the LPS between the posterior cingulate and the parahippocampal gyrus from baseline to the end of stimulation was associated with a greater reduction in PANSS negative subscale score from baseline to the end of stimulation. The regression line for the linear regression slope was shown.
Biomedicines 11 00630 g005
Table 1. Concise demographics and clinical data of the participants.
Table 1. Concise demographics and clinical data of the participants.
tACS (n = 17)Sham (n = 18)p Value
Schizophrenia/schizoaffective disorder12/515/30.44
Gender (f/m)9/88/100.62
Handedness (r/l)16/115/30.60
Age, years old42.12 ± 8.9943.17 ± 11.200.76
Years of education14.6 ± 3.212.7 ± 2.80.07
Years since diagnosis15.7 ± 10.617.3 ± 10.60.65
Olanzapine equivalent dose, mg/day a19.59 ± 11.8319.03 ± 13.460.90
PANSS total score71.82 ± 9.6474.11 ± 7.300.43
PANSS negative subscale score19.00 ± 3.8619.83 ± 3.630.52
PANSS positive subscale score15.71 ± 5.0616.28 ± 4.080.72
PANSS general subscale score37.12 ± 5.2738.99 ± 3.910.58
SANS score50.76 ± 11.1052.61 ± 10.050.61
Abbreviations: tACS, online theta transcranial alternating current stimulation; PANSS, Positive and Negative Syndrome Scale; SANS, Scale for the Assessment of Negative Symptoms. Notes: data are presented as means ± standard deviations unless otherwise stated. a The daily dose of antipsychotic medications was converted to olanzapine equivalent.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yeh, T.-C.; Huang, C.C.-Y.; Chung, Y.-A.; Park, S.Y.; Im, J.J.; Lin, Y.-Y.; Ma, C.-C.; Tzeng, N.-S.; Chang, H.-A. Online Left-Hemispheric In-Phase Frontoparietal Theta tACS Modulates Theta-Band EEG Source-Based Large-Scale Functional Network Connectivity in Patients with Schizophrenia: A Randomized, Double-Blind, Sham-Controlled Clinical Trial. Biomedicines 2023, 11, 630. https://doi.org/10.3390/biomedicines11020630

AMA Style

Yeh T-C, Huang CC-Y, Chung Y-A, Park SY, Im JJ, Lin Y-Y, Ma C-C, Tzeng N-S, Chang H-A. Online Left-Hemispheric In-Phase Frontoparietal Theta tACS Modulates Theta-Band EEG Source-Based Large-Scale Functional Network Connectivity in Patients with Schizophrenia: A Randomized, Double-Blind, Sham-Controlled Clinical Trial. Biomedicines. 2023; 11(2):630. https://doi.org/10.3390/biomedicines11020630

Chicago/Turabian Style

Yeh, Ta-Chuan, Cathy Chia-Yu Huang, Yong-An Chung, Sonya Youngju Park, Jooyeon Jamie Im, Yen-Yue Lin, Chin-Chao Ma, Nian-Sheng Tzeng, and Hsin-An Chang. 2023. "Online Left-Hemispheric In-Phase Frontoparietal Theta tACS Modulates Theta-Band EEG Source-Based Large-Scale Functional Network Connectivity in Patients with Schizophrenia: A Randomized, Double-Blind, Sham-Controlled Clinical Trial" Biomedicines 11, no. 2: 630. https://doi.org/10.3390/biomedicines11020630

APA Style

Yeh, T. -C., Huang, C. C. -Y., Chung, Y. -A., Park, S. Y., Im, J. J., Lin, Y. -Y., Ma, C. -C., Tzeng, N. -S., & Chang, H. -A. (2023). Online Left-Hemispheric In-Phase Frontoparietal Theta tACS Modulates Theta-Band EEG Source-Based Large-Scale Functional Network Connectivity in Patients with Schizophrenia: A Randomized, Double-Blind, Sham-Controlled Clinical Trial. Biomedicines, 11(2), 630. https://doi.org/10.3390/biomedicines11020630

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop