# Evaluation of Cardiovascular Risk in Hidradenitis Suppurativa Patients Using Heart Rate Variability (HRV) Analysis

Authors: Nevena Skroza; Alessandra Mambrin; Ilaria Proietti; Veronica Balduzzi; Nicoletta Bernardini; Anna Marchesiello; Simone Michelini; Ersilia Tolino; Salvatore Volpe; Patrizia Maddalena; Antonella Ciaramella; Gianfranco Raimondi; Concetta Potenza
Journal: Cardiovascular Therapeutics
 (2020)

---

## Abstract

Hidradenitis suppurativa (HS) is a chronic inflammatory skin disease associated with elevated prevalence of comorbidities, especially metabolic and cardiovascular diseases. We used a tool called Heart Rate Variability (HRV) in order to assess the correlation between HS and alterations of the sympathetic-vagal equilibrium in the autonomic cardiovascular regulation system. We found increased sympathetic activity, associated with a higher risk of cardiovascular disease. HS, according to our results, is an independent cardiovascular risk factor.


---

## Body

## 1. Introduction

Hidradenitis suppurativa (HS) is a chronic inflammatory skin disease characterized by recurrent skin nodules, abscesses, and draining sinus tracts affecting primarily the apocrine gland-rich regions of the skin, especially the axillae and groin [1, 2].HS has a chronic evolution with unpredictable periodic exacerbations. The lesions show a tendency to turn into fistulous tunnels, draining externally a serum-purulent material healing with the formation of retracting scars. Prevalence is estimated from <1% to 4% [3, 4], and the incidence is 11.4 cases per 100,000 patients, which is steadily increasing. Previous studies revealed a female predominance of 3 : 1 in the general population [5, 6]. Genetic factors play a fundamental role in the pathogenesis of HS, and one-third of patients report a positive family history of the disease. However, although the HS pathogenetic mechanism is still unknown, follicular hyperkeratinization probably characterizes the first step of the inflammatory process [7]. Histological examination revealed follicular hyperkeratosis and lymphocytic perifollicular infiltration followed by the disruption of the follicular unit.Furthermore, the immunohistochemistry studies revealed increased production of proinflammatory molecules, like tumor necrosis factor- (TNF-)α, interleukin- (IL-) 1b, IL-12, IL-17, and IL-23.In addition, TNF-α is one of the main cytokines that can be associated with the development of insulin resistance and hyperlipidemia, characteristics associated with metabolic syndrome.Recent works suggest that HS should be viewed as a systemic inflammatory disease linked most commonly to metabolic, gastrointestinal, rheumatologic, psychiatric, and cardiovascular (CV) comorbidities.The prevalence of the sympathetic branch of the autonomous cardiovascular system in the population reflects an increased CV risk. In this study, we used Heart Rate Variability (HRV) as a noninvasive tool to assess the sympathetic-vagal equilibrium in the autonomous CV regulation system [8–10]. Our aim was to evaluate whether autonomic function is affected in young HS patients and whether to consider HS as an independent cardiovascular risk factor [11].


## 2. Material and Methods

### 2.1. Subjects

Subjects are those affected by moderate to severe HS who attended our Outpatient Service of Dermatology from October 2018 to April 2019.The heart rate (HR), RR trait (in ms), PQ interval (v.n.: 120-220 ms), QRS (v.n.: <100 ms), and QTc (v.n.: <430 ms) were measured.Patients were informed not to assume any modifiers of autonomic nervous system functions such as psychoactive drugs and caffeine within 4 hours prior to the electrocardiogram (ECG). Informed consent was collected from each patient.


### 2.2. Assessments

We used the Cardiolab CE pocket PC ECG system in order to obtain an 8 min digital ECG at rest in supine position and during 8 minutes of orthostatic position obtained by the use of a passive tilt table.With dedicated software (Cardiolab—Xai Medica), we analyzed the HRV indexes in the time and frequency domains. We will also analyze the nonlinear dynamics of HRV by Kubios software [12].


### 2.3. Linear Methods

We used linear methods in order to obtain a direct estimation of HRV through two types of analysis: time-domain analysis and frequency-domain analysis.The time-domain method is based on NN intervals which are related to RR intervals in the ECG. They are analyzed to acquire variables such as SDNN (Standard Deviation of all NN intervals) and RMSSD (Root Mean Square of Successive Differences between adjacent NNs); both indexes represent parasympathetic tone.Frequency-domain analysis, which is based on the power spectral density of the heart rate time series, highlights the issue of the underlying rhythms of the mechanisms controlling the heart rate (HR) and identified three major spectral peaks: high frequency (HF: 0.15-0.4 Hz), low frequency (LF: 0.04-0.15 Hz), and very low frequency (VLF: below 0.04 Hz) in the adult HR spectrum.LF is an indicator of both sympathetic and parasympathetic activities. HF reflects parasympathetic activity. The LF/HF ratio is an index representing overall balance between sympathetic and parasympathetic systems. Higher values reflect dominance of the sympathetic branch, while lower ones reflect the predominance of the parasympathetic action.


### 2.4. Nonlinear Methods

The nonlinear methods (Poincarè plots and Detrended fluctuation analysis) are not modified by environmental features and measure the complex mechanisms regulating the signal [13].Poincarè plots is a two-dimensional vector analysis that was used to quantify the shape of the plots (Figure1). In this quantitative method, short-term (SD1) and long-term (SD2) RR interval variability and the ellipse area of the plot are quantified separately. SD1 is an indicator of vagal modulation of the sinus node. The HR correlations were defined separately for short-term (<11 beats, α1) and longer-term (>11 beats, α2) RR interval data. The α1 index was positively correlated with LF in normalized units. The other nonlinear method is the detrended fluctuation analysis (DFA) that consists of a procedure utilized to measure the fractal scaling assets of short-term and intermediate-term variability in the RR ECG intervals. The decrease of the α1 index is considered a mortality prognostic factor in patients affected by severe cardiovascular conditions such as ischemic events and cardiac insufficiency.Figure 1
An example of Poincarè plot and DFA during a Tilt Test (TT) in a patient with HS.


### 2.5. Data Analysis

Our study results were examined through two types of software: Cardiolab, Xai Medica (used to investigate HRV linear methods), and Kubios HRV 2.0 (for the HRV nonlinear analysis examination). We expressed our data as themean±standarderror. We calculated the t-test in order to estimate the quantitative variables, and the chi-squared test was performed to estimate the qualitative ones.We performed a statistical analysis using the software SigmaStat 3.5 (Systat Software Inc., Point Richmond, CA, USA). We considered aP value with a significance level<0.05.


## 2.1. Subjects

Subjects are those affected by moderate to severe HS who attended our Outpatient Service of Dermatology from October 2018 to April 2019.The heart rate (HR), RR trait (in ms), PQ interval (v.n.: 120-220 ms), QRS (v.n.: <100 ms), and QTc (v.n.: <430 ms) were measured.Patients were informed not to assume any modifiers of autonomic nervous system functions such as psychoactive drugs and caffeine within 4 hours prior to the electrocardiogram (ECG). Informed consent was collected from each patient.


## 2.2. Assessments

We used the Cardiolab CE pocket PC ECG system in order to obtain an 8 min digital ECG at rest in supine position and during 8 minutes of orthostatic position obtained by the use of a passive tilt table.With dedicated software (Cardiolab—Xai Medica), we analyzed the HRV indexes in the time and frequency domains. We will also analyze the nonlinear dynamics of HRV by Kubios software [12].


## 2.3. Linear Methods

We used linear methods in order to obtain a direct estimation of HRV through two types of analysis: time-domain analysis and frequency-domain analysis.The time-domain method is based on NN intervals which are related to RR intervals in the ECG. They are analyzed to acquire variables such as SDNN (Standard Deviation of all NN intervals) and RMSSD (Root Mean Square of Successive Differences between adjacent NNs); both indexes represent parasympathetic tone.Frequency-domain analysis, which is based on the power spectral density of the heart rate time series, highlights the issue of the underlying rhythms of the mechanisms controlling the heart rate (HR) and identified three major spectral peaks: high frequency (HF: 0.15-0.4 Hz), low frequency (LF: 0.04-0.15 Hz), and very low frequency (VLF: below 0.04 Hz) in the adult HR spectrum.LF is an indicator of both sympathetic and parasympathetic activities. HF reflects parasympathetic activity. The LF/HF ratio is an index representing overall balance between sympathetic and parasympathetic systems. Higher values reflect dominance of the sympathetic branch, while lower ones reflect the predominance of the parasympathetic action.


## 2.4. Nonlinear Methods

The nonlinear methods (Poincarè plots and Detrended fluctuation analysis) are not modified by environmental features and measure the complex mechanisms regulating the signal [13].Poincarè plots is a two-dimensional vector analysis that was used to quantify the shape of the plots (Figure1). In this quantitative method, short-term (SD1) and long-term (SD2) RR interval variability and the ellipse area of the plot are quantified separately. SD1 is an indicator of vagal modulation of the sinus node. The HR correlations were defined separately for short-term (<11 beats, α1) and longer-term (>11 beats, α2) RR interval data. The α1 index was positively correlated with LF in normalized units. The other nonlinear method is the detrended fluctuation analysis (DFA) that consists of a procedure utilized to measure the fractal scaling assets of short-term and intermediate-term variability in the RR ECG intervals. The decrease of the α1 index is considered a mortality prognostic factor in patients affected by severe cardiovascular conditions such as ischemic events and cardiac insufficiency.Figure 1
An example of Poincarè plot and DFA during a Tilt Test (TT) in a patient with HS.


## 2.5. Data Analysis

Our study results were examined through two types of software: Cardiolab, Xai Medica (used to investigate HRV linear methods), and Kubios HRV 2.0 (for the HRV nonlinear analysis examination). We expressed our data as themean±standarderror. We calculated the t-test in order to estimate the quantitative variables, and the chi-squared test was performed to estimate the qualitative ones.We performed a statistical analysis using the software SigmaStat 3.5 (Systat Software Inc., Point Richmond, CA, USA). We considered aP value with a significance level<0.05.


## 3. Results

Sixteen HS patients were recruited, including 9 females and 7 males aged between 16 and 48 years. Four of sixteen patients have normal weight, 6/16 were overweight, 3/16 patients presented with grade I obesity, 2/16 patients presented grade II obesity, and one patient presented with grade III obesity. In addition, eight of sixteen patients were smokers.Information about gender, age, weight, height, BMI, and smoking habits is reported in Table1.Table 1
HS patients included in the study and assessed with ECG.
Age (years)SexHeight (cm)Weight (kg)BMISmokeHR (bpm)RR (ms)PQ (ms)QRS (ms)QTc (ms)41F1.688329.4Yes7876811488324.0016F1.657226.4No9364214858296.0031M1.79532.9Yes8174014092344.0015M1.629335.4No8769014078332.0048M1.89027.8Yes60100417084410.0013F1.427939.2No8967213480303.0028F1.78529.4Yes748169870361.0016M1.858524.8No59102015288404.0036M1.667727.9No7085614888370.0019F1.767223.2No7876812276324.0032F1.658330.5Yes7975814886322.0028F1.577630.8Yes6789214484349.0024F1.713045Yes9364415688297.0016F1.676021.5No83720172108314.0022M1.859527.8Yes74806104110359.0021M701.8520.5No6986612088372.0025.386.2379.8029.5377.125791.375138.12585.375342We did not observe any significative difference when we catalogued patients according to BMI, though an increase of the heart rate (RR) was observed in subjects with aBMI>25 (Table 2).Table 2
HS patients included in the study and assessed with ECG, divided by BMI.
(a)
Normal weight
Physical conditionSmokeHRRRPQQRSQTcNormal weightNo59102015288404.00Normal weightNo7876812276324.00Normal weightNo83720172108314.00Normal weightNo6986612088372.0072.25843.5141.590353.50(b)
Overweight and obesity
Physical conditionSmokeHRRRPQQRSQTcObesity IYes8174014092344.00Obesity IYes7975814886322.00Obesity IYes67892144.084349.00Obesity IINo8769014078332.00Obesity IINo8967213480303.00Obesity IIIYes9364415688297.00OverweightYes7876811488324OverweightNo93.064214858296.00OverweightYes60100417084401.00OverweightYes748169870361.00OverweightNo7085614888370.00OverweightYes74806104110359.0078.7577413785.38338.17On the difficulty of having clean data, we analyzed 8/16 patients in detail, as preliminary data: 3 males and 5 females aged between 16 and 48 years. Of these subjects, only 2 patients presented with a normal weight while 6 patients had a BMI greater than 24.99 kg/m2 (2 patients were overweight, 3 patients had a grade I obesity, and 1 patient had a grade III obesity). Five patients were smokers, and 3 patients were nonsmokers. In all patients, the age of onset was between 10 and 28 years while the age of diagnosis was between 15 and 45 years. The ECG recorded did not show any alteration, except in one patient, where it showed a right conduction disorder (Table 3).Table 3
Characteristics of HS patients evaluated with HRV (Heart Rate Variability).
AgeWeight (kg)Height (m)BMIPhysical conditionSmokeAge of onsetAge of diagnosisECG16721.6526.4OverweightNo1015Normal31951.732.9Obesity IYes2030Normal48901.827.8OverweightYes2045dx conduction disorder19721.7623.2NormoweightNo1418Normal32831.6530.5Obesity IYes2828Normal28761.5730.8Obesity IYes1525Normal241301.745Obesity IIIYes1323Normal21701.8520.5NormoweightNo1420NormalAs expected, from the supine to orthostatic position, RR interval showed a significant decrease (from796.3±110.8 to 674.68±56.95 msec—P<0.001); systolic and diastolic arterial pressure showed a slight increase as in normal subjects (Tables 4 and 5).Table 4
HRV values obtained using both linear and nonlinear methods in clinostatic position.
ClinoRR (ms)STD RR (ms)LF (%)HF (%)LF/HF (-)SD₁ (ms)SD₂ (ms)SD₂/SD₁α₁α₂α₁/α₂67524.57123.23.0612.632.30.391.4110.3024.67757.116.967.323.82.8311.520.90.551.2720.462.77965.534.333.362.10.5433.335.40.940.6270.3721.69795.836.640.951.50.7927.943.70.640.8390.4361.92766.226.138.849.70.7819.131.70.61.0150.4432.29860.642.329.662.30.4819.628.50.690.690.391.77642.810.455.6311.79712.90.541.0240.5251.95907.673.635.961.40.5857.187.10.660.8290.2034.08Mean796.333.146.645.61.423.536.60.610.42.6SD110.819.415.917.11.116.122.40.20.30.11.1Table 5
HRV values obtained using both linear and nonlinear methods in orthostatic position.
OrthoRR (ms)STD RR (ms)LF (%)HF (%)LF/HF (-)SD₁ (ms)SD₂ (ms)SD₁/SD₂α₁α₂α₁/α₂599.527.979.17.510.55938.40.231.7140.4863.53724.518.563.223.42.709.624.40.391.4010.5582.51746.730.357.634.41.6725.834.20.750.9750.3732.61741.825.347.746.91.0214.332.80.441.1260.4222.67669.820.263.920.73.091126.30.421.2990.632.06646.823.2865.316.236.832.10.211.7590.5223.37614.313.754.6361.527.617.80.431.0800.5761.8865441.854.439.31.3823.354.40.431.1980.3293.64Mean674.6825.1163.3126.694.7713.4332.550.411.320.492.78SD56.958.5813.0715.055.577.2610.940.160.290.100.66We also compared the HRV indexes of subjects with HS with those observed in a group of 14 normal subjects from our previous work. We observed a reduction of total HRV in the time and frequency domain: the STD showed a significant decrease in HS patients (40.32±5.79 to 33.1±19.4 msec). The sympathetic indexes in HS patients are higher than those in normal subjects even if not all of them reached statistical significance (LF/HF ratio 1.13±0.24 vs. 1.42±1.1, DFA α1 0.85±006 vs. 1.04±0.29).


## 4. Discussion

HS is an inflammatory skin disease of the terminal hair follicle. It manifests clinically with the onset of painful nodules, abscesses, and draining sinus tracts eventually resulting in scars, involving multiple regions of the body.Some recent studies have suggested a connection between HS and a significantly increased risk of ischemic cardiovascular events, such as myocardial infarction and stroke. It was also suggested that the risk of ischemic events is higher in subjects affected by HS when compared to severe psoriasis patients [14, 15]. The association between inflammatory markers and sympathetic hyperactivation indexes remained significant after correcting important risk factors and confounding factors such as age, Body Mass Index, heart rate, and smoking. Sympathetic activation can also be triggered by reflex mechanisms (arterial baroreceptor impairment), psychological stress, oxidative stress, obstructive sleep apnea, inflammation, and metabolic factors as I.R. and dysregulated production and secretion of adipokines from visceral fat with a particular important role of leptin [16].Interestingly, a Danish cross-sectional study showed that the mean heart rate in severe HS patients in resting conditions was suggestively higher when compared with that of controls [17]. Both HS and atherosclerosis have a chronic inflammatory subset, and cardiovascular and metabolic diseases are the most common comorbidities found in HS [18, 19]. Atherosclerosis, a chronic inflammatory disease affecting the blood vessels, is considered a major cause of cardiovascular diseases. Markers of inflammation are associated with both cardiovascular risk and the severity of HS inflammatory manifestations. It is known that inflammation represents the core in the bridge between HS and the pathogenesis of cardiovascular diseases, though many doubts remain regarding its mechanism of action. An important role is played by interleukin- (IL-) 17, which is found to be upregulated in lesional and perilesional skin of HS patients [20]. Overexpression of various other cytokines, such as IL-1, TNF-α, IL-10, and IL-11, as well as antimicrobial peptides (AMPs), beta-defensin 2, psoriasin, and cathelicidin, has been observed in the lesional setting.Interleukin- (IL-) 32 expression is significantly enhanced in plaque of atherosclerosis, and its proinflammatory function could help increase the risk of CV events in patients affected by chronic inflammatory diseases [21, 22]. Moreover, HS patients generally present a higher incidence of other cardiovascular risk factors such as obesity, smoking history, diabetes, hypertriglyceridemia, and metabolic syndrome [23].The present study was based on the analysis of HRV as a noninvasive tool to evaluate the association between HS and an increased CV risk. The mechanism leading to depressed HRV in heart failure is complex and not perfectly described, but it could partially be related to the alteration of neurohumoral activity [24].The results indicate that the sympathetic activity increases during the Tilt Test more than in normal subjects as shown by the tachycardia and by the increase of the arterial pressure and LF/HF ratio and the decrease of the STD. The Tilt Test induces, also, an increase in the properties of short-term fractal correlations of heart rate dynamics (DFA) accompanied with a decrease in all the nonlinear indexes of HRV, confirming the thesis that these indexes are an expression of the cardiac sympathetic-vagal balance.HS patients, if compared to the rest of the population, tend to present a higher incidence of carotid atherosclerosis. The role of smoking in the severity of HS is still not clear because it is not known whether smoking cessation improves the course of the disease.A higher resting heart rate and a lower heart rate variability are linked to subclinical chronic inflammation in adult and elderly patients. The higher death rate that was described in this kind of situations could therefore manifest a shared etiology. A disproportion of the sympathetic branch of the autonomic nervous system could correlate with inflammatory events to play a major role in atherosclerosis.It has been previously demonstrated that the use of HRV allows us to evidence an upregulated CV risk in other inflammatory dermatological disorders, such as psoriasis [25]. As reported in patients with psoriasis, in patients with HS, the sympathetic activation indexes (LF/HF, SD2/SD1, and α1/α2) are higher than those in normal subjects too.Studies conducted on patients with known HS indicate that many of them are overweight. Obesity is one of the main risk factors associated with the development of HS, although its role is yet to be defined. We did not observe any significant difference when reorganizing patients according to the BMI, although an increase of heart rate (RR) was observed in subjects with aBodyMassIndexBMI>25. Furthermore, we observed that the sympathetic activation indexes during passive orthostatism had increased compared to clinostatism.


## 5. Conclusions

Our data show a balanced reduction of the parasympathetic influence on the sinus node in patients with moderate to severe HS. These results could correlate with an increased risk of cardiovascular disease; thus, HS should be considered as an independent CV risk factor.Our data seems to indicate an increase in the sympathetic hyperactivity indices in HS patients, but since our study was based on a limited number of cases, we consider it to be preliminary and more studies are necessary in order to obtain statistical significance. Identification of individuals at risk for subsequent morbidity and nonrisk groups requires future prospective studies to determine the sensitivity, specificity, and predictive values of HRV.Moreover, our data demonstrates that biological therapy and other anti-inflammatory therapies could alter the autonomous cardiac function in patients with moderate to severe HS, as reported in a previous work on psoriasis [26]. The important connection between HS and cardiovascular diseases underlines the necessity of a cardiovascular screening in HS patients, particularly if other risk factors are present.Although the exact role of smoking in the pathogenesis of HS remains to be determined, HS could give us another opportunity to encourage our patients to change their lifestyle.Dermatologists must be conscious of the comorbidities related to HS. The aim is to create a condition in which patients are managed by the proper specialists. This may lead to an improvement in the disease course and in the quality of life of patients affected by this disease.



---

*Source: 1321782-2020-07-01.xml*