The effect of transcutaneous auricular vagus nerve stimulation on cardiovascular function in subarachnoid hemorrhage patients: A randomized trial
eLife Assessment
The authors provide a solid set of data supporting the safety of transcutaneous auricular vagal nerve stimulation on cardiovascular parameters in the acute setting of critically ill patients presenting with subarachnoid hemorrhage. This important study also suggests a promising effect on autonomic balance.
https://doi.org/10.7554/eLife.100088.3.sa0Important: Findings that have theoretical or practical implications beyond a single subfield
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Solid: Methods, data and analyses broadly support the claims with only minor weaknesses
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Abstract
Background:
Subarachnoid hemorrhage (SAH) is characterized by intense central inflammation, leading to substantial post-hemorrhagic complications such as vasospasm and delayed cerebral ischemia. Given the anti-inflammatory effect of transcutaneous auricular vagus nerve stimulation (taVNS) and its ability to promote brain plasticity, taVNS has emerged as a promising therapeutic option for SAH patients. However, the effects of taVNS on cardiovascular dynamics in critically ill patients, like those with SAH, have not yet been investigated. Given the association between cardiac complications and elevated risk of poor clinical outcomes after SAH, it is essential to characterize the cardiovascular effects of taVNS to ensure this approach is safe in this fragile population. Therefore, this study assessed the impact of both acute and repetitive taVNS on cardiovascular function.
Methods:
In this randomized clinical trial, 24 SAH patients were assigned to either a taVNS treatment or a sham treatment group. During their stay in the intensive care unit, we monitored patient electrocardiogram readings and vital signs. We compared long-term changes in heart rate, heart rate variability (HRV), QT interval, and blood pressure between the two groups. Additionally, we assessed the effects of acute taVNS by comparing cardiovascular metrics before, during, and after the intervention. We also explored acute cardiovascular biomarkers in patients exhibiting clinical improvement.
Results:
We found that repetitive taVNS did not significantly alter heart rate, QT interval, blood pressure, or intracranial pressure (ICP). However, repetitive taVNS increased overall HRV and parasympathetic activity compared to the sham treatment. The increase in parasympathetic activity was most pronounced from 2 to 4 days after initial treatment (Cohen’s d = 0.50). Acutely, taVNS increased heart rate, blood pressure, and peripheral perfusion index without affecting the corrected QT interval, ICP, or HRV. The acute post-treatment elevation in heart rate was more pronounced in patients who experienced a decrease of more than one point in their modified Rankin Score at the time of discharge.
Conclusions:
Our study found that taVNS treatment did not induce adverse cardiovascular effects, such as bradycardia or QT prolongation, supporting its development as a safe immunomodulatory treatment approach for SAH patients. The observed acute increase in heart rate after taVNS treatment may serve as a biomarker for SAH patients who could derive greater benefit from this treatment.
Funding:
The American Association of Neurological Surgeons (ALH), The Aneurysm and AVM Foundation (ALH), The National Institutes of Health R01-EB026439, P41-EB018783, U24-NS109103, R21-NS128307 (ECL, PB), McDonnell Center for Systems Neuroscience (ECL, PB), and Fondazione Neurone (PB).
Clinical trial number:
Introduction
Subarachnoid hemorrhage (SAH) is a devastating subtype of stroke that represents a significant global health burden and causes permanent disability in approximately 30% of survivors (D’Souza, 2015; Weir, 2002; Zhang et al., 2023). Early brain injury can occur within the first 24–48 hr after ictus, which involves a cascade of elevated intracranial pressure (ICP) and a subsequent drop of cerebral perfusion (Schneider et al., 2018). Systemic and local inflammation, cerebral edema, blood–brain barrier disruption, sympathetic nervous system activation, autoregulatory failure, microthrombosis, spreading depolarizations, and inflammation have all been observed during this period (Macdonald et al., 2012; Budohoski et al., 2013). These biological processes result in the inability of cerebral perfusion to match metabolic demands, leading to secondary brain injury and delayed cerebral ischemia that typically occurs between 5 and 14 days after the SAH (Provencio, 2013; van Gijn et al., 2007; Tracey, 2002; Lv et al., 2018). Delayed cerebral ischemia and deleterious inflammation are major predictors of poor outcomes and morbidity. The autonomic nervous system (ANS), comprising the sympathetic and the parasympathetic nervous system, plays a critical role in maintaining physiological homeostasis. SAH is believed to cause sympathetic predominance, which plays a key role in the development of cerebral vasospasm, renders patients more susceptible to non-neurological complications, and exacerbates deleterious inflammatory processes (Naredi et al., 2000).
Numerous interventions have been explored to address the complex pathologies of SAH that contribute to secondary brain injury, aiming to improve patient outcomes (Provencio and Vora, 2005). Transcutaneous auricular vagus nerve stimulation (taVNS) is one of the most promising therapeutic options, as recent studies have demonstrated its efficacy in reducing inflammation, improving autonomic balance, and enhancing brain plasticity (Wu et al., 2023; Bonaz et al., 2016; Meyers et al., 2018; Provencio and Vora, 2005; Huguenard et al., 2023). The auricular branch of the vagus nerve is a sensory nerve that innervates the external ear, including the cymba concha. Stimulating the auricular branch of the vagus nerve has been shown to activate the same brain regions as cervical vagus nerve stimulation (Frangos et al., 2015). Specifically, taVNS mediates cholinergic signaling and regulates pro-inflammatory responses via the inflammatory reflex (Figure 1A; Sahn et al., 2023; Pavlov and Tracey, 2012) In this reflex, inflammatory mediators such as cytokines trigger afferent vagus nerve signaling. This afferent signal then prompts an efferent response from the vagus nerve that acts to reduce the production of pro-inflammatory cytokines (Tynan et al., 2022).
The vagus nerve also mediates cardiovascular function by regulating the autonomic system and metabolic homeostasis (Figure 1A; Keute et al., 2021). Theoretically, taVNS increases parasympathetic activity, which can be measured as increased heart rate variability (HRV). While some animal studies have reported a potential risk of bradycardia and decreased blood pressure associated with vagus nerve stimulation, two reviews of human studies have considered the cardiovascular effects of taVNS generally safe, with adverse effects reported only in patients with pre-existing heart diseases (Naggar et al., 2014; Hua et al., 2023; Kim et al., 2022). However, its cardiovascular effect in SAH patients is largely unknown. Given that critically ill patients, such as those with SAH, are extremely vulnerable to cardiovascular complications, it is essential to thoroughly examine the cardiovascular implications of taVNS to ensure its safety in this fragile population. This is particularly notable as cardiovascular abnormalities following SAH, such as prolonged elevated heart rate and QT prolongation, are associated with an increased risk of poor outcomes (Norberg et al., 2018; Schmidt et al., 2014; Zhang and Qi, 2016; van der Bilt et al., 2009). However, our limited understanding of these effects constitutes a significant barrier, preventing the advancement of taVNS from a promising therapeutic approach to an established clinical treatment for SAH. To address this gap, we assessed the effects of acute and repetitive taVNS on cardiovascular function based on electrocardiogram (ECG) and other monitored vital signs from SAH patients in the intensive care unit (ICU). The current study is part of the NAVSaH trial (NCT04557618) and focuses on the trial’s secondary outcomes, including heart rate, QT interval, HRV, and blood pressure (Huguenard et al., 2024a). This interim analysis aims to evaluate the cardiovascular safety of the taVNS protocol and to provide insights that will inform the application of taVNS in SAH patients. The primary outcomes of this trial, including change in the inflammatory cytokine tumor necrosis factor-alpha (TNF-α) and rate of radiographic vasospasm, are available as a pre-print and currently under review (Huguenard et al., 2024b). Based on a meta-analysis, most aversive effects were seen in repeated sessions lasting 60 min or more; therefore, we hypothesized that repetitive taVNS increased HRV but did not cause bradycardia and QT prolongation (Kim et al., 2022). To test this hypothesis, we compared changes in cardiovascular metrics at the phase of early brain injury (within 72 hr) and at the phase when delayed cerebral ischemia develops (after Day 4) between patients receiving taVNS and sham treatments. Root mean square of successive differences (RMSSD) and the standard deviation of normal RR intervals (SDNN) are two commonly used HRV metrics. RMSSD indicates parasympathetic activity, while lower SDNN is associated with increased cardiac risk (Kleiger et al., 2005; Shaffer and Ginsberg, 2017). Providing effective taVNS treatment modulates the autonomic system, we propose that heart rate or HRV following acute taVNS could inform which SAH patients may experience the most clinical benefit from taVNS treatment. To explore this possibility, we correlated the changes in heart rate and HRV following acute taVNS treatment and changes in the modified Rankin Score (mRS), which measures the degree of disability or dependence in the daily activities of people suffering from neurological disability.
Results
Twenty-four participants were randomized to receive the taVNS (N = 11) or Sham (N = 13) treatment (Table 1, Figure 2—figure supplement 2). Appendix 3—table 1 shows the clinical characteristics of the two treatment groups. The participants, the medical team who dictated all management decisions for the patient’s SAH, and the outcomes assessors who assigned mRS at admission and discharge were blinded to the treatment. The structure of this study is shown in Figure 1B. Following randomization, enrolled patients underwent 20 min of either taVNS or sham stimulation twice daily during their stay in the ICU. This treatment schedule was informed by findings from Addorisio et al., where a 5-min taVNS protocol was administered twice daily to patients with rheumatoid arthritis for 2 days (Addorisio et al., 2019). Their study found that circulating C-reactive protein levels significantly reduced after 2 days of treatment but returned to baseline at the second clinical assessment by Day 7. Given the high inflammatory state associated with SAH and our intention to maintain a steady reduction in inflammation, we decided to extend the treatment duration to 20 min per session. During treatment periods, a portable transcutaneous electrical nerve stimulation (TENS) device (TENS 7000 Digital TENS Unit, Compass Health Brands, OH, USA) was connected to the patient’s left ear using two ear clips (Figure 1C, D). For taVNS treatments, these ear clips were placed along the concha of the ear, while for sham treatments, the clips were placed along the earlobe to avoid stimulation of the auricular vagus nerve from tactile pressure (Figure 1C). For the taVNS group, stimulation parameters were selected based on values reported in prior studies that sought to maximize vagus somatosensory evoked potentials while avoiding the perception of pain: 20 Hz frequency, 250 µs pulse width, and 0.4 mA intensity (de Gurtubay et al., 2021). The stimulation parameters were designed to be imperceptible to the patients, and there were no reports of detection of taVNS, suggesting the success of the blinding. No electrical current was delivered during sham treatments. Please see Huguenard et al., 2024a for a detailed protocol of this study.
Effects of repetitive taVNS on cardiac function
A study has shown that 15 min of taVNS reduced sympathetic activity in healthy individuals, with effects that persist during the recovery period (Clancy et al., 2014). This finding suggests that taVNS may exert long-term effects on cardiovascular function. Therefore, we investigated whether repetitive taVNS treatment affects heart rate and QT interval, key indicators of bradycardia or QT prolongation, using 24-hr ECG recording. We found no significant differences in heart rate between groups (Figure 2D, Mann–Whitney U test, N(taVNS) = 94, N(Sham) = 95, p-value = 0.69, Cohen’s d = −0.01, W-statistics = 4317, power = 0.93). Changes in heart rate from Day 1 were equivalent between groups (two-tailed equivalence tests, p[lower threshold] = 0.006, test statistics[lower threshold] = 2.53; p[lower threshold] = 0.004, test statistics[lower threshold] = −2.72, N(VNS) = 94, N(VNS) = 95). We further confirmed that changes in heart rate were similar between treatment groups following SAH (Figure 2D, |Cohen’s d|<0.2 for Days 2–4, 5–8, 8–10, and 11–13). Moreover, changes in corrected QT interval from Day 1 were significantly higher in the Sham group compared to the taVNS group (Figure 2E, Mann–Whitney U test, N(taVNS) = 94, N(Sham) = 95, p-value <0.001, Cohen’s d = −0.57). Similarly, uncorrected QT intervals from Day 1 were higher in the Sham group (Figure 2—figure supplement 1, Cohen’s d = −0.42). After the phase of early brain injury, the mean and median corrected QT interval were lower in the taVNS group with large effect sizes (Figure 2E, |Cohen’s d| >0.5). To ensure that repetitive taVNS did not lead to QT prolongation outside the stimulation period, we calculated the percentage of prolonged QT intervals. Prolonged QT intervals were defined as corrected QT interval ≥500 ms. We found that changes in prolonged QT intervals percentage from Day 1 were higher in the Sham group (Figure 2F, Mann–Whitney U test, N(taVNS) = 94, N(Sham) = 95, p-value <0.001, Cohen’s d = −0.72).
Subsequently, we investigated the effect of taVNS treatment on RMSSD and SDNN. We found that changes in SDNN using Day 1 as baseline were not significantly different between the treatment groups (T-test, N(taVNS) = 94, N(Sham) = 95, p = 0.479, Cohen’s d = 0.10, t statistics = 0.71, Figure 3A). Changes in RMSSD using Day 1 as baseline were significantly higher in the taVNS treatment group (T-test, N(taVNS) = 94, N(Sham) = 95, Bonferroni-corrected p = 0.025, Cohen’s d = 0.42, t statistics = 2.91, Figure 3B). We further studied the effects of taVNS in different phases following SAH. Figure 3A-B show the changes in SDNN and RMSSD in bins of 3 days for the two treatment groups. The taVNS treatment increased RMSSD over the course of the treatment, with a smaller effect size (Cohen’s d = 0.29) observed between Days 2 and 4, corresponding to the early brain injury phase, and large effect sizes at the later phases (Cohen’s d = 0.41 for Days 5–7, Cohen’s d = 0.54 for Days 8–10, Cohen’s d = 0.66 for Days 11–13). We further tested if the RMSSD reduction rate was greater in the Sham treatment group with a linear regression model: RMSSD change ~Day * Treatment. The results show that the RMSSD reduced slower in the taVNS treatment group when compared to the sham treatment, but this trend did not reach significance (coefficient of taVNS * Day interaction effect = 2.00, p = 0.21, Figure 3—figure supplement 1).
RMSSD and SDNN are two of the most commonly used methods for quantifying HRV. Bartlett’s test indicated that there are significant correlations among these measures (p < 0.01, Figure 3C). To analyze the effect of taVNS treatment on the autonomic system, we used factor analysis to identify the underlying factors. We focused on the two factors with the greatest eigenvalue. Figure 3D shows the factor loading, that is, the variance explained by HRV metrics on the two factors. The first factor correlates positively with metrics representing variability, including RMSSD, SDNN, pNNI_50, and total power, and therefore has been termed overall HRV. The second factor correlated positively with RMSSD and normalized high-frequency power, representing parasympathetic activity, and negatively correlated with the cardiac sympathetic index (CSI). Hence, it is termed parasympathetic activity. Overall HRV change from Day 1 was significantly higher in the VNS group (Figure 3E, Mann–Whitney U test, N(taVNS) = 94, N(Sham) = 95, p-value = 0.04, Cohen’s d = 0.37). The effect size was trivial between Days 2 and 4 and increased over the course of treatment. The parasympathetic activity was also significantly higher in the VNS treatment group, and we observed the largest effect size between Days 2 and 4 (Cohen’s d = 0.50, Figure 3F).
We also investigated the potential association between clinical outcomes, as measured by changes in the mRS from admission to discharge, and HRV metrics. We found that heart rate was lower in patients with improved mRS (i.e., <0) (Mann–Whitney U test, N(mRS < 0) = 122, N(mRS > 0) = 98, p-value <0.01, Cohen’s d = −0.54). Parasympathetic activity, overall HRV, and corrected QT interval did not differ significantly between patients with improved mRS and patients with worsened mRS (Figure 3—figure supplement 4).
Effects of repetitive taVNS on vascular function
Elevated blood pressure is a common occurrence in SAH and is linked with a higher risk of re-rupture of cerebral aneurysms and vasospasm (Minhas et al., 2022; Hosmann et al., 2020). In this study, patients in both treatment groups received medical treatment determined by the medical team, including vasopressors and medication for blood pressure management. We investigated whether taVNS induced any additional blood pressure changes beyond those managed by the medical team. We found that the median and mean blood pressure change from the first hospitalized day were greater than 0 for both treatment groups (Figure 4B). No significant differences were detected in changes in blood pressure and ICP between the treatment groups (Figure 4B, C). Equivalence testing confirmed that the ICP changes from the first hospitalization day were not significantly different between treatment groups, with a 2-mmHg equivalence margin (two one-sided t-tests, p[lower threshold] = 3.66 × 10−13, t[lower threshold] = 8.07; p[upper threshold] = 3.33 × 10−10, t[upper threshold] = −6.73). Equivalence testing also indicated that there were no significant different changes in blood pressure between treatment groups (two one-sided t-tests, p[lower threshold] = 0.07, t[lower threshold] = 1.51; p[upper threshold] = 0.002, t[upper threshold] = −3.00). We further verified that there were no significant changes in arterial line blood pressure obtained via continuous invasive monitoring between treatment groups (Figure 4—figure supplement 1). We subsequently compared the plethysmography peripheral perfusion index (PPI) between the treatment groups as it is a proxy metric for cardiac stroke volume and vascular tone (Elshal et al., 2021; Coutrot et al., 2021). We found that PPI change was significantly lower (Mann–Whitney U test, N(taVNS) = 83, N(Sham) = 95, Bonferroni corrected p < 0.01, Cohen’s d = −0.49). In addition, respiration rate change was significantly higher (Mann–Whitney U test, N(taVNS) = 94, N(Sham) = 95, Bonferroni corrected p = 0.02, Cohen’s d = 0.37) in the taVNS group, as compared to the Sham group (Figure 4D, E). We hypothesized that the increase in respiratory rate was a compensatory mechanism to ensure similar oxygen delivery. We found a significant negative correlation between changes in PPI and changes in respiration rate only for the taVNS treatment group (Pearson correlation coefficient = −0.37, p < 0.001, t-test, Figure 4—figure supplement 1D). The Pearson correlation coefficient for the sham treatment group is −0.08 (p = 0.36).
Acute effects of taVNS on cardiovascular function
Assessing the acute effect of taVNS on cardiovascular is crucial for its safe translation into clinical practice. We compared the acute change of heart rate, corrected QT interval, and HRV between treatment groups, as abrupt changes in the pacing cycle may increase the risk of arrhythmias (Zaniboni, 2024). The change in heart rate from treatment onset is shown in Figure 5B. We subsequently tested whether taVNS affects changes in heart rate between post- and pre-treatment. We found that the changes in heart rate were not significantly different between treatment groups although heart rate increased in the taVNS group (Wilcoxon rank-sum test, N = 188, Bonferroni corrected p = 0.03, Cohen’s d = 0.11) but not in the Sham group (Wilcoxon signed-rank test, N = 199, Bonferroni corrected p = 0.72, Cohen’s d = 0.00) (Figure 5C). However, the increase in heart rate after taVNS was within 0.5 standard deviations of daily heart rate. There were no significant differences in changes in corrected QT interval or HRV, as measured by RMSSD, SDNN, and relative power of high-frequency band between treatment groups (Figure 5D, E, Figure 5—figure supplements 1 and 3). Figure 2—figure supplement 1C shows the acute changes in uncorrected QT interval. Appendix 3—table 1 summarizes the absolute changes in cardiovascular metrics for the treatment groups. We further asked whether heart rate can serve as a biomarker that indicates which SAH patients would receive the greatest benefit from continuing taVNS treatment. We investigated the relationship between changes in heart rate from pre- to post-taVNS treatment and changes in mRS between admission and discharge using a linear mixed-effects model. In this model, the treatment group, mRS change, and their interaction were included as fixed effects, while subject was included as a random effect. Our analysis revealed that the slope between changes in heart rate and changes in mRS is significantly more negative for the taVNS treatment group (Table 2). This finding suggests that an increase in heart rate following acute taVNS treatment is associated with improved clinical outcomes (Figure 5F). Post hoc analysis showed that patients in the taVNS treatment group who had an improvement in mRS of –2 or greater compared to other patients had significantly greater increases in heart rate (Mann–Whitney U test, p = 0.02, N(mRS change < −1) = 53, N(mRS change ≥ −1) = 135, Cohen’s d = 0.34, Figure 5—figure supplement 1C). Conversely, HRV change, represented by RMSSD, was not significantly different based on mRS in SAH patients (Figure 5—figure supplement 1D).
Subsequently, we compared changes in blood pressure, PPI, ICP, and respiration rate from pre- to post-treatment periods between treatment groups. We found that changes in PPI and blood pressure were significantly higher in the taVNS group, as compared to the Sham group (Mann–Whitney U test, blood pressure: p = 0.03, Cohen’s d = 0.22, N = 180 for Sham and 159 for taVNS; PPI: p < 0.01, Cohen’s d = 0.19, N = 227 for Sham and 186 for taVNS, Figure 5—figure supplement 2). Only PPI remained significantly different between treatment groups after Bonferroni correction. The acute changes in PPI and blood pressure remained within the daily standard deviation. No significant differences in post-treatment changes in ICP or respiration rate were observed between treatment groups.
Discussion
This study examined the effects of taVNS on cardiovascular function in patients with SAH. We investigated both the cumulative and acute impacts of taVNS. The findings in our study indicate that repetitive taVNS is not associated with previously suggested risks, such as bradycardia and QT prolongation. Furthermore, repetitive taVNS treatment increased overall HRV and parasympathetic activity, which are indicators of a healthy cardiovascular system. When looking at the acute effects, taVNS only significantly increased the PPI but not heart rate, HRV, corrected QT interval, blood pressure, or ICP. The findings are summarized in Table 3. Interestingly, we found that heart rate can serve as a biomarker for identifying SAH patients who are most likely to benefit from taVNS treatment. Collectively, this study substantiates the safety of treating SAH patients with taVNS and provides foundational data for future efforts to optimize and translate taVNS therapy toward clinical use.
taVNS and autonomic system
The ANS, comprising the sympathetic and the parasympathetic nervous system, plays a critical role in maintaining physiological homeostasis. These two systems work synergistically to mediate interactions between the nervous and immune systems, which is thought to be the underlying mechanism for the immunomodulatory effect of taVNS. Our study is aligned with the finding that the autonomic balance is shifted toward sympathetic dominance following SAH (Figure 3, Figure 3—figure supplement 2; Chiu et al., 2012; Bai et al., 2023). In addition, we found that dysregulation of sympathovagal balance toward sympathetic dominance could be restored by taVNS treatment.
A key metric that reflects this restored sympathovagal balance is HRV (Figure 3F). Specifically, factor analysis based on HRV metrics showed that the parasympathetic activity was significantly higher in the taVNS treatment group. This difference was most pronounced during the early phase, between Days 2 and 4 following SAH. In addition to analyzing the correlation between the parasympathetic activity factor and established HRV measures that reflect parasympathetic activity, such as RMSSD and pNNI_50 (Figure 3C), we also examined changes in a frequency-domain HRV measure—the relative power of the high-frequency band (0.15–0.4 Hz)—to validate the accuracy of the factor analysis. The relative power of the high-frequency band is widely used to indicate respiratory sinus arrhythmia, a process primarily driven by the parasympathetic nervous system. We found that both the change in parasympathetic activity factor and relative high-frequency power were higher in the taVNS group at the early phase (Days 2–4, Figure 3—figure supplement 2). Conversely, we observed higher high-frequency power in the Sham group during the later phase. If the factor analysis successfully isolates the parasympathetic activity, there should be other factors than the parasympathetic activity affecting the relative power of the high-frequency band. One such factor is the respiration rate. The high-frequency range is between 0.15 and 0.4 Hz, corresponding to respiration’s frequency range of approximately 9–24 breaths per minute. If the respiration rate increases and exceeds 24 breaths per minute, the respiratory-driven HRV might occur at a frequency higher than the typical high-frequency band. Given that the respiration rate was higher in the taVNS treatment group, a compensatory mechanism to ensure oxygen delivery (Figure 4E), we hypothesized that the observed lower high-frequency power in the taVNS treatment group compared to sham at later phases was a result of increased respiration rate. Indeed, we found the normalized high-frequency power was higher when RR was less than 25 bpm compared to when RR >25 bpm (Cohen’s d = 0.85, Figure 3—figure supplement 3A). Moreover, an increase in RR in the taVNS treatment group was associated with a decrease in high-frequency power (Figure 3—figure supplement 3B). These control analyses underscored the necessity of performing factor analysis to robustly measure parasympathetic activities and confirm that taVNS treatment mitigated the sympathetic overactivation during the early phase.
Additionally, taVNS led to a decreased QTc without a significant change in heart rate, mimicking the effects observed with propranolol administration, a beta-blocker that reduces sympathetic activity (Solti et al., 1989). This finding suggests that repetitive taVNS reduces sympathetic overactivation and influences ventricular repolarization processes. Age affects sinus node function and is potentially associated with a higher risk of poor outcomes. To control for individual differences, including those due to age, our study compared the change in cardiovascular parameters from Day 1 within each subject across treatment groups. To further verify if age influences autonomic changes following SAH, we performed ANCOVA on autonomic function parameters with age included as a covariate. This analysis showed that age was negatively correlated with changes in heart rate, SDNN, and RMSSD from Day 1 but not with changes in QT intervals. After adjusting for age, we found that RMSSD and SDNN changes were significantly higher, while QTc changes were significantly lower in the taVNS treatment group (Appendix 4—table 1). These results align with the conclusion that repetitive taVNS treatment increased HRV and was unlikely to cause bradycardia or QT prolongation. In addition, autonomic changes following SAH may be influenced by age. Specifically, lower RMSSD and SDNN in older individuals suggest a greater shift toward sympathetic predominance following SAH (Appendix 4—table 1).
PPI is primarily influenced by cardiac output and vascular tone. Elevated PPI is associated with vasodilation and/or increased stroke volume. In the Sham group, increases in both PPI and blood pressure were observed when compared to Day 1 values (Figure 3). This effect may be due to higher stroke volume resulting from sympathetic activation following SAH. Alternatively, this could represent the heightened need for vasopressor interventions to improve cerebral perfusion due to more robust sympathetically driven cerebral vasospasm. The increase in PPI was less for the taVNS treatment group (Figure 4), suggesting a restored autonomic balance in the taVNS treatment group. However, the effects of taVNS on blood pressure require further investigation as more than half of the patients were on vasopressor and ionotropic drugs. Intuitively, sympathetic activation is associated with increases in both PPI and blood pressure. The blood pressure management might lead to similar blood pressure changes between the two treatment groups. Although repetitive taVNS increases HRV days after initiation of the treatment, this effect is not seen acutely. Also, while repetitive taVNS was associated with a reduced PPI and no change in heart rate and blood pressure, there were small acute increases in PPI, heart rate, and blood pressure. All patients who were capable of verbal communication were asked if they felt any prickling or pain during all sessions. We confirmed that the current taVNS protocol is below the perception threshold for all trialed patients. Altogether, successful activation of the afferent vagal pathway by taVNS increased arousal, resulting in increased heart rate (Naredi et al., 2000; Skora et al., 2024Sharon et al., 2021; Skora et al., 2024; Tan et al., 2024). These speculative mechanisms warrant further validation through animal or pharmacological studies directly investigating the effects of taVNS on autonomic function and vascular tone.
Considerations for applying taVNS on SAH patients
Blood pressure management and cardiac function monitoring are crucial in patients following SAH (Minhas et al., 2022). This study shows that blood pressure, QT interval, and heart rate over days were not significantly different between taVNS and sham treatment groups. This suggests that adding taVNS in treatments for SAH patients is unlikely to cause adverse blood pressure alterations or cardiac complications. Our findings suggest that repetitive taVNS could enhance parasympathetic tone following SAH. This effect could lead to favorable clinical outcomes as lower HRV was found to be associated with neurocardiogenic injury (Bai et al., 2023; Megjhani et al., 2020). Given the negative association between pro-inflammatory markers and HRV, our finding that HRV was higher in the taVNS treatment group aligns with the findings of primary outcomes of this clinical trial, which showed that taVNS treatment reduced pro-inflammatory cytokines, including TNF-α and interleukin-6 (Huguenard et al., 2024b; Williams et al., 2019). The consistency between these findings strengthens the evidence supporting the anti-inflammatory effects of taVNS. In addition, the sympathetic predominance following SAH is implicated in an increased risk of delayed cerebral vasospasm, which is most commonly detected 5–7 days after SAH (Budohoski et al., 2013). Given that taVNS treatment mitigated the sympathetic overactivation before the typical onset of cerebral vasospasm, it could potentially reduce the severity of this complication. Additionally, reduced PPI was associated with increased respiration rate only for the taVNS treatment group, suggesting that the autonomic system self-regulates to maintain cardiovascular homeostasis. Thus, it is important to consider autonomic system self-regulation when studying the therapeutic effects of taVNS (Bjerkne Wenneberg et al., 2020). Also, while acute cardiovascular changes were noted after taVNS, these changes were within normal daily variations in this study, making them unlikely to pose a risk to the patient. That said, the observed acute increases in PPI following taVNS necessitate caution when considering taVNS treatment for patients to whom peripheral vasodilatation is not desired.
As we pioneer the application of taVNS as an immunomodulation technique in SAH patients, we adopt parameters (20 Hz, 0.4 mA) reported in similar studies (Jelinek et al., 2024). The current study provides a basis for future preclinical and clinical studies of taVNS in this patient population. To build on our findings, a systematic evaluation of the relationship between parameters such as frequency, intensity, and duration and taVNS’s effects on the immune system and cardiovascular function is necessary to establish taVNS as an effective therapeutic option for SAH patients (Dusi et al., 2022).
Limitations and outlook
While this study supports the safety of taVNS treatment in SAH patients, we should be cautious when generalizing these findings to broader clinical populations. The current study did not explore the effects of taVNS on less commonly used cardiovascular metrics, such as QTc dispersion. Our study considers each day as an independent sample for the following considerations: (1) heart rate and HRV metrics exhibited great daily variations. Their value on 1 day was not predictive of the metrics on another day, which could be due to medications, interventions, or individualized SAH recovery process during the patient’s stay in the ICU. (2) SAH patients in the ICU often experience daily changes in clinical status, including fluctuations in ICP, blood pressure, neurological status, and other vital signs. (3) Day-to-day cardiovascular function changes varied as the patient recovered or encountered setbacks. To conclusively establish that there is no significant cardiovascular effect of repetitive taVNS on any given day following SAH, we would need to perform statistical tests between treatment groups for each day. In this context, 64 subjects per treatment group are required to achieve 80% power, assuming a medium effect size (Cohen’s d = 0.5)and 0.05 type I error probability (two-sample t-test).
Mild cardiac abnormalities are common in SAH patients (Norberg et al., 2018), complicating the precise calculation of cardiovascular metrics from ECG signals and the interpretation of the results. Systematic verification of methods for calculating cardiovascular metrics to ensure their applicability in SAH patients is crucial. We noticed a high variance of change in heart rate for Days 5–7, 8–10, and 11–13 for both treatment groups (Figure 2D). This may be due to the small sample size in the later days, given that the mean duration of hospitalization for the 24 subjects included in this study was 11.3 days with a standard deviation of 6.4. Differences in medical history and clinical outcomes during hospitalization may also explain the variance of change in heart rate for the later days. For example. heart rate was lower in patients with improved mRS scores (Figure 3—figure supplement 4B). Understanding the association between cardiovascular metrics and clinical assessments, such as vasospasm and inflammation, could help decide whether future taVNS trials should control for these factors when evaluating the effects of taVNS on cardiovascular function. Additional care should be paid when interpreting the results of blood pressure, as hypertension was intentionally induced for some patients being treated for vasospasm. Patient medical histories are summarized in Table 1.
Kulkarni et al. showed that the response to low-level tragus stimulation varied among patients with atrial fibrillation (Kulkarni et al., 2021). Similarly, in our study, not all patients in the taVNS treatment group showed a reduction in mRS scores (improved degree of disability or dependence). This differential response may be inherent to taVNS and potentially influenced by factors such as anatomical variations in the distribution of the vagus nerve at the outer ear. These findings underscore the importance of using acute biomarkers to guide patient selection and optimize stimulation parameters. Furthermore, we found that increased heart rate was a potential acute biomarker for identifying SAH patients who are most likely to respond favorably to taVNS treatment. Translating this finding into clinical practice will require further research to elucidate the mechanisms by which an acute increase in heart rate may predict the outcomes of patients receiving taVNS, including its relationship with neurological evaluations, vasospasm, echocardiography, and inflammatory markers.
Conclusions
Utilizing taVNS as a neuromodulation technique in SAH patients is safe without inducing bradycardia or QT prolongation. Repetitive taVNS treatment increased parasympathetic activity. Acute taVNS elevated heart rate, which might be an acute biomarker to identify SAH patients who are likely to respond favorably to taVNS treatment.
Methods
Study participants
Participants in this study were recruited from adult patients who were admitted to the ICU at Barnes Jewish Hospital, St. Louis, MO, following an acute, spontaneous, aneurysmal SAH. Inclusion criteria were: (1) Patients with SAH confirmed by CT scan; (2) Age >18; (3) Patients or their legally authorized representative are able to give consent. Exclusion criteria were: (1) Age <18; (2) Use of immunosuppressive medications; (3) Receiving ongoing cancer therapy; (4) Implanted electrical device; (5) Sustained bradycardia on admission with a heart rate <50 beats per minute for >5 min; (6) Considered moribund/at risk of imminent death. Participants were randomized to receive either the taVNS (N = 11) or Sham (N = 13) treatment. Patients were enrolled prior to randomization by a member of the research team who went through the informed consent process with the patient or their legally authorized representative. Treatment group assignment was via a computer-generated randomization sequence, with the next number obscured until patient enrollment. Research team members who applied the ear clips and set stimulation parameters were not blinded to the treatment. The participants, the medical team who dictated all management decisions for the patient’s SAH, and the outcomes assessors who assigned mRS at admission and discharge were blinded to the treatment. The structure of this study is shown in Figure 1B. This study was approved by the Washington University School of Medicine Review Board and was conducted in accordance with institutional and national ethics guidelines and the Declaration of Helsinki (Clinical trial number: NCT04557618).
taVNS protocol
Following randomization, enrolled patients underwent 20 min of either taVNS or sham stimulation twice daily during their stay in the ICU. During treatment periods, a portable TENS device (TENS 7000 Digital TENS Unit, Compass Health Brands, OH, USA) was connected to the patient’s left ear using two ear clips (Figure 1C, D). For taVNS treatments, these ear clips were placed along the concha of the ear, while for sham treatments, the clips were placed along the earlobe to avoid stimulation of the auricular vagus nerve from tactile pressure (Figure 1C). For the taVNS group, stimulation parameters were selected based on values reported in prior studies that sought to maximize vagus somatosensory evoked potentials while avoiding the perception of pain: 20 Hz frequency, 250 µs pulse width, and 0.4 mA intensity (de Gurtubay et al., 2021). The stimulation was not perceptible for the patients. No electrical current was delivered during sham treatments. For both groups, the TENS device was connected to the patient and a bedside recording computer. The computer recorded continuous ECG and vital signs, including blood pressure, temperature, respiration rate, PPI, ICP, and arterial blood pressure. The collection of ICP and arterial blood pressure data varied, being dependent on the treatment protocol assigned by the clinical team, and thus was not uniformly available for all patients throughout the study. Please see Huguenard et al., 2024a for a detailed protocol of this study.
Data processing
A 3-lead system was used for ECG. ECG signals, sampled at 500 Hz, and other vital signs, such as blood pressure, sampled at 1 Hz, were recorded from the Intellivue patient monitor (Philips, Netherlands) using vitalDB software (Lee and Jung, 2018).
To calculate cardiac metrics, we first applied a 0.5-Hz fifth-order high-pass Butterworth filter and a 60-Hz powerline filter on ECG data to reduce artifacts (Makowski et al., 2021). We detected QRS complexes based on the steepness of the absolute gradient of the ECG signal using the Neurokit2 software package (Makowski et al., 2021). R-peaks were detected as local maxima in the QRS complexes. P waves, T waves, and QRS complexes were delineated based on the wavelet transform of the ECG signals proposed by Figure 2A–C; Martínez et al., 2004. This algorithm identifies the QRS complex by searching for modulus maxima, which are peaks in the wavelet transform coefficients that exceed specific thresholds. The onset of the QRS complex is determined as the beginning of the first modulus maximum before the modulus maximum pair created by the R wave. To identify the T wave, the algorithm searches for local maxima in the absolute wavelet transform in a search window defined relative to the QRS complex. Thresholding is used to identify the offset of the T wave. RR intervals were preprocessed to exclude outliers, defined as RR intervals greater than 2 s or less than 300 ms. RR intervals with >20% relative difference to the previous interval were considered ectopic beats and excluded from analyses. After preprocessing, RR intervals were used to calculate heart rate, HRV, and corrected QT (QTc) based on Bazett’s formula: (Bazett, 1920). The corrected QT interval (QTc) estimates the QT interval at a standard heart rate of 60 bpm. HRV measures included the root mean square of successive difference of normal RR intervals (RMSSD), indicating parasympathetic activity, and the standard deviation of normal RR intervals (SDNN), a clinical measure of cardiac risk (Kleiger et al., 2005; Shaffer and Ginsberg, 2017). HRV calculations are detailed in Appendix 1.
To investigate the effect of repetitive taVNS on cardiovascular function, we compared HRV, heart rate, corrected QT intervals, blood pressure, and ICP calculated over 24 hr between patients receiving taVNS and sham treatment. In addition, we compared the mean PPI and respiration rate over 24 hr between treatment groups to determine the effects of repetitive taVNS on the autonomic system. Data collection commenced on the first day of each patient’s ICU admission. The average duration of continuous data recording was 11.1 days, with a standard deviation of 6.8 days. To analyze the effects of taVNS treatment more granularly, we segmented the changes in these metrics from the initial day at 3-day intervals, facilitating comparison between the taVNS and sham treatment groups over the course of their ICU stay.
To study the effects of acute treatment over time, we focused on blood pressure, HRV, heart rate, and corrected QT intervals 20 min before treatment (pre-treatment), during the 20 min treatment (during-treatment), and 20 min after treatment (post-treatment). The treatment event signals were rectified and binarized based on their half-maximum value to identify the treatment onset and offset (Figure 2A). We calculated metrics using 6 min sliding windows over ECG data starting from treatment onset/offset and moving bi-directionally with a 3-min step. To correct daily and between-subject variation, we applied the same sliding window strategy to calculate the mean and standard deviation of these cardiac metrics for each patient each day as a reference. Subsequently, HRV, heart rate, and corrected QT interval around treatment onset/offset were normalized based on the reference. In addition, we calculated the difference in blood pressure, HRV, and heart rate, and corrected QT intervals between during- and pre-treatment, as well as the difference between post- and pre-treatment for each patient and for each treatment. To study the effects of acute taVNS, we compared the two differences between the treatment groups.
Factor analysis
We performed an exploratory factor analysis to identify the factors underlying autonomic system activity. Besides RMSSD and SDNN, variables derived from preprocessed RR intervals and used to perform factor analysis included the percentage of successive normal-to-normal (NN) Intervals that differ by more than 50 ms (pNNI_50), total power (below 0.4 Hz), normalized high-frequency power (0.15–0.4 Hz), cardiac vagal index, and CSI. The total power is thought to represent the overall HRV, while normalized high-frequency power primarily reflects parasympathetic activity (Kleiger et al., 2005). These variables were normalized using a z-score method based on individual daily means and standard deviations before factor analysis. Factor analysis was performed using the factor_analyzer Python package. The number of factors was set to 2 based on the Scree plot. The factor loadings were calculated using the Minimum Residual Method. After factor extraction, a Varimax rotation was applied for better interpretability so that each factor had high loadings for a smaller number of variables and low loadings for the remaining variables.
Statistical analyses
To investigate the effect of taVNS at the phase of early brain injury and later phases, we grouped the change of HRV, heart rate, and corrected QT interval from the first hospitalized day in bins of 3 days. The change in blood pressure, ICP, respiration rate, and PPI from the first hospitalization day were also compared between treatment groups. We used t-tests for comparisons between treatment groups when the data were normally distributed, as determined by the Shapiro–Wilk test. We employed Mann–Whitney U tests for non-normally distributed data. We used Wilcoxon signed-rank tests to compare the difference in heart rate between post- and during-treatment against 0. To control the familywise error rate, we applied Bonferroni correction. Specifically, when investigating the cardiac effects of taVNS, we compared six metrics between treatment groups, including heart rate, corrected QT interval, RMSSD, SDNN, and two factors representing HRV. Consequently, the p-values were corrected by a factor of six. In this study, we reported the statistical power achieved for tests that yielded non-significant results. The achieved power is calculated based on a two-sample t-test assuming a medium effect size (Cohen’s d of 0.5) and a type I error probability (a) of 0.05. We used two one-sided tests to confirm that taVNS did not induce long-term changes in heart rate, corrected QT interval, or blood pressure, with equivalency test margins set to 5 bpm for heart rate, 50 ms for QT interval, and 2 mmHg for blood pressure. A summary of statistical tests is provided in Appendix 2—table 1.
Appendix 1
Calculation of cardiovascular metrics
1. Heart rate (HR) over a period is the mean of over this period.
2. The percentage of prolonged QT intervals is defined as the number of intervals with a corrected QT interval (QTc) ≥500 ms, divided by the total number of QT intervals.
3. Root mean square of successive difference of normal RR intervals (RMSSD) for a given period:
, where N is the number of NN intervals within this period.
4. Standard deviation of normal RR intervals (SDNN) for a given period:
, where represents the mean of NN intervals within this period.
5. To calculate total power and normalized high-frequency power, power spectral density is calculated from NN intervals with the help of Python package hrv-analysis. Linear interpolation is used to make NN intervals evenly spaced, assuming a sampling frequency of 4 Hz. Welch’s method is used to calculate power spectral density .
6. Cardiac vagal index (CVI) and cardiac sympathetic index (CSI):
The width (SD1) and length (SD2) of the Poincaré plot, a graphical representation of the correlation between consecutive NN intervals, are used to calculate CVI and CSI.
Appendix 2
Appendix 3
Appendix 4
Data availability
All relevant data has been made publicly available on Dryad at https://doi.org/10.5061/dryad.rfj6q57kw. Codes for generating the figures are available at https://github.com/GanshengT/taVNS_SAH (copy archived at Tan, 2024).
-
Dryad Digital RepositoryElectrocardiography and vital signs recorded in patients with subarachnoid hemorrhage (NCT04557618).https://doi.org/10.5061/dryad.rfj6q57kw
References
-
An analysis of the time-relations of electrocardiogramsHeart (British Cardiac Society) 7:353–370.
-
Heart rate variability monitoring for the detection of delayed cerebral ischemia after aneurysmal subarachnoid hemorrhageActa Anaesthesiologica Scandinavica 64:945–952.https://doi.org/10.1111/aas.13582
-
Chronic vagus nerve stimulation in Crohn’s disease: a 6‐month follow‐up pilot studyNeurogastroenterology & Motility 28:948–953.https://doi.org/10.1111/nmo.12792
-
Clinical relevance of cerebral autoregulation following subarachnoid haemorrhageNature Reviews. Neurology 9:152–163.https://doi.org/10.1038/nrneurol.2013.11
-
Depressed sympathovagal balance predicts mortality in patients with subarachnoid hemorrhageThe American Journal of Emergency Medicine 30:651–656.https://doi.org/10.1016/j.ajem.2011.02.037
-
Perfusion index: Physical principles, physiological meanings and clinical implications in anaesthesia and critical careAnaesthesia Critical Care & Pain Medicine 40:100964.https://doi.org/10.1016/j.accpm.2021.100964
-
Aneurysmal Subarachnoid HemorrhageJournal of Neurosurgical Anesthesiology 27:222–240.https://doi.org/10.1097/ANA.0000000000000130
-
Neuromodulation devices for heart failureEuropean Heart Journal Supplements 24:E12–E27.https://doi.org/10.1093/eurheartjsupp/suac036
-
Plethysmographic peripheral perfusion index: could it be a new vital sign?Frontiers in Medicine 8:651909.https://doi.org/10.3389/fmed.2021.651909
-
Endogenous arterial blood pressure increase after aneurysmal subarachnoid hemorrhageClinical Neurology and Neurosurgery 190:105639.https://doi.org/10.1016/j.clineuro.2019.105639
-
Heart rate variability: measurement and clinical utilityAnnals of Noninvasive Electrocardiology 10:88–101.https://doi.org/10.1111/j.1542-474X.2005.10101.x
-
Low-level tragus stimulation modulates atrial alternans and fibrillation burden in patients with paroxysmal atrial fibrillationJournal of the American Heart Association 10:e020865.https://doi.org/10.1161/JAHA.120.020865
-
NeuroKit2: A Python toolbox for neurophysiological signal processingBehavior Research Methods 53:1689–1696.https://doi.org/10.3758/s13428-020-01516-y
-
A wavelet-based ECG delineator: evaluation on standard databasesIEEE Transactions on Bio-Medical Engineering 51:570–581.https://doi.org/10.1109/TBME.2003.821031
-
The vagus nerve and the inflammatory reflex--linking immunity and metabolismNature Reviews. Endocrinology 8:743–754.https://doi.org/10.1038/nrendo.2012.189
-
Subarachnoid hemorrhage and inflammation: bench to bedside and backSeminars in Neurology 25:435–444.https://doi.org/10.1055/s-2005-923537
-
Inflammation in subarachnoid hemorrhage and delayed deterioration associated with vasospasm: A reviewActa Neurochirurgica. Supplement 115:233–238.https://doi.org/10.1007/978-3-7091-1192-5_42
-
Inflammatory events following subarachnoid hemorrhage (SAH)Current Neuropharmacology 16:1385–1395.https://doi.org/10.2174/1570159X16666180412110919
-
An overview of heart rate variability metrics and normsFrontiers in Public Health 5:258.https://doi.org/10.3389/fpubh.2017.00258
-
Transcutaneous vagus nerve stimulation in humans induces pupil dilation and attenuates alpha oscillationsThe Journal of Neuroscience 41:320–330.https://doi.org/10.1523/JNEUROSCI.1361-20.2020
-
The effect of sympathetic and parasympathetic activity on QT duration. Clinical study in patients with normal and prolonged QT timeCor et Vasa 31:9–15.
-
Control of inflammation using non-invasive neuromodulation: past, present and promiseInternational Immunology 34:119–128.https://doi.org/10.1093/intimm/dxab073
-
Unruptured intracranial aneurysms: a reviewJournal of Neurosurgery 96:3–42.https://doi.org/10.3171/jns.2002.96.1.0003
-
Heart rate variability and inflammation: A meta-analysis of human studiesBrain, Behavior, and Immunity 80:219–226.https://doi.org/10.1016/j.bbi.2019.03.009
-
Electrocardiographic abnormalities predict adverse clinical outcomes in patients with subarachnoid hemorrhageJournal of Stroke and Cerebrovascular Diseases 25:2653–2659.https://doi.org/10.1016/j.jstrokecerebrovasdis.2016.07.011
Article and author information
Author details
Funding
American Association of Neurological Surgeons
- Anna L Huguenard
Aneurysm and AVM Foundation
- Anna L Huguenard
National Institutes of Health (R01-EB026439)
- Peter Brunner
National Institutes of Health (P41-EB018783)
- Peter Brunner
National Institutes of Health (U24-NS109103)
- Peter Brunner
National Institutes of Health (R21-NS128307)
- Eric C Leuthardt
McDonnell Center for Systems Neuroscience
- Peter Brunner
- Eric C Leuthardt
Fondazione Neurone
- Peter Brunner
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors acknowledge physicians/nurses for helping administer treatment. The authors thank Dr. Paul Cassidy for his contributions to the scientific editing of this manuscript, supported by the Institute of Clinical and Translational Sciences grant UL1TR002345 from the National Center for Advancing Translational Sciences (NCATS).
Ethics
Clinical trial registration NCT04557618.
This study was approved by the Washington University School of Medicine Review Board and was conducted in accordance with institutional and national ethics guidelines and the Declaration of Helsinki (Clinical trial number: NCT04557618). Participants or their legally authorized representatives provided written informed consent prior to enrollment in the study.
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.100088. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2024, Tan, Huguenard et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 259
- views
-
- 16
- downloads
-
- 0
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Immunology and Inflammation
- Medicine
Background:
Individuals with Down syndrome (DS), the genetic condition caused by trisomy 21 (T21), display clear signs of immune dysregulation, including high rates of autoimmunity and severe complications from infections. Although it is well established that T21 causes increased interferon responses and JAK/STAT signaling, elevated autoantibodies, global immune remodeling, and hypercytokinemia, the interplay between these processes, the clinical manifestations of DS, and potential therapeutic interventions remain ill defined.
Methods:
We report a comprehensive analysis of immune dysregulation at the clinical, cellular, and molecular level in hundreds of individuals with DS, including autoantibody profiling, cytokine analysis, and deep immune mapping. We also report the interim analysis of a Phase II clinical trial investigating the safety and efficacy of the JAK inhibitor tofacitinib through multiple clinical and molecular endpoints.
Results:
We demonstrate multi-organ autoimmunity of pediatric onset concurrent with unexpected autoantibody-phenotype associations in DS. Importantly, constitutive immune remodeling and hypercytokinemia occur from an early age prior to autoimmune diagnoses or autoantibody production. Analysis of the first 10 participants to complete 16 weeks of tofacitinib treatment shows a good safety profile and no serious adverse events. Treatment reduced skin pathology in alopecia areata, psoriasis, and atopic dermatitis, while decreasing interferon scores, cytokine scores, and levels of pathogenic autoantibodies without overt immune suppression.
Conclusions:
JAK inhibition is a valid strategy to treat autoimmune conditions in DS. Additional research is needed to define the effects of JAK inhibition on the broader developmental and clinical hallmarks of DS.
Funding:
NIAMS, Global Down Syndrome Foundation.
Clinical trial number:
-
- Immunology and Inflammation
- Medicine
Metabolic abnormalities associated with liver disease have a significant impact on the risk and prognosis of cholecystitis. However, the underlying mechanism remains to be elucidated. Here, we investigated this issue using Wilson’s disease (WD) as a model, which is a genetic disorder characterized by impaired mitochondrial function and copper metabolism. Our retrospective clinical study found that WD patients have a significantly higher incidence of cholecystitis and a poorer prognosis. The hepatic immune cell landscape using single-cell RNA sequencing showed that the tissue immune microenvironment is altered in WD, mainly a major change in the constitution and function of the innate immune system. Exhaustion of natural killer (NK) cells is the fundamental factor, supported by the upregulated expression of inhibitory receptors and the downregulated expression of cytotoxic molecules, which was verified in clinical samples. Further bioinformatic analysis confirmed a positive correlation between NK cell exhaustion and poor prognosis in cholecystitis and other inflammatory diseases. The study demonstrated dysfunction of liver immune cells triggered by specific metabolic abnormalities in WD, with a focus on the correlation between NK cell exhaustion and poor healing of cholecystitis, providing new insights into the improvement of inflammatory diseases by assessing immune cell function.