Advertisement
Research Article|Articles in Press

Prefrontal tDCS modulates autonomic responses in COVID-19 inpatients

Open AccessPublished:March 18, 2023DOI:https://doi.org/10.1016/j.brs.2023.03.001

      Highlights

      • Bihemispheric prefrontal tDCS is safe and feasible in hospitalized patients with COVID-19.
      • Bihemispheric prefrontal tDCS showed a large effect size on heart rate variability in COVID-19 inpatients.
      • Oxygen saturation may be improved through prefrontal tDCS.
      • Bihemispheric prefrontal tDCS seems promising to enhance clinical outcomes after severe infections.

      Abstract

      Background

      maladaptive changes in the autonomic nervous system (ANS) have been observed in short and long-term phases of COVID-19 infection. Identifying effective treatments to modulate autonomic imbalance could be a strategy for preventing and reducing disease severity and induced complications.

      Objective

      to investigate the efficacy, safety, and feasibility of a single session of bihemispheric prefrontal tDCS on indicators of cardiac autonomic regulation and mood of COVID-19 inpatients.

      Methods

      patients were randomized to receive a single 30-minute session of bihemispheric active tDCS over the dorsolateral prefrontal cortex (2mA; n = 20) or sham (n = 20). Changes in time [post-pre intervention] in heart rate variability (HRV), mood, heart rate, respiratory rate, and oxygen saturation were compared between groups. Additionally, clinical worsening indicators and the occurrence of falls and skin injuries were evaluated. The Brunoni Adverse Effects Questionary was employed after the intervention.

      Results

      there was a large effect size (Hedges’ g = 0.7) of intervention on HRV frequency parameters, suggesting alterations in cardiac autonomic regulation. An increment in oxygen saturation was observed in the active group but not in the sham after the intervention (P = 0.045). There were no group differences regarding mood, incidence and intensity of adverse effects, no occurrence of skin lesions, falls, or clinical worsening.

      Conclusions

      a single prefrontal tDCS session is safe and feasible to modulate indicators of cardiac autonomic regulation in acute COVID-19 inpatients. Further research comprising a thorough assessment of autonomic function and inflammatory biomarkers is required to verify its potential to manage autonomic dysfunctions, mitigate inflammatory responses and enhance clinical outcomes.

      Keywords

      Abbreviations:

      tDCS (transcranial direct current stimulation), DLPFC (dorsolateral prefrontal cortex), ANS (autonomic nervous system), HRV (heart rate variability), SDNN (standard deviation of the NN intervals), RMSSD (the root mean square of successive differences between RR intervals), VLF (very-low-frequency), LF (low-frequency), HF (high-frequency), LFnu (low-frequency in normalized units), HFnu (high-frequency in normalized units), BDI-II (Beck Depression Inventory-II), BAI (Beck Anxiety Inventory), ICU (Intensive Care Unit), sgACC (subgenual anterior cingulate cortex), ChAP (cholinergic anti-inflammatory pathway), HPAA (hypothalamic pituitary adrenal axis)

      1. Introduction

      The complications of coronavirus disease 2019 (COVID-19) have become a challenge for the global health and scientific community. Besides the respiratory system, the disease also affects the nervous system through different mechanisms, such as inflammatory, direct viral infiltration, cerebrovascular dysfunctions, and hypoxia [
      • Boldrini M.
      • Canoll P.D.
      • Klein R.S.
      How COVID-19 affects the brain.
      ,
      • Solomon T.
      Neurological infection with SARS-CoV-2 — the story so far.
      ,
      • Mao L.
      • Jin H.
      • Wang M.
      • Hu Y.
      • Chen S.
      • He Q.
      • et al.
      Neurologic manifestations of hospitalized patients with coronavirus disease 2019 in wuhan, China.
      ]. Maladaptive changes in the autonomic nervous system (ANS) may occur in the early and long stages of COVID-19 [
      • Milovanovic B.
      • Djajic V.
      • Bajic D.
      • Djokovic A.
      Assessment of autonomic nervous system dysfunction in the early phase of infection with SARS-CoV-2 virus.
      ,
      • Buoite Stella A.
      • Furlanis G.
      • Frezza N.A.
      • Valentinotti R.
      • Ajcevic M.
      • Manganotti P.
      Autonomic dysfunction in post-COVID patients with and without neurological symptoms: a prospective multidomain observational study.
      ,
      • Scala I.
      • Rizzo P.A.
      • Bellavia S.
      • Brunetti V.
      • Col F.
      • Broccolini A.
      • et al.
      Autonomic dysfunction during acute SARS-CoV-2 infection : a systematic review.
      ] and have been strongly associated with the overproduction of pro-inflammatory cytokines, the so-called “cytokine storm” [
      • Fudim M.
      • Qadri Y.J.
      • Ghadimi K.
      • MacLeod D.B.
      • Molinger J.
      • Piccini J.P.
      • et al.
      Implications for neuromodulation therapy to control inflammation and related organ dysfunction in COVID-19.
      ]. Therefore, identifying effective interventions to restore autonomic function and minimize the inflammatory process may be a strategy to prevent and reduce clinical worsening following COVID-19 and other diseases.
      Techniques targeting neuromodulation may have therapeutic potential for managing symptoms and complications associated with severe infections like COVID-19. A recent randomized control trial investigated the effects of high-definition transcranial direct current stimulation (HD-tDCS) applied over the diaphragmatic motor cortex on the recovery of critically ill COVID-19 patients requiring pulmonary rehabilitation. The authors reported that HD-tDCS was a safe and effective add-on intervention in enhancing the clinical condition, reducing hospitalization length and the number of days of mechanical ventilation in patients who received active stimulation compared to sham [
      • Andrade S.M.
      • Cecília de Araújo Silvestre M.
      • Tenório de França E.É.
      • Bezerra Sales Queiroz M.H.
      • de Jesus Santana K.
      • Lima Holmes Madruga M.L.
      • et al.
      Efficacy and safety of HD-tDCS and respiratory rehabilitation for critically ill patients with COVID-19 the HD-RECOVERY randomized clinical trial.
      ]. However, the amount of patients with COVID-19 requiring intensive care and mechanical ventilation is relatively low (about 5% at the initial phase of the pandemic) [
      • Wu Z.
      • McGoogan J.M.
      Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China.
      ]. Thus, a distinct strategy could be identifying a tDCS configuration to prevent and treat clinical worsening instead of focusing on the diaphragm function.
      In this sense, stimulating brain regions associated with the regulation of ANS may be a promising approach to treating patients with COVID-19. As this system plays an important role in the inflammatory process [
      • Marvar P.J.
      • Harrison D.G.
      Inflammation, immunity and the autonomic nervous system. Prim. Auton. Nerv. Syst.
      ], its modulation could reduce the disease's inflammatory response besides restoring autonomic function. A recent study indeed showed improvements in the levels of inflammatory markers of hospitalized patients with COVID-19 after receiving a non-invasive vagus nerve stimulation therapy [
      • Tornero C.
      • Pastor E.
      • Garzando M del M.
      • Orduña J.
      • Forner M.J.
      • Bocigas I.
      • et al.
      Non-invasive vagus nerve stimulation for COVID-19: results from a randomized controlled trial (SAVIOR I).
      ]. In tDCS and transcranial magnetic stimulation studies, the prefrontal cortex has been suggested as a target area in the modulation of the ANS [
      • Iseger T.A.
      • van Bueren N.E.R.
      • Kenemans J.L.
      • Gevirtz R.
      • Arns M.
      A frontal-vagal network theory for Major Depressive Disorder: implications for optimizing neuromodulation techniques.
      ,
      • Makovac E.
      • Thayer J.F.
      • Ottaviani C.
      A meta-analysis of non-invasive brain stimulation and autonomic functioning: implications for brain-heart pathways to cardiovascular disease.
      ], with primary evidence of increased vagal activation in distinct conditions and pathologies [
      • Casanova M.F.
      • Hensley M.K.
      • Sokhadze E.M.
      • El-Baz A.S.
      • Wang Y.
      • Li X.
      • et al.
      Effects of weekly low-frequency rTMS on autonomic measures in children with autism spectrum disorder.
      ,
      • Nikolin S.
      • Boonstra T.W.
      • Loo C.K.
      • Martin D.
      Combined effect of prefrontal transcranial direct current stimulation and a working memory task on heart rate variability.
      ,
      • Brunoni A.R.
      • Vanderhasselt M.-A.
      • Boggio P.S.
      • Fregni F.
      • Dantas E.M.
      • Mill J.G.
      • et al.
      Polarity- and valence-dependent effects of prefrontal transcranial direct current stimulation on heart rate variability and salivary cortisol.
      ,
      • Lin Y.-Y.
      • Chang C.-C.
      • Huang C.C.-Y.
      • Tzeng N.-S.
      • Kao Y.-C.
      • Chang H.-A.
      Efficacy and neurophysiological predictors of treatment response of adjunct bifrontal transcranial direct current stimulation (tDCS) in treating unipolar and bipolar depression.
      ]. Previous accounts indeed reported alterations in indicators of cardiac autonomic modulation, through heart rate variability (HRV) analysis, after a single session of prefrontal tDCS in healthy individuals [
      • Nikolin S.
      • Boonstra T.W.
      • Loo C.K.
      • Martin D.
      Combined effect of prefrontal transcranial direct current stimulation and a working memory task on heart rate variability.
      ,
      • Brunoni A.R.
      • Vanderhasselt M.-A.
      • Boggio P.S.
      • Fregni F.
      • Dantas E.M.
      • Mill J.G.
      • et al.
      Polarity- and valence-dependent effects of prefrontal transcranial direct current stimulation on heart rate variability and salivary cortisol.
      ]. These findings corroborate the potential benefit of neuromodulation in restoring autonomic dysfunctions observed in many diseases and COVID-19 survivors, besides preventing exacerbated inflammatory responses, as abovementioned. However, evidence about early tDCS efficacy and safety is still scarce in hospital settings, particularly in COVID-19.
      The present is an exploratory, randomized, sham-controlled, double-blinded study aiming to provide insights into the safety and feasibility of a single session of bihemispheric prefrontal tDCS in COVID-19 inpatients and its acute effects on indicators of cardiac autonomic modulation assessed through the HRV. Additionally, since the relationship between HRV and mood alterations has been well established in the literature [
      • Kemp A.H.
      • Quintana D.S.
      • Gray M.A.
      • Felmingham K.L.
      • Brown K.
      • Gatt J.M.
      Impact of depression and antidepressant treatment on heart rate variability: a review and meta-analysis.
      ], the patient's mood was also an outcome investigated in this study.

      2. Material and methods

      2.1 Study design, randomization, and participants

      The present study consisted of a randomized, sham-controlled, double-blinded design conducted at the semi-intensive care unit of the Hospital Quinta D’Or (Rio de Janeiro, Brazil). Given the exploratory nature of the study and the limited resources for data collection [
      • Lakens D.
      Sample size justification.
      ], a sample of 40 participants was enrolled. This study is part of a protocol previously registered with clinicaltrials.gov (NCT04808284) and executed with no major deviations. It was approved by the institutional ethics committee (No. 4.156.081) and reported according to CONSORT guidelines [
      • Schulz K.F.
      CONSORT 2010 statement: updated guidelines for reporting parallel group randomized trials.
      ]. All patients provided informed written consent before enrolment according to the Declaration of Helsinki.
      Patients were admitted consecutively and randomly assigned to receive either sham or active-tDCS. Randomization was performed through a web-based system, and group assignments were concealed from patients and investigators who conducted the evaluation protocol.
      Recruitments were based on medical records or a physician referral according to the following inclusion criteria: 1) individuals of both sex, aged between 18-80 years; 2) confirmed infection by SARS-CoV-2; 3) ability to understand and execute the proposed protocol; 4) stable vital signs (body temperature <38 °C; respiratory rate between 12 and 30rpm; and blood pressure between 90/60mmHg and 140/90mmHg; heart rate between 50 and 120 bpm; SpO2 ≥ 90%). Exclusion criteria were: 1) dyspnoea or signs of respiratory effort; 2) hemodynamic instability; 3) deep vein thrombosis, active bleeding, use of cardiac pacemaker; 4) injury, pain, or metallic implants in the skull or scalp; 5) seizure history; 6) suspected or confirmed pregnancy; 7) concomitant or previous rheumatic or neurological diseases; 8) severe psychiatric diseases (schizophrenia, bipolar disorder, intellectual disability); 9) severe musculoskeletal and integumentary disorders; 10) severe liver or kidney disease. COVID-19 severity was classified using the Ordinal Scale for Clinical Improvement proposed by a special WHO committee [
      World Health Organization
      COVID-19 therapeutic trial synopsis.
      ].

      2.2 Experimental protocol

      To investigate the safety, feasibility, and acute effects of bihemispheric prefrontal tDCS on indicators of autonomic modulation of COVID-19 inpatients, the following evaluations were performed: clinical data collection, assessment of adverse effects of tDCS, mood assessment, and HRV recording.
      The entire protocol was conducted in the patient's semi-intensive care room. One investigator was aware of group randomization for the patients and responsible for setting up the stimulator according to the protocol for sham or active-tDCS. This investigator was not involved in any of the evaluations. Two other investigators were responsible for applying the evaluations and, as for the patients, they were unaware of the group assignment.

      2.3 Stimulation protocol

      Patients were asked to rest in a sitting position on the bed or chair to feel comfortable during the entire stimulation session. A single tDCS session was applied for both groups by a DC stimulator (NeuroConn, Germany), with a pair of saline-soaked sponge electrodes (surface 35cm2) positioned and fixed to the scalp using elastic bands. For active-tDCS, electrodes were placed bilaterally over the dorsolateral prefrontal cortex (DLPFC), specifically on F3 (anode) and F4 (cathode), according to the 10–20 EEG system [
      • Nuwer M.R.
      • Comi G.
      • Emerson R.
      • Fuglsang-Frederiksen A.
      • Guerit M.
      • Hinrichs H.
      • et al.
      IFCN standards for digital recording of clinical EEG.
      ]. This montage was chosen due to evidence of its associations with the regulation of autonomic function, comprising HRV indexes [
      • Nikolin S.
      • Boonstra T.W.
      • Loo C.K.
      • Martin D.
      Combined effect of prefrontal transcranial direct current stimulation and a working memory task on heart rate variability.
      ,
      • Brunoni A.R.
      • Vanderhasselt M.-A.
      • Boggio P.S.
      • Fregni F.
      • Dantas E.M.
      • Mill J.G.
      • et al.
      Polarity- and valence-dependent effects of prefrontal transcranial direct current stimulation on heart rate variability and salivary cortisol.
      ,
      • Cerqueira J.J.
      • Almeida O.F.X.
      • Sousa N.
      The stressed prefrontal cortex.
      ]. The current intensity was ramped up for 30s until it reached 2mA, and after 30min, the intensity was ramped down for 30s [
      • Bikson M.
      • Grossman P.
      • Thomas C.
      • Zannou A.L.
      • Jiang J.
      • Adnan T.
      • et al.
      Safety of transcranial direct current stimulation: evidence based update 2016.
      ]. For sham-tDCS, the current intensity was ramped up for 30s until it reached 2mA, and after 30s, the intensity was ramped down for 15s. To ensure the placebo effect during sham-tDCS, patients remained with the electrodes attached to the scalp for 30min. Investigators and patients were unaware of when stimulation was turned off [
      • Gandiga P.C.
      • Hummel F.C.
      • Cohen L.G.
      Transcranial DC stimulation (tDCS): a tool for double-blind sham-controlled clinical studies in brain stimulation.
      ].
      Each patient used disposable electrodes, sponges, and elastic bands to reduce the risk of infections, and all instruments were disinfected with alcohol 70° after each experimental session. In addition, strict adherence to the appropriate collection, biosafety, and personal protective equipment protocols were adopted [
      • Bikson M.
      • Hanlon C.A.
      • Woods A.J.
      • Gillick B.T.
      • Charvet L.
      • Lamm C.
      • et al.
      Guidelines for TMS/tES clinical services and research through the COVID-19 pandemic.
      ]. One of the investigators monitored patients' vital signs throughout the evaluation and intervention process. The protocol would be immediately interrupted if any exclusion criteria were identified.

      2.4 Clinical data

      Clinical data were extracted from patients' medical records, including comorbidities, clinical signs and symptoms at hospital admission, and drugs used during hospitalization. The following parameters were collected before and after tDCS intervention during the rest period that heart rate variability (HRV) was recorded: heart rate, respiratory rate, and oxygen saturation measured by pulse oximetry. The presence of severe skin lesions, fall events, and acute clinical worsening, including transfer to the Intensive Care Unit (ICU), orotracheal intubation, and death, were evaluated within one hour after the intervention.

      2.5 Adverse effects of tDCS

      The adverse effects of tDCS were assessed through the Brunoni Adverse Effects Questionnaire [
      • Brunoni A.R.
      • Amadera J.
      • Berbel B.
      • Volz M.S.
      • Rizzerio B.G.
      • Fregni F.
      A systematic review on reporting and assessment of adverse effects associated with transcranial direct current stimulation.
      ]. This questionnaire contains 11 items in which patients must score the intensity of specific symptoms (headache, neck ache, scalp pain, itching, tingling, burning sensation, redness on the skin, sleepiness, difficulty concentrating, sudden change of mood) during stimulation, on a rating scale of 1–4 (1, absent; 2, mild; 3, moderate; 4, severe). In addition, patients quantified how their experienced adverse effects were related to tDCS on a Likert-type items (1, none; 2, remote; 3, possible; 4, probable; 5, definite). They were also asked about the presence of unexpected adverse effects. The parameters considered for analysis were the occurrence, intensity, and perceived relationship of adverse effects with tDCS. Patients were also questioned whether they believed the stimulation was real or fictitious.

      2.6 Mood assessment

      Depressive and anxious symptoms were assessed before and after the intervention through the Brazilian-validated versions of the Beck Depression Inventory-II (BDI-II) [
      • Gomes-Oliveira M.H.
      • Gorenstein C.
      • Neto F.L.
      • Andrade L.H.
      • Wang Y.P.
      Validação da versão Brasileira em Português do Inventário de Depressão de Beck-II numa amostra da comunidade.
      ,

      Beck AT, Beamesderfer A. Assessment of Depression: The Depression Inventory, n.d, pp. 151–169. https://doi.org/10.1159/000395074.

      ] and the Beck Anxiety Inventory (BAI) [
      • Beck A.T.
      • Epstein N.
      • Brown G.
      • Steer R.A.
      An inventory for measuring clinical anxiety: psychometric properties.
      ] respectively. Both BDI-II and BAI were composed of 21 questions; each item was scored from 0 to 3.

      2.7 Heart rate variability (HRV)

      In the present study, indicators of cardiac autonomic modulation were defined through parameters of HRV. A heart rate monitor (Polar® RS800, Polar Electro Oy Inc., Kempele, Finland) was secured with a strap around the patient's chest and positioned at the height of the xiphoid process. A conductive gel (Mercur, Rio Grande do Sul, Brasil) was placed in the electrode areas. RR intervals were recorded at a sampling frequency of 1000Hz for six minutes before and after the tDCS intervention. The patients remained at rest during and at least 15 min before data acquisition. The patient's room was also controlled regarding any external disturbances, like noise or light changes [
      • Hartikainen J.E.K.
      • Tahvanainen K.U.O.
      • Kuusela T.A.
      Short-term measurement of heart rate variability. Clin. Guid. To card. Auton. Tests.
      ]. One investigator remained outside the room to prevent anyone from entering during the acquisition. In case of external disturbances, the HRV record was restarted.
      Signal processing was undertaken by the Kubious HRV software (Kubios Oy, Kuopio, Finland) [
      • Tarvainen M.P.
      • Niskanen J.-P.
      • Lipponen J.A.
      • Ranta-aho P.O.
      • Karjalainen P.A.
      Kubios HRV – heart rate variability analysis software.
      ] after excluding the first minute of data collected in both resting periods (before and after intervention). After visual inspection of the time series of RR interval (tachogram), HRV data from four patients (three sham and one active-tDCS) were excluded due to low-quality signals. The HRV parameters computed followed the recommendations in HRV literature for short-term recordings [
      • Hartikainen J.E.K.
      • Tahvanainen K.U.O.
      • Kuusela T.A.
      Short-term measurement of heart rate variability. Clin. Guid. To card. Auton. Tests.
      ,
      • European Society of Cardiology H.
      Rate variability, of measurement, physiological interpretation and clinical use, American TF of TES of C and TN, Electrophysiology S of P and. Guidelines Heart rate variability.
      ]. Temporal parameters included the standard deviation of the NN intervals (SDNN), and the root mean square of successive differences between RR intervals (RMSSD). Power spectral analysis of RR interval variability was performed using fast Fourier transformation and quantified into the following standard frequency-domain measurements: very-low-frequency (VLF; frequency range: ≤0.04Hz), low-frequency power (LF; frequency range: 0.04–0.15Hz) and high-frequency power (HF; frequency range: 0.15–0.40Hz), recommended by “Task force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology” [
      • European Society of Cardiology H.
      Rate variability, of measurement, physiological interpretation and clinical use, American TF of TES of C and TN, Electrophysiology S of P and. Guidelines Heart rate variability.
      ]. LF and HF powers were also computed in normalized units (LFnu and HFnu), which corresponds to the relative value of each power component (LF and HF) in proportion to the total power minus the VLF component [
      • Tarvainen M.P.
      • Niskanen J.-P.
      • Lipponen J.A.
      • Ranta-aho P.O.
      • Karjalainen P.A.
      Kubios HRV – heart rate variability analysis software.
      ]. The normalization process allows a more reliable comparison of frequency components across studies, making normalized parameters less varied by methodological differences and less sensitive to variations in total power [
      • European Society of Cardiology H.
      Rate variability, of measurement, physiological interpretation and clinical use, American TF of TES of C and TN, Electrophysiology S of P and. Guidelines Heart rate variability.
      ]. The LF/HF ratio was also computed since it is usually considered to assess the sympathetic and parasympathetic modulation of the ANS and sympathovagal balance [
      • European Society of Cardiology H.
      Rate variability, of measurement, physiological interpretation and clinical use, American TF of TES of C and TN, Electrophysiology S of P and. Guidelines Heart rate variability.
      ].

      2.8 Statistical analysis

      Parametric, inferential statistics were used when data followed a normal distribution (Kolmogorov-Smirnov test; p > 0.05), and the homogeneity of variance was confirmed (Levene's test; p > 0.05). Otherwise, non-parametric tests were applied. Differences between groups (sham and active-tDCS) related to adverse effects were analyzed using the Fishers exact test. The independent sample t-test was applied to test for group differences in age, symptoms' duration, and BDI-II score. Mann-Whitney U test was used to test for between-group differences related to disease severity, hospitalization length, the intensity of each adverse effect of tDCS, and BAI score. Missing HRV data (4 out of 40 patients) were replaced through the mean imputation method. Differences between assessments [change = post-pre intervention] between groups for HRV parameters, heart rate, respiratory rate, and oxygen saturation were assessed by one-way analysis of covariance (ANCOVA) [
      • O'Connell N.S.
      • Dai L.
      • Jiang Y.
      • Speiser J.L.
      • Ward R.
      • Wei W.
      • et al.
      Methods for analysis of pre-post data in clinical research: a comparison of five common methods.
      ]. We added age as a covariate due to its known relationship with neuromodulation outcomes [
      • Müller-Dahlhaus J.F.M.
      • Orekhov Y.
      • Liu Y.
      • Ziemann U.
      Interindividual variability and age-dependency of motor cortical plasticity induced by paired associative stimulation.
      ,
      • Fathi D.
      • Ueki Y.
      • Mima T.
      • Koganemaru S.
      • Nagamine T.
      • Tawfik A.
      • et al.
      Effects of aging on the human motor cortical plasticity studied by paired associative stimulation.
      ]. Although the randomization process aims at balancing the distribution of patients' characteristics between treatment groups, considering the variability usually observed for HRV results during short-term recordings [
      • Nunan D.
      • Sandercock G.R.H.
      • Brodie D.A.
      A quantitative systematic review of normal values for short-term heart rate variability in healthy adults.
      ], we also performed an ANCOVA including HRV baseline values and age as covariates. The effect size was computed using the Hedges' g, and statistical power was performed through a post hoc power analysis [
      • Faul F.
      • Erdfelder E.
      • Lang A.-G.
      • Buchner A.
      G*Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences.
      ]. Statistical analyses were done using MATLAB (Mathworks, USA) and the software STATISTICA 7.0 (StatSoft, USA). The significance level was set at p < 0.05.

      3. Results

      We screened 254 patients, of which 144 did not meet the eligibility criteria, 18 refused to participate, 20 were included in another trial, and 32 were either discharged from the hospital or transferred to other units. The study was conducted in 2021, and no SARS-CoV-2 vaccine was administered to the recruited patients. Forty patients were randomized into sham or active-tDCS groups (Fig. 1). There were no dropouts during the experimental protocol.
      Fig. 1
      Fig. 1Flow diagram of the study giving an overview of the hospitalized patients included in the study and analyzed. The main inclusion criteria patients excluded did not meet (N = 144) were: age greater than 80 years; unstable vital signs; and negative COVID-19 diagnosis. Abbreviation: tDCS, transcranial direct current stimulation; HRV, heart rate variability.
      All patients were hospitalized in a semi-intensive care unit due to COVID-19 diagnosis, determined by reverse-transcriptase-polymerase-chain-reaction (RT-PCR) positive for SARS-CoV-2 infection from nasopharyngeal or nasal sampling. They were evaluated within 4–23 days from the onset of symptoms, and in most cases (58 out of 60 patients), COVID-19 severity was classified as mild [
      World Health Organization
      COVID-19 therapeutic trial synopsis.
      ]. Statistical analysis showed no differences between groups for sex, age, disease severity, symptoms duration, and hospitalization length (n = 40; p > 0.49 for all cases; Table 1). Clinical data for each group are shown in Supplementary Table 1 and include comorbidities, clinical signs, and symptoms at hospital admission. Additionally, details concerning the reasons for hospitalization, the stage of the infection (defined as from the onset of symptoms), and medications (including names and dosages) are described in Supplementary Table 2.
      Table 1Comparison of clinical profile and mood between sham and active-tDCS group.
      sham-tDCSactive-tDCSp-valueHedges' gF value (p-value)
      Group factorAge factor
      Sex [male], N° (%)17(85)16(80)1.00
      Fisher's exact test.
      Age [years]49.6(12.8)50.6(10.1)0.785
      Student's t-test.
      Disease severity [score]3.4(0.7)3.2(0.4)0.534
      Mann-Whitney U test.
      Symptoms' length [days]12(5.2)12.3(4)0.812
      Student's t-test.
      Hospitalization length [days]4.6(5.3)5.1(4.3)0.490
      Mann-Whitney U test.
      BDI-II, Post-Pre tDCS [score]−0.85(1.50)−1.25(1.74)0.44
      Student's t-test.
      0.241
      BAI, Post-Pre tDCS [score]−1.45(2.11)−1.52(3.56)0.58
      Mann-Whitney U test.
      0.368
      Heart rate, Post-Pre tDCS [bpm]−1.51(4.55)−1.52(3.62)0.0040.00(0.980)0.08(0.775)
      Respiratory rate, Post-Pre tDCS [breaths/min]−0.85(3.05)−0.45(3.10)0.1270.12(0.729)2.47(0.124)
      Oxygen saturation, Post-Pre tDCS [%]−0.30(1.78)0.75(1.58)0.6114.30(0.045)*2.43(0.127)
      Data are expressed as mean (standard deviation). Abbreviation: tDCS, transcranial direct current stimulation; BDI-II, Beck Depression Inventory-II; BAI, Beck Anxiety Inventory.
      *p < 0.05.
      a Fisher's exact test.
      b Student's t-test.
      c Mann-Whitney U test.
      Table 1 reports results for mood assessments, heart rate, respiratory rate, and oxygen saturation. Positive changes were observed for oxygen saturation after the neuromodulation in the active-tDCS group with respect to negative changes in the sham group (Fig. 2A). No differences were found between groups for changes in BDI-II and BAI scores, heart rate, and respiratory rate (p > 0.05 for all cases).
      Fig. 2
      Fig. 2Changes in oxygen saturation (A) and HRV normalized parameters, HFnu (B) and LFnu (C) from baseline to after intervention within sham and active-tDCS groups. Data are expressed as mean and standard error of the mean. Light and dark gray bars represent results from sham and active-tDCS groups, respectively. The asterisk indicates a statistically significant (p < 0.05) difference between groups. Abbreviations: tDCS, transcranial direct current stimulation; n.u., normalized units.
      Table 2 shows the results for both ANCOVA analyses performed for HRV parameters. When only age was considered as a covariate, HRV normalized parameters were significantly different between assessments [post-pre intervention] among groups. HFnu, a parameter suggested to be related to the activation of the parasympathetic branch of the ANS, increased after bihemispheric prefrontal stimulation, whereas it decreased in the sham group after intervention. Likewise, LFnu, a parameter suggested to be associated with both ANS branches’ activity, decreased in the active-tDCS group and increased in the sham group. When the baseline was included as an additional covariate, there was a significant baseline effect and no group and age effects for changes in time [post-pre intervention] in the parameters LFnu and HFnu. However, there was a large effect size of tDCS intervention on these HRV frequency parameters (Table 2; Fig. 2B and C). No significant effects were found for changes in SDNN, RMSSD, absolute values of LF and HF, and LH/HF ratio (p > 0.05 for all cases).
      Table 2Comparison of heart rate variability between sham and active-tDCS group.
      Outcomes, Post-Pre tDCSsham-tDCSactive-tDCSHedges' gF value (p-value)F value (p-value)
      Group factorAge factorGroup factorAge factorBaseline factor
      SDNN [ms]1.93(8.01)1.08(4.60)0.1300.11(0.741)1.66(0.206)0.23(0.632)2.52(0.122)1.61(0.213)
      RMSSD [ms]1.61(10.93)1.87(6.35)0.0290.02(0.883)1.22(0.276)0.001(0.974)1.40(0.245)0.28(0.598)
      LF [ms2]45.38(146.73)−2.94(79.82)0.4071.43(0.240)0.56(0.458)1.06(0.311)1.46(0.235)3.03(0.091)
      HF [ms2]55.92(328.95)21.12(54.84)0.1490.18(0.675)0.24(0.627)0.20(0.656)0.26(0.612)0.03(0.857)
      LFnu5.54(9.53)−1.57(9.93)0.7144.85(0.034)*0.74(0.395)1.47(0.233)2.48(0.125)18.89(0.000)*
      HFnu−5.41(9.34)1.58(9.95)0.7074.76(0.036)*0.72(0.401)1.43(0.241)2.40(0.131)18.39(0.000)*
      Ratio LF/HF2.31(9.09)−1.43(6.43)0.4702.00(0.166)0.00(0.973)2.24(0.145)0.00(0.980)0.51(0.478)
      Data are expressed as mean (standard deviation). Abbreviations: tDCS, transcranial direct current stimulation; SDNN, standard deviation of normal to normal NN intervals; RMSSD, root mean square of successive differences between RR intervals; LF, low-frequeny component; HF, high-frequency component; nu, normalized unit;
      *p < 0.05.
      Based on the effect size estimated from our data, a single tDCS session ensured 70% statistical power. Supplementary Table 3 shows group means for each outcome evaluated before (pre-tDCS) and after (pos-tDCS) the intervention, with respective coefficients of variation and confidence intervals.
      The protocol was safe. No skin lesions, fall events, transfers to the ICU, orotracheal intubations, or deaths were found for any included patients up to one hour after the intervention. There were no statistical differences between groups for both occurrence and intensity of all the symptoms evaluated using the Brunoni tDCS Adverse Effects Questionnaire (p > 0.34 for all cases; Fig. 3). Furthermore, no patient reported severe symptoms. For the sham-tDCS group, headache and sleepiness were the most reported symptoms (20% for each symptom), while for the active-tDCS group, tingling was the most frequent symptom (30%). Other symptoms were presented in a maximum of 15% (n = 3) of patients per group. In addition, even when comparing the occurrence of adverse effects among groups separately for each intensity level (i.e., absence, mild, moderate), we still did not find statistical differences between groups for any reported symptoms (p > 0.34 for all cases).
      Fig. 3
      Fig. 3Summary of patients' reported adverse effects as frequency and intensity. Abbreviation: tDCS, Direct Current Transcranial Stimulation.
      Regarding the perceived relationship of adverse effects with tDCS (Table 3), the symptoms itching and burning sensation were among those with higher median scores for all groups. Finally, most patients in the sham-tDCS group (17 out of 20) considered stimulation was real. Regarding patients in active tDCS groups, 18 responded correctly that stimulation was real, while one patient said it was fictitious. Four patients (N sham-tDCS = 3; N active-tDCS = 1) answered they did not know.
      Table 3Perceived relationship of adverse effects reported by patients with tDCS.
      Adverse Effectsham-tDCSactive-tDCS
      Headache3.5(2–5)4(4–4)
      Neck pain4(1–4)1(1–1)
      Scalp pain5(5–5)
      Itching5(5–5)4(3–5)
      Tingling4(2–4)3(2–5)
      Burning sensation5(5–5)4(3–5)
      Sleepiness2(1–4)1(1–2)
      Trouble concentrating2(2–2)1(1–1)
      Spatial disorientation2(2–2)
      Abbreviation: tDCS, transcranial direct current stimulation. Data are expressed as median (minimum-maximum).

      4. Discussion

      In the present study, we explored the efficacy, safety, and feasibility of a single session of bihemispheric tDCS over the DLPFC on indicators of cardiac autonomic regulation assessed through HRV and mood of COVID-19 inpatients. Patients were randomized into two groups, sham-tDCS and active-tDCS. Our findings showed the intervention was safe and feasible since all patients completed the proposed protocol, and there were no mood changes, clinical worsening, presence of skin lesions, or fall events associated with the stimulation session. Furthermore, there were no differences between groups regarding the occurrence and intensity of symptoms investigated with the Brunoni tDCS Adverse Effects Questionnaire. An improvement of oxygen saturation was observed in COVID-19 patients who received active stimulation. There was a large effect size of intervention on normalized spectral HRV indexes.

      4.1 Safe and feasible bihemispheric prefrontal tDCS in hospitalized COVID-19 patients

      Despite the absence of large-scale studies with tDCS in hospitalized patients and the novelty of treating patients with COVID-19, methodological issues were carefully considered to ensure a stimulation protocol as safest and tolerable as possible. The current intensity (2mA) and electrodes’ size (35cm2) were chosen in accordance with previous studies that applied tDCS protocols over the prefrontal cortex and did not report serious adverse effects [
      • Salehinejad M.A.
      • Ghanavai E.
      • Rostami R.
      • Nejati V.
      Cognitive control dysfunction in emotion dysregulation and psychopathology of major depression (MD): evidence from transcranial brain stimulation of the dorsolateral prefrontal cortex (DLPFC).
      ,
      • Alizadehgoradel J.
      • Nejati V.
      • Sadeghi Movahed F.
      • Imani S.
      • Taherifard M.
      • Mosayebi-Samani M.
      • et al.
      Repeated stimulation of the dorsolateral-prefrontal cortex improves executive dysfunctions and craving in drug addiction: a randomized, double-blind, parallel-group study.
      ]. In addition, we used sponge electrodes dampened with saline (NaCl) solution as the method of conductivity with the scalp, frequently used in tDCS studies [
      • Dundas J.E.
      • Thickbroom G.W.
      • Mastaglia F.L.
      Perception of comfort during transcranial DC stimulation: effect of NaCl solution concentration applied to sponge electrodes.
      ].
      This study had no dropouts or interruptions throughout the intervention protocol and no occurrence of serious adverse effects, as discussed as follows. The intensity and frequency of adverse effects reported by the patients after the tDCS protocol were not group-dependent. Overall, the symptoms presented were classified as mild intensity and occurred at a similar rate between active and sham groups. No patient reported severe symptoms. Headache was mainly reported by patients in the sham group and may have been a symptom of the disease itself [
      • Mao L.
      • Jin H.
      • Wang M.
      • Hu Y.
      • Chen S.
      • He Q.
      • et al.
      Neurologic manifestations of hospitalized patients with coronavirus disease 2019 in wuhan, China.
      ]. The symptom of tingling, reported primarily in the active stimulation group, is expected during transcutaneous electrical stimulations, likely due to the activation of specific cutaneous receptors [
      • Garnsworthy R.K.
      • Gully R.L.
      • Kenins P.
      • Westerman R.A.
      Transcutaneous electrical stimulation and the sensation of prickle.
      ]. Our findings agree with evidence reporting headache and tingling as common mild symptoms during tDCS protocols applied in healthy and pathological individuals [
      • Brunoni A.R.
      • Amadera J.
      • Berbel B.
      • Volz M.S.
      • Rizzerio B.G.
      • Fregni F.
      A systematic review on reporting and assessment of adverse effects associated with transcranial direct current stimulation.
      ,
      • Poreisz C.
      • Boros K.
      • Antal A.
      • Paulus W.
      Safety aspects of transcranial direct current stimulation concerning healthy subjects and patients.
      ]. Other adverse effects were observed in the abovementioned studies, e.g., pain, fatigue, nervousness, and discomfort. These symptoms, however, were not presented in our sample.
      In addition to the adverse effects, aspects that could contraindicate the tDCS application, such as worsening health conditions, were considered. Those outcome measures are not commonly reported in tDCS safety studies, likely because most of the available protocols were performed on outpatients. Even if the present study included critically ill patients from a semi-intensive care unit, there were no episodes of falls, and no patient was transferred to ICU, intubated, or died within one hour after the intervention applied in our research.

      4.2 Prefrontal tDCS induced HRV and oxygen saturation changes in COVID-19 patients

      HRV results indicated a single bihemispheric prefrontal tDCS session may be a potential approach to modulate the cardiac autonomic responses of hospitalized patients with COVID-19. The literature well established that autonomic regulation is a key regulator of the HRV power spectrum. The high-frequency (HF) component of HRV is strongly associated with cardiovagal activity [
      • Pagani M.
      • Lombardi F.
      • Guzzetti S.
      • Rimoldi O.
      • Furlan R.
      • Pizzinelli P.
      • et al.
      Power spectral analysis of heart rate and arterial pressure variabilities as a marker of sympatho-vagal interaction in man and conscious dog.
      ,
      • Malliani A.
      • Pagani M.
      • Lombardi F.
      • Cerutti S.
      Cardiovascular neural regulation explored in the frequency domain.
      ,
      • Pomeranz B.
      • Macaulay R.J.
      • Caudill M.A.
      • Kutz I.
      • Adam D.
      • Gordon D.
      • et al.
      Assessment of autonomic function in humans by heart rate spectral analysis.
      ,
      • Akselrod S.
      • Gordon D.
      • Ubel F.
      • Shannon D.C.
      • Barger A.C.
      • Cohen R.J.
      Power spectrum analysis of heart rate fluctuation: a quantitative probe of beat-to-beat cardiovascular control.
      ], while the interpretation of low-frequency (LF) is still questionable due to its relation with both ANS branches [
      • Malliani A.
      • Pagani M.
      • Lombardi F.
      • Cerutti S.
      Cardiovascular neural regulation explored in the frequency domain.
      ,
      • Pomeranz B.
      • Macaulay R.J.
      • Caudill M.A.
      • Kutz I.
      • Adam D.
      • Gordon D.
      • et al.
      Assessment of autonomic function in humans by heart rate spectral analysis.
      ,
      • Akselrod S.
      • Gordon D.
      • Ubel F.
      • Shannon D.C.
      • Barger A.C.
      • Cohen R.J.
      Power spectrum analysis of heart rate fluctuation: a quantitative probe of beat-to-beat cardiovascular control.
      ] or sympathetic activity [
      • Pagani M.
      • Lombardi F.
      • Guzzetti S.
      • Rimoldi O.
      • Furlan R.
      • Pizzinelli P.
      • et al.
      Power spectral analysis of heart rate and arterial pressure variabilities as a marker of sympatho-vagal interaction in man and conscious dog.
      ]. The ultra-low and very-high frequency components were not considered for analysis because they are not recommended for short-term recordings [
      • Hartikainen J.E.K.
      • Tahvanainen K.U.O.
      • Kuusela T.A.
      Short-term measurement of heart rate variability. Clin. Guid. To card. Auton. Tests.
      ,
      Electrophysiology TF of the ES of C the NA
      Heart rate variability.
      ]. When we performed the ANCOVA analysis considering age as a covariate, our results showed an increase in the HFnu parameter and a decrease in the LFnu parameter after the intervention in the active-tDCS group, compared to the sham group. Nevertheless, the inclusion of baseline as a covariate in addition to age showed it had a significant effect on the changes in the parameters LFnu and HFnu. These findings indicate that large baseline differences for these parameters contributed to the group differences found when ANCOVA analysis considered only age as a covariate. Despite the lack of a statistically significant group effect, we would like to point out the large effect size of tDCS intervention on both HRV frequency parameters (Hedges’ g: LFnu = 0.714 and HFnu = 0.707 [
      • Cohen J.
      Statistical power analysis for the behavioral sciences.
      ]). The other HRV parameters evaluated did not show statistical changes between groups or large effect sizes. By observing their coefficient of variation (CV) values in Supplementary Table 3, these data show a larger variability when compared to LFnu and HFnu, decreasing the probability of finding different responses between conditions. Moreover, the normalization process is less sensitive to variations in total power, allowing a more reliable comparison of frequency components [
      • Hartikainen J.E.K.
      • Tahvanainen K.U.O.
      • Kuusela T.A.
      Short-term measurement of heart rate variability. Clin. Guid. To card. Auton. Tests.
      ,
      • European Society of Cardiology H.
      Rate variability, of measurement, physiological interpretation and clinical use, American TF of TES of C and TN, Electrophysiology S of P and. Guidelines Heart rate variability.
      ]. Such variability in HRV parameters was indeed reported elsewhere [
      • Nunan D.
      • Sandercock G.R.H.
      • Brodie D.A.
      A quantitative systematic review of normal values for short-term heart rate variability in healthy adults.
      ]. Therefore, as observed in our findings, the consideration of baseline in the statistical model was fundamental to show the need for planning studies with larger samples to explore the effects of bihemispheric prefrontal tDCS on indicators of cardiac autonomic responses. The observed effects here can be useful for such planning.
      Overall, the large effect size found for the normalized spectral parameters indicated that the bihemispheric dorsolateral prefrontal cortex (DLPFC) may be a promising tool to modulate the cardiac autonomic response of patients. The magnitude of the effect and the validity of the results may be more important than the statistical significance since they emphasize the clinical relevance of the treatment [
      • Wasserstein R.L.
      • Schirm A.L.
      • Lazar N.A.
      Moving to a world beyond “ p < 0.05.”.
      ]. These results corroborate previous accounts that observed alterations in HRV parameters following a single session of tDCS over the DLPFC of healthy individuals [
      • Nikolin S.
      • Boonstra T.W.
      • Loo C.K.
      • Martin D.
      Combined effect of prefrontal transcranial direct current stimulation and a working memory task on heart rate variability.
      ,
      • Brunoni A.R.
      • Vanderhasselt M.-A.
      • Boggio P.S.
      • Fregni F.
      • Dantas E.M.
      • Mill J.G.
      • et al.
      Polarity- and valence-dependent effects of prefrontal transcranial direct current stimulation on heart rate variability and salivary cortisol.
      ]. There is evidence that the HFnu and LFnu changes found in the present study suggest an increase in patients’ vagal activity and attenuation of sympathetic activity [
      • Pagani M.
      • Lombardi F.
      • Guzzetti S.
      • Rimoldi O.
      • Furlan R.
      • Pizzinelli P.
      • et al.
      Power spectral analysis of heart rate and arterial pressure variabilities as a marker of sympatho-vagal interaction in man and conscious dog.
      ,
      • Malliani A.
      • Pagani M.
      • Lombardi F.
      • Cerutti S.
      Cardiovascular neural regulation explored in the frequency domain.
      ,
      • Pomeranz B.
      • Macaulay R.J.
      • Caudill M.A.
      • Kutz I.
      • Adam D.
      • Gordon D.
      • et al.
      Assessment of autonomic function in humans by heart rate spectral analysis.
      ,
      • Akselrod S.
      • Gordon D.
      • Ubel F.
      • Shannon D.C.
      • Barger A.C.
      • Cohen R.J.
      Power spectrum analysis of heart rate fluctuation: a quantitative probe of beat-to-beat cardiovascular control.
      ]. A thorough assessment of autonomic function [
      • Cheshire W.P.
      • Freeman R.
      • Gibbons C.H.
      • Cortelli P.
      • Wenning G.K.
      • Hilz M.J.
      • et al.
      Electrodiagnostic assessment of the autonomic nervous system: a consensus statement endorsed by the American autonomic society, American academy of neurology, and the international federation of clinical neurophysiology.
      ] would be suitable to further confirm such modulations in SNA. Lastly, BDI-II and BAI scores did not change after the intervention for both groups, suggesting HRV indexes results were not associated with mood changes [
      • Kemp A.H.
      • Quintana D.S.
      • Gray M.A.
      • Felmingham K.L.
      • Brown K.
      • Gatt J.M.
      Impact of depression and antidepressant treatment on heart rate variability: a review and meta-analysis.
      ].
      The electrical field generated by the active-tDCS may have caused a direct stimulation of cortical neurons, modulating the excitability of brain regions associated with autonomic cardiovascular control. These findings corroborate a previous meta-analysis [
      • Makovac E.
      • Thayer J.F.
      • Ottaviani C.
      A meta-analysis of non-invasive brain stimulation and autonomic functioning: implications for brain-heart pathways to cardiovascular disease.
      ] that identified the DLPFC as an effective brain region in regulating HRV, likely due to functional connectivity between the DLPFC and deeper brain regions associated with autonomic cardiovascular control, such as the subgenual anterior cingulate cortex (sgACC) [
      • Fox M.D.
      • Buckner R.L.
      • White M.P.
      • Greicius M.D.
      • Pascual-leone A.
      Efficacy of transcranial magnetic stimulation targets for depression is related to intrinsic functional connectivity with the subgenual cingulate.
      ,
      • Iseger T.A.
      • Padberg F.
      • Kenemans J.L.
      • Gevirtz R.
      • Arns M.
      Neuro-Cardiac-Guided TMS (NCG-TMS): probing DLPFC-sgACC-vagus nerve connectivity using heart rate e First results.
      ]. On the other hand, stimulation of peripheral nerves, as suggested by the patients' reports of tingling sensation [
      • Kessler S.K.
      • Turkeltaub P.E.
      • Benson J.G.
      • Hamilton R.H.
      Differences in the experience of active and sham transcranial direct current stimulation.
      ], may also have been a mechanism underpinning our findings. A tingling sensation suggests peripheral nerves, such as the trigeminal and occipital nerves, were likely stimulated in the patients evaluated [
      • van Boekholdt L.
      • Kerstens S.
      • Khatoun A.
      • Asamoah B.
      • Mc Laughlin M.
      tDCS peripheral nerve stimulation: a neglected mode of action?.
      ]. The brain circuits associated with sensory information processing are still unknown in the tDCS field. The tDCS protocol we applied in our study does not allow us to discriminate the contribution of transcranial and/or transcutaneous mechanisms. Boekholdt et al. (2020) [
      • van Boekholdt L.
      • Kerstens S.
      • Khatoun A.
      • Asamoah B.
      • Mc Laughlin M.
      tDCS peripheral nerve stimulation: a neglected mode of action?.
      ] suggest including specific control groups in future tDCS protocols, in which transcranial and transcutaneous mechanisms can be evaluated independently. The approaches proposed by the authors to investigate each mechanism are, respectively, blocking the peripherical stimulation through topical anesthetic creams and targeting peripheral sites more distant from the brain. However, more studies are required in this field.
      Studies have also reported symptoms of autonomic dysfunction in the early and long stages of COVID-19 infection, such as orthostatic intolerance [
      • Milovanovic B.
      • Djajic V.
      • Bajic D.
      • Djokovic A.
      Assessment of autonomic nervous system dysfunction in the early phase of infection with SARS-CoV-2 virus.
      ,
      • Buoite Stella A.
      • Furlanis G.
      • Frezza N.A.
      • Valentinotti R.
      • Ajcevic M.
      • Manganotti P.
      Autonomic dysfunction in post-COVID patients with and without neurological symptoms: a prospective multidomain observational study.
      ,
      • Scala I.
      • Rizzo P.A.
      • Bellavia S.
      • Brunetti V.
      • Col F.
      • Broccolini A.
      • et al.
      Autonomic dysfunction during acute SARS-CoV-2 infection : a systematic review.
      ]. Sympathetic overactivity was shown to be a mark of autonomic dysfunction in severe COVID-19 survivors [
      • Faria D.
      • Moll-Bernardes R.J.
      • Testa L.
      • Moniz C.M.V.
      • Rodrigues E.C.
      • Rodrigues A.G.
      • et al.
      Sympathetic neural overdrive, aortic stiffening, endothelial dysfunction, and impaired exercise capacity in severe COVID-19 survivors: a mid-term study of cardiovascular sequelae.
      ]. The possible mechanisms underlying the sympathetic overactivity in COVID-19 may be the overproduction of pro-inflammatory cytokines, alterations in blood gases (e.g., hypoxia), emotional distress, and dysregulation of the renin–angiotensin–aldosterone system [
      • Fudim M.
      • Qadri Y.J.
      • Ghadimi K.
      • MacLeod D.B.
      • Molinger J.
      • Piccini J.P.
      • et al.
      Implications for neuromodulation therapy to control inflammation and related organ dysfunction in COVID-19.
      ,
      • Porzionato A.
      • Emmi A.
      • Barbon S.
      • Boscolo‐Berto R.
      • Stecco C.
      • Stocco E.
      • et al.
      Sympathetic activation: a potential link between comorbidities and COVID‐19.
      ]. Additionally, Junarta et al. [
      • Junarta J.
      • Riley J.M.
      • Pavri B.B.
      Describing heart rate variability in patients with chronic atrial fibrillation during hospitalization for COVID‐19.
      ] reported reduced HRV in patients with diagnosis of chronic atrial fibrillation during COVID-19 acute phase, with respect to the period before the pandemic, suggesting a reduced cardiac parasympathetic tone. Therefore, the HRV results observed in our study prompt further investigation to verify whether bihemispheric prefrontal tDCS could be a potential therapeutic tool for restoring sympathovagal imbalance caused by this disease, preventing the development of dysautonomia and subsequent involvement of several organs [
      • Al-kuraishy H.M.
      • Al-Gareeb A.I.
      • Qusti S.
      • Alshammari E.M.
      • Gyebi G.A.
      • Batiha G.E.-S.
      Covid-19-Induced dysautonomia: a menace of sympathetic storm.
      ].
      Depression and anxiety have been associated with autonomic dysfunction, mainly characterized by decreased vagal activity [
      • Thayer J.F.
      • Siegle G.J.
      Neurovisceral integration in cardiac and emotional regulation.
      ]. They are among the most common symptoms reported in mid-term and long-term phases of COVID-19 (the second known as long COVID), reaching up to 40% of survivors and associated with risk of mortality [
      • Vai B.
      • Mazza M.G.
      • Delli Colli C.
      • Foiselle M.
      • Allen B.
      • Benedetti F.
      • et al.
      Mental disorders and risk of COVID-19-related mortality, hospitalisation, and intensive care unit admission: a systematic review and meta-analysis.
      ,
      • Premraj L.
      • Kannapadi N.V.
      • Briggs J.
      • Seal S.M.
      • Battaglini D.
      • Fanning J.
      • et al.
      Mid and long-term neurological and neuropsychiatric manifestations of post-COVID-19 syndrome: a meta-analysis.
      ,
      • Mazza M.G.
      • De Lorenzo R.
      • Conte C.
      • Poletti S.
      • Vai B.
      • Bollettini I.
      • et al.
      Anxiety and depression in COVID-19 survivors: role of inflammatory and clinical predictors.
      ]. Consequently, researchers have highlighted the necessity of developing intervention strategies to reduce the burden of such disorders due to COVID-19. Therefore, bihemispheric prefrontal tDCS can be a promising adjunctive therapy to modulate the SNA of COVID-19 survivors and help to treat neuropsychiatric manifestations in this population.
      ANS modulation could also be a strategy to regulate exacerbated inflammatory responses to severe infections, such as COVID-19. The ANS is responsible for regulating inflammatory and immune responses throughout the modulation of vagus nerve activity, the main component of the parasympathetic nervous system [
      • Pavlov V.A.
      • Tracey K.J.
      The vagus nerve and the inflammatory reflex—linking immunity and metabolism.
      ]. The activation of vagal efferent fibers induces two distinct mechanisms of the immune response: the cholinergic anti-inflammatory pathways (ChAP) and the hypothalamic pituitary adrenal axis (HPAA). The ChAP suppresses pro-inflammatory cytokines through acetylcholine release induced by activation of T-cells in the spleen. The HPAA releases anti-inflammatory glucocorticoids via modulation of cells in adrenal glands [
      • Bonaz B.
      • Picq C.
      • Sinniger V.
      • Mayol J.F.
      • Clarençon D.
      Vagus nerve stimulation: from epilepsy to the cholinergic anti-inflammatory pathway.
      ,
      • Bonaz B.
      • Sinniger V.
      • Pellissier S.
      Vagus nerve stimulation: a new promising therapeutic tool in inflammatory bowel disease.
      ,
      • Borovikova L.V.
      • Ivanova S.
      • Zhang M.
      • Yang H.
      • Botchkina G.I.
      • Watkins L.R.
      • et al.
      Vagus nerve stimulation attenuates the systemic inflammatory response to endotoxin.
      ,
      • Tracey K.J.
      Physiology and immunology of the cholinergic antiinflammatory pathway.
      ]. Both mechanisms can be activated through vagus nerve stimulation, as previous studies reported reduced inflammation in distinct pathological conditions [
      • Stavrakis S.
      • Humphrey M.B.
      • Scherlag B.
      • Iftikhar O.
      • Parwani P.
      • Abbas M.
      • et al.
      Low-level vagus nerve stimulation suppresses post-operative atrial fibrillation and inflammation.
      ,
      • Mondal B.
      • Choudhury S.
      • Simon B.
      • Baker M.R.
      • Kumar H.
      Noninvasive vagus nerve stimulation improves gait and reduces freezing of gait in Parkinson's disease.
      ,
      • Bremner J.D.
      • Gurel N.Z.
      • Jiao Y.
      • Wittbrodt M.T.
      • Levantsevych O.M.
      • Huang M.
      • et al.
      Transcutaneous vagal nerve stimulation blocks stress-induced activation of Interleukin-6 and interferon-γ in posttraumatic stress disorder: a double-blind, randomized, sham-controlled trial.
      ]. In COVID-19, encouraging results have been reported with decreases in levels of inflammatory markers following therapies with non-invasive vagus nerve stimulation [
      • Tornero C.
      • Pastor E.
      • Garzando M del M.
      • Orduña J.
      • Forner M.J.
      • Bocigas I.
      • et al.
      Non-invasive vagus nerve stimulation for COVID-19: results from a randomized controlled trial (SAVIOR I).
      ] and transcutaneous auricular vagus nerve stimulation [
      • Boezaart A.P.
      • Botha D.A.
      Treatment of stage 3 COVID-19 with transcutaneous auricular vagus nerve stimulation drastically reduces interleukin-6 blood levels: a report on two cases.
      ]. If confirmed the hypothesis that active bihemispheric prefrontal tDCS can heighten vagal activation of COVID-19 patients, it would be a potential intervention to mitigate hyperinflammatory syndrome associated with COVID-19 and other diseases.
      Our findings also showed an increase in oxygen saturation of patients in the active tDCS with respect to the sham group after the intervention, with a medium estimated effect size (Hedges’ g = 0.61 [
      • Cohen J.
      Statistical power analysis for the behavioral sciences.
      ]). These results seem very interesting since oxygen saturation has been widely compared among COVID-19 survivors [
      • Sobel J.A.
      • Levy J.
      • Almog R.
      • Reiner-Benaim A.
      • Miller A.
      • Eytan D.
      • et al.
      Descriptive characteristics of continuous oximetry measurement in moderate to severe covid-19 patients.
      ,
      • Spaccarotella C.
      • Polimeni A.
      • Mancuso C.
      • Pelaia G.
      • Esposito G.
      • Indolfi C.
      Assessment of non-invasive measurements of oxygen saturation and heart rate with an apple smartwatch: comparison with a standard pulse oximeter.
      ,
      • Motta L.P.
      • da Silva P.P.F.
      • Borguezan B.M.
      • do Amaral JLM.
      • Milagres L.G.
      • Bóia M.N.
      • et al.
      An emergency system for monitoring pulse oximetry, peak expiratory flow, and body temperature of patients with COVID-19 at home: development and preliminary application.
      ]. Although we do not have data to support this assumption, oxygen saturation and autonomic changes may be related. The sympathetic discharge has an important role in the development of pulmonary dysfunction. It may affect the pulmonary vascular bed through adrenergic receptors, leading to local inflammatory response due to activation of macrophages by catecholamines and subsequent release of pro-inflammatory cytokines [
      • Flierl M.A.
      • Rittirsch D.
      • Nadeau B.A.
      • Sarma J.V.
      • Day D.E.
      • Lentsch A.B.
      • et al.
      Upregulation of phagocyte-derived catecholamines augments the acute inflammatory response.
      ]. In addition, catecholamines release provokes increased pulmonary vascular resistance [
      • Hyman A.L.
      • Kadowitz P.J.
      Evidence for existence of postjunctional alpha 1- and alpha 2-adrenoceptors in cat pulmonary vascular bed.
      ,
      • Hyman A.L.
      • Kadowitz P.J.
      Enhancement of alpha- and beta-adrenoceptor responses by elevations in vascular tone in pulmonary circulation.
      ], reducing the distribution of blood to peripheral alveoli and hindering a proper gas exchange. On the other hand, parasympathetic activity elicits pulmonary vasodilator responses under vasoconstriction tone [
      • Nandiwada P.A.
      • Hyman A.L.
      • Kadowitz P.J.
      Pulmonary vasodilator responses to vagal stimulation and acetylcholine in the cat.
      ], besides induces anti-inflammatory effects as abovementioned. Therefore, future studies should investigate whether bihemispheric prefrontal tDCS can indeed attenuate sympathetic activation and heightened parasympathetic activity, and whether these changes can contribute to reducing pulmonary vascular resistance and local inflammation, enhancing oxygen saturation. If this assumption is confirmed, bihemispheric prefrontal tDCS may be a helpful add-on therapy to re-establish sufficient oxygen supply, usually required in about 70% of COVID-19 inpatients [
      • Yang X.
      • Yu Y.
      • Xu J.
      • Shu H.
      • Xia J.
      • Liu H.
      • et al.
      Clinical course and outcomes of critically ill patients with SARS-CoV-2 pneumonia in Wuhan, China: a single-centered, retrospective, observational study.
      ,
      • Myers L.C.
      • Parodi S.M.
      • Escobar G.J.
      • Liu V.X.
      Characteristics of hospitalized adults with COVID-19 in an integrated health care system in California.
      ].
      The hypothesis discussed above corroborates studies that observed acute improvements (i.e., until 90 min after intervention) in respiratory symptoms (e.g., dyspnea) and levels of inflammatory markers after applying non-invasive vagus nerve stimulation therapy in patients with acute bronchoconstriction and primary Sjögren's syndrome, respectively [
      • Tarn J.
      • Legg S.
      • Mitchell S.
      • Simon B.
      • Ng W.-F.
      The effects of noninvasive vagus nerve stimulation on fatigue and immune responses in patients with primary Sjögren’s syndrome.
      ,
      • Staats P.
      • Emala C.
      • Simon B.
      • Errico J.P.
      Neurostimulation for asthma. Neuromodulation.
      ]. The possibility of suppressing pro-inflammatory cytokines, restoring autonomic balance, and enhancing oxygen saturation through bihemispheric prefrontal tDCS therapy, may hinder clinical worsening during severe infections like COVID-19, which can result in acute respiratory distress syndrome (ARDS), multiple organ failure, and death [
      • Wu Z.
      • McGoogan J.M.
      Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China.
      ]. Further research investigating the safety and efficacy of more tDCS sessions on modulation and inflammatory markers would be helpful to support the clinical benefits of such intervention.

      4.3 Study limitations

      The present study was a first approach to investigate the safety and feasibility of a single bihemispheric prefrontal tDCS session in hospitalized patients with COVID-19 and its acute effects on indicators of cardiac autonomic modulation and mood. A priori power analysis (i.e., sample size calculation) was not performed because neuromodulation studies in a hospital setting are rare, and to our knowledge, no previous work investigated bilateral prefrontal stimulation. However, based on the effect size estimated from our data, a single tDCS session ensured a relatively high (70%) statistical power for the outcomes evaluated.
      Since this study was conducted during the pandemic period, with a high risk of SARS-CoV-2 virus transmission, a non-COVID-19 control group was not evaluated. Thus, it is not possible to infer whether the acute findings observed after the intervention were specific to COVID-19 disease. However, bihemispheric prefrontal tDCS showed to be safe and feasible in this population.
      In addition to HRV, other tests should be used to thoroughly assess the efficacy of bihemispheric prefrontal tDCS on the autonomic function of COVID-19 patients, including evaluation of cardiovascular sympathetic adrenergic, cardiovagal, and sudomotor (sweating) function [
      • Cheshire W.P.
      • Freeman R.
      • Gibbons C.H.
      • Cortelli P.
      • Wenning G.K.
      • Hilz M.J.
      • et al.
      Electrodiagnostic assessment of the autonomic nervous system: a consensus statement endorsed by the American autonomic society, American academy of neurology, and the international federation of clinical neurophysiology.
      ]. Considering the clinical conditions of COVID-19 inpatients and the limitations of performing an experimental protocol in a hospital setting, a short-term HRV recording was suitable to investigate the safety and acute effects of bihemispheric prefrontal tDCS on indicators of autonomic modulation. We computed and interpreted the reliable HRV parameters for short-term recordings, as recommended in HRV literature [
      • Hartikainen J.E.K.
      • Tahvanainen K.U.O.
      • Kuusela T.A.
      Short-term measurement of heart rate variability. Clin. Guid. To card. Auton. Tests.
      ,
      • European Society of Cardiology H.
      Rate variability, of measurement, physiological interpretation and clinical use, American TF of TES of C and TN, Electrophysiology S of P and. Guidelines Heart rate variability.
      ]. However, for HRV analysis, long-term recordings are considered the reference standard for clinical evaluation. It has a higher predictive power and is expected to provide a more accurate reflection of the overall responsiveness of the cardiovascular system than short-term recordings [
      • Kleiger R.E.
      • Stein P.K.
      • Bigger J.T.
      Heart rate variability: measurement and clinical utility.
      ].
      Finally, the use of concomitant medications by the patients evaluated in this study can be a confounding factor. Although questionable, some medications like beta blockers and ACE inhibitors may affect HRV [
      • Silva F.C.
      • Paiva F.A.
      • Müller-Ribeiro F.C.
      • Caldeira H.M.A.
      • Fontes M.A.P.
      • de Menezes R.C.A.
      • et al.
      Chronic treatment with ivabradine does not affect cardiovascular autonomic control in rats.
      ,
      • Dias da Silva V.J.
      • Tobaldini E.
      • Rocchetti M.
      • Wu M.A.
      • Malfatto G.
      • Montano N.
      • et al.
      Modulation of sympathetic activity and heart rate variability by ivabradine.
      ].
      To understand the clinical benefits of bihemispheric prefrontal tDCS in COVID-19 survivors, further research with greater samples is needed to investigate the efficacy and long-term effects of a therapy with more sessions, including a control group, a complete assessment of autonomic function, and analyses of inflammatory markers.

      5. Conclusions

      Our findings are the first insights into the safety and feasibility of a single bihemispheric prefrontal tDCS session in COVID-19 inpatients and its acute effects on cardiac autonomic modulation and mood. Bihemispheric prefrontal tDCS was safe and feasible since no mood impairments, clinical worsening or fall events were observed in the patients after the intervention, and adverse effects were neither severe nor tDCS-dependent. An improvement in oxygen saturation was also observed after the stimulation. A large effect size was observed for HRV frequency indexes suggesting the intervention may have attenuated sympathetic and increased vagal activities.
      Our findings also reinforced the importance of considering baseline measures when assessing intervention effects on studies outcomes, even in exploratory research with relatively small sample sizes. Further randomized clinical trials with greater sample sizes and more tDCS sessions, applying a thorough assessment of autonomic function and inflammatory biomarkers, can foster the investigation of bihemispheric prefrontal tDCS as a potential adjuvant therapeutic tool to help manage autonomic dysfunctions, mitigate inflammatory responses, and enhance clinical outcomes.

      Funding

      This work was supported by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq); the Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ); the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior(CAPES); and the Instituto D’Or de Pesquisa e Ensino (IDOR). TPP received stipends from CNPq and FAPERJ. EAF received stipends from CAPES (Finance Code 001). The funding sources had no involvement in study design, in the collection, analysis and interpretation of data, in the writing of the report, and in the decision to submit the article for publication.

      CRediT authorship contribution statement

      Talita P. Pinto: Conceptualization, Methodology, Software, Validation, Formal analysis, Investigation, Data curation, Writing – original draft. Jacqueline Cunha Inácio: Methodology, Investigation, Writing – original draft. Erivelton de Aguiar Ferreira: Methodology, Investigation, Writing – original draft. Arthur de Sá Ferreira: Formal analysis, Writing – review & editing. Felipe Kenji Sudo: Methodology, Writing – review & editing, Funding acquisition. Fernanda Tovar-Moll: Conceptualization, Resources, Writing – review & editing, Supervision, Funding acquisition. Erika Rodrigues: Conceptualization, Methodology, Resources, Writing – original draft, Supervision, Project administration, Funding acquisition.

      Acknowledgments

      We thank Dr. Thiago Lemos and Dr. Myriam Monteiro for their support during the research project development. We also thank Dr. Patricia Vigário and MSc Roberto Costa for their support in methodological procedures.

      Appendix A. Supplementary data

      The following is/are the supplementary data to this article.

      References

        • Boldrini M.
        • Canoll P.D.
        • Klein R.S.
        How COVID-19 affects the brain.
        JAMA Psychiatr. 2021; 78: 682https://doi.org/10.1001/jamapsychiatry.2021.0500
        • Solomon T.
        Neurological infection with SARS-CoV-2 — the story so far.
        Nat Rev Neurol. 2021; 17: 65-66https://doi.org/10.1038/s41582-020-00453-w
        • Mao L.
        • Jin H.
        • Wang M.
        • Hu Y.
        • Chen S.
        • He Q.
        • et al.
        Neurologic manifestations of hospitalized patients with coronavirus disease 2019 in wuhan, China.
        JAMA Neurol. 2020; 77: 683-690https://doi.org/10.1001/jamaneurol.2020.1127
        • Milovanovic B.
        • Djajic V.
        • Bajic D.
        • Djokovic A.
        Assessment of autonomic nervous system dysfunction in the early phase of infection with SARS-CoV-2 virus.
        Front Neurosci. 2021; 15: 1-14https://doi.org/10.3389/fnins.2021.640835
        • Buoite Stella A.
        • Furlanis G.
        • Frezza N.A.
        • Valentinotti R.
        • Ajcevic M.
        • Manganotti P.
        Autonomic dysfunction in post-COVID patients with and without neurological symptoms: a prospective multidomain observational study.
        J Neurol. 2022; 269: 587-596https://doi.org/10.1007/s00415-021-10735-y
        • Scala I.
        • Rizzo P.A.
        • Bellavia S.
        • Brunetti V.
        • Col F.
        • Broccolini A.
        • et al.
        Autonomic dysfunction during acute SARS-CoV-2 infection : a systematic review.
        J Clin Med. 2022; 11: 1-26https://doi.org/10.3390/jcm11133883
        • Fudim M.
        • Qadri Y.J.
        • Ghadimi K.
        • MacLeod D.B.
        • Molinger J.
        • Piccini J.P.
        • et al.
        Implications for neuromodulation therapy to control inflammation and related organ dysfunction in COVID-19.
        J Cardiovasc Transl Res. 2020; 13: 894-899https://doi.org/10.1007/s12265-020-10031-6
        • Andrade S.M.
        • Cecília de Araújo Silvestre M.
        • Tenório de França E.É.
        • Bezerra Sales Queiroz M.H.
        • de Jesus Santana K.
        • Lima Holmes Madruga M.L.
        • et al.
        Efficacy and safety of HD-tDCS and respiratory rehabilitation for critically ill patients with COVID-19 the HD-RECOVERY randomized clinical trial.
        Brain Stimul. 2022; 15: 780-788https://doi.org/10.1016/j.brs.2022.05.006
        • Wu Z.
        • McGoogan J.M.
        Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China.
        JAMA. 2020; 323: 1239https://doi.org/10.1001/jama.2020.2648
        • Marvar P.J.
        • Harrison D.G.
        Inflammation, immunity and the autonomic nervous system. Prim. Auton. Nerv. Syst.
        Elsevier, 2012: 325-329https://doi.org/10.1016/B978-0-12-386525-0.00067-6
        • Tornero C.
        • Pastor E.
        • Garzando M del M.
        • Orduña J.
        • Forner M.J.
        • Bocigas I.
        • et al.
        Non-invasive vagus nerve stimulation for COVID-19: results from a randomized controlled trial (SAVIOR I).
        Front Neurol. 2022; 13https://doi.org/10.3389/fneur.2022.820864
        • Iseger T.A.
        • van Bueren N.E.R.
        • Kenemans J.L.
        • Gevirtz R.
        • Arns M.
        A frontal-vagal network theory for Major Depressive Disorder: implications for optimizing neuromodulation techniques.
        Brain Stimul. 2020; 13: 1-9https://doi.org/10.1016/j.brs.2019.10.006
        • Makovac E.
        • Thayer J.F.
        • Ottaviani C.
        A meta-analysis of non-invasive brain stimulation and autonomic functioning: implications for brain-heart pathways to cardiovascular disease.
        Neurosci Biobehav Rev. 2017; 74: 330-341https://doi.org/10.1016/j.neubiorev.2016.05.001
        • Casanova M.F.
        • Hensley M.K.
        • Sokhadze E.M.
        • El-Baz A.S.
        • Wang Y.
        • Li X.
        • et al.
        Effects of weekly low-frequency rTMS on autonomic measures in children with autism spectrum disorder.
        Front Hum Neurosci. 2014; 8https://doi.org/10.3389/fnhum.2014.00851
        • Nikolin S.
        • Boonstra T.W.
        • Loo C.K.
        • Martin D.
        Combined effect of prefrontal transcranial direct current stimulation and a working memory task on heart rate variability.
        PLoS One. 2017; 12e0181833https://doi.org/10.1371/journal.pone.0181833
        • Brunoni A.R.
        • Vanderhasselt M.-A.
        • Boggio P.S.
        • Fregni F.
        • Dantas E.M.
        • Mill J.G.
        • et al.
        Polarity- and valence-dependent effects of prefrontal transcranial direct current stimulation on heart rate variability and salivary cortisol.
        Psychoneuroendocrinology. 2013; 38: 58-66https://doi.org/10.1016/j.psyneuen.2012.04.020
        • Lin Y.-Y.
        • Chang C.-C.
        • Huang C.C.-Y.
        • Tzeng N.-S.
        • Kao Y.-C.
        • Chang H.-A.
        Efficacy and neurophysiological predictors of treatment response of adjunct bifrontal transcranial direct current stimulation (tDCS) in treating unipolar and bipolar depression.
        J Affect Disord. 2021; 280: 295-304https://doi.org/10.1016/j.jad.2020.11.030
        • Kemp A.H.
        • Quintana D.S.
        • Gray M.A.
        • Felmingham K.L.
        • Brown K.
        • Gatt J.M.
        Impact of depression and antidepressant treatment on heart rate variability: a review and meta-analysis.
        Biol Psychiatr. 2010; 67: 1067-1074https://doi.org/10.1016/j.biopsych.2009.12.012
        • Lakens D.
        Sample size justification.
        Collabra Psychol. 2022; 8https://doi.org/10.1525/collabra.33267
        • Schulz K.F.
        CONSORT 2010 statement: updated guidelines for reporting parallel group randomized trials.
        Ann Intern Med. 2010; 152: 726https://doi.org/10.7326/0003-4819-152-11-201006010-00232
        • World Health Organization
        COVID-19 therapeutic trial synopsis.
        2020
        • Nuwer M.R.
        • Comi G.
        • Emerson R.
        • Fuglsang-Frederiksen A.
        • Guerit M.
        • Hinrichs H.
        • et al.
        IFCN standards for digital recording of clinical EEG.
        Electroencephalogr Clin Neurophysiol. 1998; 106: 259-261https://doi.org/10.1016/S0013-4694(97)00106-5
        • Cerqueira J.J.
        • Almeida O.F.X.
        • Sousa N.
        The stressed prefrontal cortex.
        Left? Right! Brain Behav Immun. 2008; 22: 630-638https://doi.org/10.1016/j.bbi.2008.01.005
        • Bikson M.
        • Grossman P.
        • Thomas C.
        • Zannou A.L.
        • Jiang J.
        • Adnan T.
        • et al.
        Safety of transcranial direct current stimulation: evidence based update 2016.
        Brain Stimul. 2016; 9: 641-661https://doi.org/10.1016/j.brs.2016.06.004
        • Gandiga P.C.
        • Hummel F.C.
        • Cohen L.G.
        Transcranial DC stimulation (tDCS): a tool for double-blind sham-controlled clinical studies in brain stimulation.
        Clin Neurophysiol. 2006; 117: 845-850https://doi.org/10.1016/j.clinph.2005.12.003
        • Bikson M.
        • Hanlon C.A.
        • Woods A.J.
        • Gillick B.T.
        • Charvet L.
        • Lamm C.
        • et al.
        Guidelines for TMS/tES clinical services and research through the COVID-19 pandemic.
        Brain Stimul. 2020; 13: 1124-1149https://doi.org/10.1016/j.brs.2020.05.010
        • Brunoni A.R.
        • Amadera J.
        • Berbel B.
        • Volz M.S.
        • Rizzerio B.G.
        • Fregni F.
        A systematic review on reporting and assessment of adverse effects associated with transcranial direct current stimulation.
        Int J Neuropsychopharmacol. 2011; 14: 1133-1145https://doi.org/10.1017/S1461145710001690
        • Gomes-Oliveira M.H.
        • Gorenstein C.
        • Neto F.L.
        • Andrade L.H.
        • Wang Y.P.
        Validação da versão Brasileira em Português do Inventário de Depressão de Beck-II numa amostra da comunidade.
        Rev Bras Psiquiatr. 2012; 34: 389-394https://doi.org/10.1016/j.rbp.2012.03.005
      1. Beck AT, Beamesderfer A. Assessment of Depression: The Depression Inventory, n.d, pp. 151–169. https://doi.org/10.1159/000395074.

        • Beck A.T.
        • Epstein N.
        • Brown G.
        • Steer R.A.
        An inventory for measuring clinical anxiety: psychometric properties.
        J Consult Clin Psychol. 1988; 56: 893-897https://doi.org/10.1037/0022-006X.56.6.893
        • Hartikainen J.E.K.
        • Tahvanainen K.U.O.
        • Kuusela T.A.
        Short-term measurement of heart rate variability. Clin. Guid. To card. Auton. Tests.
        Springer Netherlands, Dordrecht1998: 149-176https://doi.org/10.1007/978-94-017-1057-2_6
        • Tarvainen M.P.
        • Niskanen J.-P.
        • Lipponen J.A.
        • Ranta-aho P.O.
        • Karjalainen P.A.
        Kubios HRV – heart rate variability analysis software.
        Comput Methods Progr Biomed. 2014; 113: 210-220https://doi.org/10.1016/j.cmpb.2013.07.024
        • European Society of Cardiology H.
        Rate variability, of measurement, physiological interpretation and clinical use, American TF of TES of C and TN, Electrophysiology S of P and. Guidelines Heart rate variability.
        Eur Heart J. 1996; 17: 354-381https://doi.org/10.1161/01.CIR.93.5.1043
        • O'Connell N.S.
        • Dai L.
        • Jiang Y.
        • Speiser J.L.
        • Ward R.
        • Wei W.
        • et al.
        Methods for analysis of pre-post data in clinical research: a comparison of five common methods.
        J Biometrics Biostat. 2017; 8https://doi.org/10.4172/2155-6180.1000334
        • Müller-Dahlhaus J.F.M.
        • Orekhov Y.
        • Liu Y.
        • Ziemann U.
        Interindividual variability and age-dependency of motor cortical plasticity induced by paired associative stimulation.
        Exp Brain Res. 2008; 187: 467-475https://doi.org/10.1007/s00221-008-1319-7
        • Fathi D.
        • Ueki Y.
        • Mima T.
        • Koganemaru S.
        • Nagamine T.
        • Tawfik A.
        • et al.
        Effects of aging on the human motor cortical plasticity studied by paired associative stimulation.
        Clin Neurophysiol. 2010; 121: 90-93https://doi.org/10.1016/j.clinph.2009.07.048
        • Nunan D.
        • Sandercock G.R.H.
        • Brodie D.A.
        A quantitative systematic review of normal values for short-term heart rate variability in healthy adults.
        PACE - Pacing Clin Electrophysiol. 2010; 33: 1407-1417https://doi.org/10.1111/j.1540-8159.2010.02841.x
        • Faul F.
        • Erdfelder E.
        • Lang A.-G.
        • Buchner A.
        G*Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences.
        Behav Res Methods. 2007; 39: 175-191https://doi.org/10.3758/BF03193146
        • Salehinejad M.A.
        • Ghanavai E.
        • Rostami R.
        • Nejati V.
        Cognitive control dysfunction in emotion dysregulation and psychopathology of major depression (MD): evidence from transcranial brain stimulation of the dorsolateral prefrontal cortex (DLPFC).
        J Affect Disord. 2017; 210: 241-248https://doi.org/10.1016/j.jad.2016.12.036
        • Alizadehgoradel J.
        • Nejati V.
        • Sadeghi Movahed F.
        • Imani S.
        • Taherifard M.
        • Mosayebi-Samani M.
        • et al.
        Repeated stimulation of the dorsolateral-prefrontal cortex improves executive dysfunctions and craving in drug addiction: a randomized, double-blind, parallel-group study.
        Brain Stimul. 2020; 13: 582-593https://doi.org/10.1016/j.brs.2019.12.028
        • Dundas J.E.
        • Thickbroom G.W.
        • Mastaglia F.L.
        Perception of comfort during transcranial DC stimulation: effect of NaCl solution concentration applied to sponge electrodes.
        Clin Neurophysiol. 2007; 118: 1166-1170https://doi.org/10.1016/j.clinph.2007.01.010
        • Garnsworthy R.K.
        • Gully R.L.
        • Kenins P.
        • Westerman R.A.
        Transcutaneous electrical stimulation and the sensation of prickle.
        J Neurophysiol. 1988; 59: 1116-1127https://doi.org/10.1152/jn.1988.59.4.1116
        • Poreisz C.
        • Boros K.
        • Antal A.
        • Paulus W.
        Safety aspects of transcranial direct current stimulation concerning healthy subjects and patients.
        Brain Res Bull. 2007; 72: 208-214https://doi.org/10.1016/j.brainresbull.2007.01.004
        • Pagani M.
        • Lombardi F.
        • Guzzetti S.
        • Rimoldi O.
        • Furlan R.
        • Pizzinelli P.
        • et al.
        Power spectral analysis of heart rate and arterial pressure variabilities as a marker of sympatho-vagal interaction in man and conscious dog.
        Circ Res. 1986; 59: 178-193https://doi.org/10.1161/01.RES.59.2.178
        • Malliani A.
        • Pagani M.
        • Lombardi F.
        • Cerutti S.
        Cardiovascular neural regulation explored in the frequency domain.
        Circulation. 1991; 84: 482-492https://doi.org/10.1161/01.CIR.84.2.482
        • Pomeranz B.
        • Macaulay R.J.
        • Caudill M.A.
        • Kutz I.
        • Adam D.
        • Gordon D.
        • et al.
        Assessment of autonomic function in humans by heart rate spectral analysis.
        Am J Physiol. 1985; 248: H151-H153
        • Akselrod S.
        • Gordon D.
        • Ubel F.
        • Shannon D.C.
        • Barger A.C.
        • Cohen R.J.
        Power spectrum analysis of heart rate fluctuation: a quantitative probe of beat-to-beat cardiovascular control.
        Science. 1981; 213: 220-222
        • Electrophysiology TF of the ES of C the NA
        Heart rate variability.
        Circulation. 1996; 93: 1043-1065https://doi.org/10.1161/01.CIR.93.5.1043
        • Cohen J.
        Statistical power analysis for the behavioral sciences.
        second ed. Lawrence Erlbaum Associates, Publishers, Hillsdale, Nj1988
        • Wasserstein R.L.
        • Schirm A.L.
        • Lazar N.A.
        Moving to a world beyond “ p < 0.05.”.
        Am Statistician. 2019; 73: 1-19https://doi.org/10.1080/00031305.2019.1583913
        • Cheshire W.P.
        • Freeman R.
        • Gibbons C.H.
        • Cortelli P.
        • Wenning G.K.
        • Hilz M.J.
        • et al.
        Electrodiagnostic assessment of the autonomic nervous system: a consensus statement endorsed by the American autonomic society, American academy of neurology, and the international federation of clinical neurophysiology.
        Clin Neurophysiol. 2021; 132: 666-682https://doi.org/10.1016/j.clinph.2020.11.024
        • Fox M.D.
        • Buckner R.L.
        • White M.P.
        • Greicius M.D.
        • Pascual-leone A.
        Efficacy of transcranial magnetic stimulation targets for depression is related to intrinsic functional connectivity with the subgenual cingulate.
        Biol Psychiatr. 2012; 72: 595-603https://doi.org/10.1016/j.biopsych.2012.04.028
        • Iseger T.A.
        • Padberg F.
        • Kenemans J.L.
        • Gevirtz R.
        • Arns M.
        Neuro-Cardiac-Guided TMS (NCG-TMS): probing DLPFC-sgACC-vagus nerve connectivity using heart rate e First results.
        Brain Stimul. 2017; 10: 1006-1008https://doi.org/10.1016/j.brs.2017.05.002
        • Kessler S.K.
        • Turkeltaub P.E.
        • Benson J.G.
        • Hamilton R.H.
        Differences in the experience of active and sham transcranial direct current stimulation.
        Brain Stimul. 2012; 5: 155-162https://doi.org/10.1016/j.brs.2011.02.007
        • van Boekholdt L.
        • Kerstens S.
        • Khatoun A.
        • Asamoah B.
        • Mc Laughlin M.
        tDCS peripheral nerve stimulation: a neglected mode of action?.
        Mol Psychiatr. 2021; 26: 456-461https://doi.org/10.1038/s41380-020-00962-6
        • Faria D.
        • Moll-Bernardes R.J.
        • Testa L.
        • Moniz C.M.V.
        • Rodrigues E.C.
        • Rodrigues A.G.
        • et al.
        Sympathetic neural overdrive, aortic stiffening, endothelial dysfunction, and impaired exercise capacity in severe COVID-19 survivors: a mid-term study of cardiovascular sequelae.
        Hypertension. 2023; 80: 470-481https://doi.org/10.1161/HYPERTENSIONAHA.122.19958
        • Porzionato A.
        • Emmi A.
        • Barbon S.
        • Boscolo‐Berto R.
        • Stecco C.
        • Stocco E.
        • et al.
        Sympathetic activation: a potential link between comorbidities and COVID‐19.
        FEBS J. 2020; 287: 3681-3688https://doi.org/10.1111/febs.15481
        • Junarta J.
        • Riley J.M.
        • Pavri B.B.
        Describing heart rate variability in patients with chronic atrial fibrillation during hospitalization for COVID‐19.
        J Arrhythmia. 2021; 37: 893-898https://doi.org/10.1002/joa3.12569
        • Al-kuraishy H.M.
        • Al-Gareeb A.I.
        • Qusti S.
        • Alshammari E.M.
        • Gyebi G.A.
        • Batiha G.E.-S.
        Covid-19-Induced dysautonomia: a menace of sympathetic storm.
        ASN Neuro. 2021; 13175909142110576https://doi.org/10.1177/17590914211057635
        • Thayer J.F.
        • Siegle G.J.
        Neurovisceral integration in cardiac and emotional regulation.
        IEEE Eng Med Biol Mag. 2002; 21: 24-29https://doi.org/10.1109/MEMB.2002.1032635
        • Vai B.
        • Mazza M.G.
        • Delli Colli C.
        • Foiselle M.
        • Allen B.
        • Benedetti F.
        • et al.
        Mental disorders and risk of COVID-19-related mortality, hospitalisation, and intensive care unit admission: a systematic review and meta-analysis.
        Lancet Psychiatr. 2021; 8: 797-812https://doi.org/10.1016/S2215-0366(21)00232-7
        • Premraj L.
        • Kannapadi N.V.
        • Briggs J.
        • Seal S.M.
        • Battaglini D.
        • Fanning J.
        • et al.
        Mid and long-term neurological and neuropsychiatric manifestations of post-COVID-19 syndrome: a meta-analysis.
        J Neurol Sci. 2022; 434120162https://doi.org/10.1016/j.jns.2022.120162
        • Mazza M.G.
        • De Lorenzo R.
        • Conte C.
        • Poletti S.
        • Vai B.
        • Bollettini I.
        • et al.
        Anxiety and depression in COVID-19 survivors: role of inflammatory and clinical predictors.
        Brain Behav Immun. 2020; 89: 594-600https://doi.org/10.1016/j.bbi.2020.07.037
        • Pavlov V.A.
        • Tracey K.J.
        The vagus nerve and the inflammatory reflex—linking immunity and metabolism.
        Nat Rev Endocrinol. 2012; 8: 743-754https://doi.org/10.1038/nrendo.2012.189
        • Bonaz B.
        • Picq C.
        • Sinniger V.
        • Mayol J.F.
        • Clarençon D.
        Vagus nerve stimulation: from epilepsy to the cholinergic anti-inflammatory pathway.
        Neuro Gastroenterol Motil. 2013; 25: 208-221https://doi.org/10.1111/nmo.12076
        • Bonaz B.
        • Sinniger V.
        • Pellissier S.
        Vagus nerve stimulation: a new promising therapeutic tool in inflammatory bowel disease.
        J Intern Med. 2017; 282: 46-63https://doi.org/10.1111/joim.12611
        • Borovikova L.V.
        • Ivanova S.
        • Zhang M.
        • Yang H.
        • Botchkina G.I.
        • Watkins L.R.
        • et al.
        Vagus nerve stimulation attenuates the systemic inflammatory response to endotoxin.
        Nature. 2000; 405: 458-462https://doi.org/10.1038/35013070
        • Tracey K.J.
        Physiology and immunology of the cholinergic antiinflammatory pathway.
        J Clin Invest. 2007; 117: 289-296https://doi.org/10.1172/JCI30555
        • Stavrakis S.
        • Humphrey M.B.
        • Scherlag B.
        • Iftikhar O.
        • Parwani P.
        • Abbas M.
        • et al.
        Low-level vagus nerve stimulation suppresses post-operative atrial fibrillation and inflammation.
        JACC Clin Electrophysiol. 2017; 3: 929-938https://doi.org/10.1016/j.jacep.2017.02.019
        • Mondal B.
        • Choudhury S.
        • Simon B.
        • Baker M.R.
        • Kumar H.
        Noninvasive vagus nerve stimulation improves gait and reduces freezing of gait in Parkinson's disease.
        Mov Disord. 2019; 34: 917-918https://doi.org/10.1002/mds.27662
        • Bremner J.D.
        • Gurel N.Z.
        • Jiao Y.
        • Wittbrodt M.T.
        • Levantsevych O.M.
        • Huang M.
        • et al.
        Transcutaneous vagal nerve stimulation blocks stress-induced activation of Interleukin-6 and interferon-γ in posttraumatic stress disorder: a double-blind, randomized, sham-controlled trial.
        Brain, Behav Immun - Heal. 2020; 9100138https://doi.org/10.1016/j.bbih.2020.100138
        • Boezaart A.P.
        • Botha D.A.
        Treatment of stage 3 COVID-19 with transcutaneous auricular vagus nerve stimulation drastically reduces interleukin-6 blood levels: a report on two cases.
        Neuromodulation Technol Neural Interface. 2021; 24: 166-167https://doi.org/10.1111/ner.13293
        • Sobel J.A.
        • Levy J.
        • Almog R.
        • Reiner-Benaim A.
        • Miller A.
        • Eytan D.
        • et al.
        Descriptive characteristics of continuous oximetry measurement in moderate to severe covid-19 patients.
        Sci Rep. 2023; 13: 442https://doi.org/10.1038/s41598-022-27342-0
        • Spaccarotella C.
        • Polimeni A.
        • Mancuso C.
        • Pelaia G.
        • Esposito G.
        • Indolfi C.
        Assessment of non-invasive measurements of oxygen saturation and heart rate with an apple smartwatch: comparison with a standard pulse oximeter.
        J Clin Med. 2022; 11: 1467https://doi.org/10.3390/jcm11061467
        • Motta L.P.
        • da Silva P.P.F.
        • Borguezan B.M.
        • do Amaral JLM.
        • Milagres L.G.
        • Bóia M.N.
        • et al.
        An emergency system for monitoring pulse oximetry, peak expiratory flow, and body temperature of patients with COVID-19 at home: development and preliminary application.
        PLoS One. 2021; 16e0247635https://doi.org/10.1371/journal.pone.0247635
        • Flierl M.A.
        • Rittirsch D.
        • Nadeau B.A.
        • Sarma J.V.
        • Day D.E.
        • Lentsch A.B.
        • et al.
        Upregulation of phagocyte-derived catecholamines augments the acute inflammatory response.
        PLoS One. 2009; 4e4414https://doi.org/10.1371/journal.pone.0004414
        • Hyman A.L.
        • Kadowitz P.J.
        Evidence for existence of postjunctional alpha 1- and alpha 2-adrenoceptors in cat pulmonary vascular bed.
        Am J Physiol Circ Physiol. 1985; 249: H891-H898https://doi.org/10.1152/ajpheart.1985.249.4.H891
        • Hyman A.L.
        • Kadowitz P.J.
        Enhancement of alpha- and beta-adrenoceptor responses by elevations in vascular tone in pulmonary circulation.
        Am J Physiol Circ Physiol. 1986; 250: H1109-H1116https://doi.org/10.1152/ajpheart.1986.250.6.H1109
        • Nandiwada P.A.
        • Hyman A.L.
        • Kadowitz P.J.
        Pulmonary vasodilator responses to vagal stimulation and acetylcholine in the cat.
        Circ Res. 1983; 53: 86-95https://doi.org/10.1161/01.RES.53.1.86
        • Yang X.
        • Yu Y.
        • Xu J.
        • Shu H.
        • Xia J.
        • Liu H.
        • et al.
        Clinical course and outcomes of critically ill patients with SARS-CoV-2 pneumonia in Wuhan, China: a single-centered, retrospective, observational study.
        Lancet Respir Med. 2020; 8: 475-481https://doi.org/10.1016/S2213-2600(20)30079-5
        • Myers L.C.
        • Parodi S.M.
        • Escobar G.J.
        • Liu V.X.
        Characteristics of hospitalized adults with COVID-19 in an integrated health care system in California.
        JAMA. 2020; 323: 2195https://doi.org/10.1001/jama.2020.7202
        • Tarn J.
        • Legg S.
        • Mitchell S.
        • Simon B.
        • Ng W.-F.
        The effects of noninvasive vagus nerve stimulation on fatigue and immune responses in patients with primary Sjögren’s syndrome.
        Neuromodulation Technol Neural Interface. 2019; 22: 580-585https://doi.org/10.1111/ner.12879
        • Staats P.
        • Emala C.
        • Simon B.
        • Errico J.P.
        Neurostimulation for asthma. Neuromodulation.
        Elsevier, 2018: 1339-1345https://doi.org/10.1016/B978-0-12-805353-9.00111-X
        • Kleiger R.E.
        • Stein P.K.
        • Bigger J.T.
        Heart rate variability: measurement and clinical utility.
        Ann Noninvasive Electrocardiol. 2005; 10: 88-101https://doi.org/10.1111/j.1542-474X.2005.10101.x
        • Silva F.C.
        • Paiva F.A.
        • Müller-Ribeiro F.C.
        • Caldeira H.M.A.
        • Fontes M.A.P.
        • de Menezes R.C.A.
        • et al.
        Chronic treatment with ivabradine does not affect cardiovascular autonomic control in rats.
        Front Physiol. 2016; 7https://doi.org/10.3389/fphys.2016.00305
        • Dias da Silva V.J.
        • Tobaldini E.
        • Rocchetti M.
        • Wu M.A.
        • Malfatto G.
        • Montano N.
        • et al.
        Modulation of sympathetic activity and heart rate variability by ivabradine.
        Cardiovasc Res. 2015; 108: 31-38https://doi.org/10.1093/cvr/cvv180