Autonomic nervous activity analysis based on visibility graph complex networks and skin sympathetic nerve activity

Background: Autonomic nerve system (ANS) plays an important role in regulating cardiovascular function and cerebrovascular function. Traditional heart rate variation (HRV) and emerging skin sympathetic nerve activity (SKNA) analyses from ultra-short-time (UST) data cannot fully reveal neural activity, thereby quantitatively reflect ANS intensity. Methods: Electrocardiogram and SKNA from sixteen patients (seven cerebral hemorrhage (CH) patients and nine control group (CO) patients) were recorded using a portable device. Ten derived HRV (mean, standard deviation and root mean square difference of sinus RR intervals (NNmean, SDNN and RMSSD), ultra-low frequency (<0.003 Hz, uLF), very low frequency ([0.003 Hz, 0.04 Hz), vLF), low frequency ([0.04 Hz, 0.15 Hz), LF) and high frequency power ([0.15 Hz, 0.4 Hz), HF), ratio of LF to HF (LF/HF), the standard deviation of instantaneous beat-to-beat R-R interval variability (SD1), and approximate entropy (ApEn)) and ten visibility graph (VG) features (diameter (Dia), average node degree (aND), average shortest-path length (aSPL), clustering coefficient (CC), average closeness centrality (aCC), transitivity (Trans), average degree centrality (aDC), link density (LD), sMetric (sM) and graph energy (GE) of the constructed complex network) were compared on 5-min and UST segments to verify their validity and robustness in discriminating CH and CO under different data lengths. Besides, their potential for quantifying ANS-Load were also investigated. Results: The validation results of HRV and VG features in discriminating CH from CO showed that VG features were more clearly distinguishable between the two groups than HRV features. For effectiveness evaluation of analyzing ANS on UST segment, the NNmean, SDNN, RMSSD, LF, HF and LF/HF in HRV features and the CC, Trans, Dia and GE of VG features remained stable in both activated and inactivated segments across all data lengths. The capability of HRV and VG features in quantifying ANS-Load were evaluated and compared under different ANS-Load, the results showed that most HRV features (SDNN, LFHF, RMSSD, vLF, LF and HF) and almost all VG features were correlated to sympathetic nerve activity intensity. Conclusions: The proposed autonomic nervous activity analysis method based on VG and SKNA offers a new insight into ANS assessment in UST segments and ANS-Load quantification.