Spike detection in human muscle sympathetic nerve activity using the kurtosis of stationary wavelet transform coefficients

The accurate assessment of autonomic sympathetic function is important in the diagnosis and study of various autonomic and cardiovascular disorders. Sympathetic function in humans can be assessed by recording the muscle sympathetic nerve activity, which is characterized by synchronous neuronal discharges separated by periods of neural silence dominated by colored Gaussian noise. The raw nerve activity is generally rectified, integrated, and quantified using the integrated burst rate or area. We propose an alternative quantification involving spike detection using a two-stage stationary wavelet transform (SWT) de-noising method. The SWT coefficients are first separated into noise-related and burst-related coefficients on the basis of their local kurtosis. The noise-related coefficients are then used to establish a threshold to identify spikes within the bursts. This method demonstrated better detection performance than an unsupervised amplitude discriminator and similar wavelet-based methods when confronted with simulated data of varying burst rate and signal to noise ratio. Additional validation on data acquired during a graded head-up tilt protocol revealed a strong correlation between the mean spike rate and the mean integrate burst rate (r = 0.85) and burst area rate (r = 0.91). In conclusion, the kurtosis-based wavelet de-noising technique is a potentially useful method of studying sympathetic nerve activity in humans.