- Research article
- Open Access
A method to quantify autonomic nervous system function in healthy, able-bodied individuals
Bioelectronic Medicine volume 7, Article number: 13 (2021)
The autonomic nervous system (ANS) maintains physiological homeostasis in various organ systems via parasympathetic and sympathetic branches. ANS function is altered in common diffuse and focal conditions and heralds the beginning of environmental and disease stresses. Reliable, sensitive, and quantitative biomarkers, first defined in healthy participants, could discriminate among clinically useful changes in ANS function. This framework combines controlled autonomic testing with feature extraction during physiological responses.
Twenty-one individuals were assessed in two morning and two afternoon sessions over two weeks. Each session included five standard clinical tests probing autonomic function: squat test, cold pressor test, diving reflex test, deep breathing, and Valsalva maneuver. Noninvasive sensors captured continuous electrocardiography, blood pressure, breathing, electrodermal activity, and pupil diameter. Heart rate, heart rate variability, mean arterial pressure, electrodermal activity, and pupil diameter responses to the perturbations were extracted, and averages across participants were computed. A template matching algorithm calculated scaling and stretching features that optimally fit the average to an individual response. These features were grouped based on test and modality to derive sympathetic and parasympathetic indices for this healthy population.
A significant positive correlation (p = 0.000377) was found between sympathetic amplitude response and body mass index. Additionally, longer duration and larger amplitude sympathetic and longer duration parasympathetic responses occurred in afternoon testing sessions; larger amplitude parasympathetic responses occurred in morning sessions.
These results demonstrate the robustness and sensitivity of an algorithmic approach to extract multimodal responses from standard tests. This novel method of quantifying ANS function can be used for early diagnosis, measurement of disease progression, or treatment evaluation.
This study registered with Clinicaltrials.gov, identifier NCT04100486. Registered September 24, 2019, https://www.clinicaltrials.gov/ct2/show/NCT04100486.
The autonomic nervous system (ANS) regulates and integrates the physiology of the heart, lung, spleen, intestines, and other organ systems. That physiology also includes control over blood vessels, pupils, perspiration, and salivary glands. Regulation depends on a balance between the sympathetic and parasympathetic systems, in which it is possible to monitor real-time ANS activity by recording neural activity from candidate cranial and peripheral nerves (Barman & Yates, 2017; Cracchiolo et al., 2019; Devor et al., 1994; Masi et al., 2019; Yoshimura et al., 1994; Zanos, 2019; Zanos et al., 2018). However, this daunting recording task would require implanted electrodes, a challenging prospect for animal experiments, no less clinical diagnosis and treatment. The direct measurements of ANS-dependent physiological signals offer obvious alternative and verifiable measures compared to invasive implant recording. These classes of measurements are now possible through advances in accepted, noninvasive clinical testing (Weimer, 2010). Specifically, standard techniques of autonomic testing include measuring heart rate (HR) and blood pressure (BP) during deep breathing (Coote & Chauhan, 2016; Russo et al., 2017; Shields, 2009), posture or tilt Table (Porta et al., 2007; Scheen & Philips, 2012; Yeh et al., n.d.), cold pressor (Allen et al., 1992; Doytchinova et al., 2017; Heath & Downey, 1990; Mourot et al., 2009; Wirch et al., n.d.), diving reflex (Hilz & Dütsch, 2006; Hilz et al., 1999), and the Valsalva maneuver (De Becker et al., 1998; Doytchinova et al., 2017; Gibbons et al., n.d.; Goldstein & Cheshire, 2017; Novak, 2011; Vogel et al., 2005). Sudomotor testing relies on thermoregulatory sweat testing (TST), quantitative sudomotor axon reflex testing (QSART), sympathetic skin response (SSR), silicone impressions, the acetylcholine sweat-spot test, and quantitative direct and indirect reflex testing (QDIRT), which also belong to the standard range of provocative autonomic tests (Illigens & Gibbons, 2008; Low et al., 2006; Pittenger et al., 2005; Sumner et al., 2003).
The potential data acquired from this range of tests, however, will fast exceed the size and homogeneity requirement for standard statistical evaluation. One early study developed a composite scoring scale that was able to detect generalized autonomic dysfunction, but it failed on discriminating by disorder or diagnosis (Low, 1993). While much attention has focused on measuring heart rate variability (HRV) as a proxy for vagal tone (Akselrod et al., 1981; Beckers et al., 2001; Ducla-Soares et al., 2007; Mainardi, 2009; Pagani et al., 1984), there is less evidence of the usefulness of this relationship (Billman, 2013; Ernst, 2017). HRV measures quantify fluctuations between inter-beat intervals (IBI); different time- and frequency-domain indices have been linked to short term (~five minutes) and 24-h metrics of sympathetic or parasympathetic activity (Shaffer & Ginsberg, 2017; Stavrakis et al., 2020). Although ultra-short term (less than five minutes) HRV measures are not common, it is reported that the root mean square of successive R-R interval differences (RMSSD), in as low as 30 s samples, may be a reliable surrogate for parasympathetic activity (Baek et al., 2015; Kang et al., 2016; Munoz et al., 2015; Salahuddin et al., 2007). Perhaps because HRV depends on underlying heart and respiratory rates that can be easily altered by diet, exercise, psychological stress, diurnal factors, and medication, as well as age and gender, its reliability as a surrogate marker has remained controversial (Kyriakou et al., 2019; Perini & Veicsteinas, 2003; Stein et al., 1997; Umetani et al., 1998; Yamamoto et al., 1991; Young & Benton, 2018; Zhang, 2007). Besides HRV, other ANS assessment techniques require drug induced response (baroreflex sensitivity), additional equipment (imaging or tilt tables), or invasive intraneural microelectrodes (direct muscle and sympathetic nerve activity measurement), making these assessments more difficult to administer (Goldberger et al., 2019; Stavrakis et al., 2020).
Thus, the dearth of effective, reliable, and reproducible data-driven approaches to quantify non-invasive recording modalities while administering a full battery of tests provokes an application of signal processing, machine learning techniques, and decoding algorithms on a set of commonly used clinical measurement of physiological signals in an attempt to derive a better understanding of autonomic function and dysregulation. Herein, we test whether these metrics of autonomic function and responses are sensitive and can reliably identify subtle ANS deviations in health, able bodied control participants. The findings demonstrated significant deviations of ANS responses correlated with body mass index (BMI), and also showed trends related to circadian rhythm. While BMI is not a gold standard index for autonomic activity, it is related to hemodynamic alterations that are linked with the ANS dysfunction (Guarino et al., 2017). If these noninvasive, quantitative autonomic function metrics can be used to retrieve or predict BMI from deviations in cardiovascular activity in a small sample size, this method might well enable objective measures of other disease states linked with ANS dysfunction and provide a useful tool for diagnosis and disease management.
This study recruited and enrolled 21 healthy, able-bodied participants between the ages of 18–60 years and a BMI less than 30. The mean age (±SD) of the participants was 29.9 (±6.5) years, with sixteen males and five females. The mean BMI (±SD) was 24.4 (±2.9). Exclusion criteria were: history of cardiac arrhythmia, coronary artery disease, autoimmune disease, chronic inflammatory disease, anemia, malignancy, depression, neurologic disease, diabetes mellitus, renal disease, dementia, psychiatric illness including active psychosis, or any other chronic medical condition, treatment with anti-cholinergic medication, current tobacco, nicotine, or other recreational drug use, pre-existing neurological disease, pregnancy, and implantable electronic devices. Participants were asked to fast and refrain from caffeine for at least four hours prior to testing. This study was approved by the Northwell Health Institutional Review Board, IRB #19–0461 and registered with Clinicaltrials.gov, identifier NCT04100486.
Autonomic testing sessions
Testing sessions occurred in a laboratory setting with shielded from external noise, light, or distractions, with an average humidity of 22% and average temperature of 22.1 °C. Lighting was set such that significant pupil changes were detected, averaging approximately 22 lm (Lux Light Meter Pro, developed by Marina Polyanskaya). Each participant attended four one-hour testing sessions: two in the morning (starting between 7 and 11 am) and two in the afternoon (starting between, 1 and 6 pm) over the course of two weeks. Five autonomic tests were performed during each session, with the order randomized for each participant and for each testing session (Fig. 1a); at least 1 min elapsed between tests, and if a test was incorrectly performed or there was a potential loss of data or signal, that test was repeated. In the squat test, the participant actively stood still and calm for one minute, sat in a single deep squat for one minute, and then stood again for the final minute, all in succession. For the cold pressor test, the participant immersed his or her right hand in ice water (< 5 °C) for at least 30 s and up to three minutes; after the first 30 s, the participant was allowed to remove their hand at their own will. During deep breathing, the participant followed a visual cue on a dimly lit computer monitor to time their respiratory rate to six breaths per minute for seven minutes. The diving reflex test was administered by placing a refrigerated gel-filled compress on the participant’s face for one minute, followed by one minute of recovery. Lastly, the Valsalva maneuver was performed by a forceful attempted exhalation, by bearing down and “expelling” air while keeping the mouth and nose shut for 15 s, followed by one minute of recovery. Each session also included a five-minute period of continuous transcutaneous auricular vagus nerve stimulation (taVNS), a noninvasive method of applying stimulation at the cymba conchae of the ear to target the auricular branch of the vagus nerve, applied during rest randomly within the series of the five autonomic tests; while physiological signals were collected for this stimulation period, the results are not reported as a part of this work (Debnath et al., n.d.). No patients experienced any vagal syncope events during testing sessions.
In each session, the participant’s cardiovascular data were continuously captured by noninvasive sensors transmitted through a data acquisition system (PowerLab 16/35, ADInstruments, Sydney, Australia) (Fig. 1b): six lead electrocardiography (ECG) (Lead Wires and Octal BioAmp, ADInstruments, Sydney, Australia), respiration (Respiratory Belt Transducer, ADInstruments, Sydney, Australia), and blood pressure (Human NIBP Nano System, Finapres Medical Systems, Enschede, The Netherlands). The six lead ECG wires were attached to four foam adhesive electrodes placed at each shoulder and each ankle, while the respiratory belt was wrapped and tightened around the torso. BP was recorded by a small non-invasive and inflatable cuff around the middle phalanx of the middle finger on the left hand. Additional sensors were attached to capture electrodermal activity (EDA) (dry, metal electrodes on two fingers) and pupil diameter (Tobii Pro Glasses 2, Tobii Pro AB, Stockholm, Sweden). All of the recordings were transmitted through that data acquisition system and software (LabChart, ADInstruments, Sydney, Australia) at a sampling rate of 1 kHz. All signals were marked and synchronized with experimental cues aligned with each test. Each participant’s BMI was recorded and monitored as the main characteristic to be retrieved from autonomic signals. Continuous autonomic data from each participant and multiple non-invasive sensors, including ECG, BP, respiration, EDA, and pupil diameter, were measured at multiple time points relative to rest and onset of the autonomic tests. Individual responses were characterized by an approach based on extracting unique features from recorded data to identify significant responses and compute physiologically relevant and plausible indices that are representative of deviations from resting ANS function; these methods were not linked with previously developed indices, tools, or features.
All signal processing, analysis, and statistical evaluation were performed using MATLAB R2019a (MathWorks, Natick, MA). A 60 Hz notch filter was applied to remove line noise from raw cardiovascular data signals. Pupil diameter data was calibrated based on gaze location (Supplementary Fig. 1); for all tests, participants were reminded to keep their eyes open and continue looking forward. Since illumination and perceived brightness can influence pupil size, measured pupil diameter was normalized to extract the autonomic response. Participants were not focused on any near objects, so pupil reactions due to accommodation were expected to be minimal and removed through averaging. The RGB values of the gaze location were converted to relative luminance values, and this was linearly fit to pupil diameter. By dividing by the slope of the line, the pupil sizes were normalized to account for effects of luminance.
Only modalities with a significant response to an autonomic test were considered for analysis. The mean μ and standard deviation σ were calculated for the baseline, taken in the same sitting posture before each autonomic test. A threshold for significance was determined as 1σ over baseline, a common way to define a threshold in signal detection theory (Merfeld, 2011); only peaks of the average response for a modality and test above this threshold were considered for analysis.
Each participant’s HR, mean arterial pressure (MAP), and RMSSD responses to the squat test, cold pressor test, diving reflex test, and Valsalva maneuver were characterized by a template matching method (Fig. 2). The comparisons between individual responses and the average response template in each epoch were quantified to reflect the autonomic response. An individual response may be delayed or occur sooner than the average response. Individual responses can also be shorter or longer in duration and smaller or larger in amplitude than an average response. To quantify these deviations, the parameters that minimize the following objective function were estimated:
In equation (1), the normalized sum of the squared error is minimized. f(x) is an average response template centered at time x0 and y is the trimmed individual response. The parameter a scales the duration of the response, and the parameter b delays or advances the response. The parameter c scales the amplitude of the template and will be reported as Vs, and d vertically offsets the template and will be reported as Vo. Since the average response of the template is normalized by the baseline period before each autonomic test, Vs and Vo can be combined and reported as a single parameter V = Vs + Vo that describes the net scaling of the individual response with respect to the template for each physiological signal for each test. The parameter H is reported to represent the parameter a to define the duration scaling for each individual response. Simulations were run based on the feature extraction process to ensure robustness to both initial conditions and additive noise and confirm curve fit performance.
For the deep breathing test, a robust parasympathetic test (Shields, 2009), only HR changes were considered; this is due to limited parasympathetic innervation to blood vessels that affect total resistance, and any BP response is primarily due to changes in HR and stroke volume. The differences between maximum and minimum HR for each breathing cycle in the first two minutes of the test were averaged for each individual. This value was divided by the average peak-to-peak HR over all individuals to calculate a scale-like feature.
Once features were calculated, they were separated by the type of autonomic response: sympathetic or parasympathetic. Based on literature (Kinoshita et al., 2006; Marfella et al., 1994; Sandroni et al., 1991; Shields, 2009; Victor et al., 1987) and observed data, responses were categorized by modality for each test, shown in Table 1. While both branches of the autonomic nervous system are active during each autonomic test, as reflected by changes in HR and BP, features were sorted as primarily driven by one branch. For example, responses that led to an increase in HR were considered as driven by a sympathetic nervous system increase (as opposed to a parasympathetic withdrawal), while slowed HRs were attributed as primarily induced by the parasympathetic nervous system. In the squat test, during the phase from squat to stand, only scale (V) was used to represent sympathetic activity while only stretch (H) was used to represent parasympathetic activity (Droguett & Santos, 2015; Du et al., 2005). During deep breathing, the calculated scale is used to represent parasympathetic activity. For all other tests and modalities, both scale and stretch features were used in tandem as surrogates for sympathetic or parasympathetic responses.
Modeling and validation
Classified signals were correlated with measured characteristics of the patient cohort, particularly participant BMI. Linear regression models were applied to fit the average of collected features from all modalities and testing epochs for each patient. Significant trends were validated by 10 repeats of 7-fold cross-validation; data from three patients were left out of the modeling for each fold, and samples were reshuffled for every repetition. The p-value was determined by a t-statistic; values less than 0.05 were considered significant.
Monitoring physiological signals to compute average responses
Cardiovascular, pupil dilation, and EDA signals were synchronously recorded for each participant, with raw recordings for ECG and BP used to calculate HR, HRV (RMSSD), and MAP (Fig. 3). The individual responses of calculated signals during each test were averaged to determine a response template (Fig. 4, Supplementary Fig. 2). For each test and each modality, baseline-normalized responses were averaged across all participants and sessions; taVNS did not have any significant effect on physiological responses (Debnath et al., n.d.). For HRV, we calculated the RMSSD feature for the three tests with intervals longer than 60 s, since this HRV measure requires 60 s of activity to be calculated accurately (Baek et al., 2015; Salahuddin et al., 2007); HR, MAP, pupil dilation and EDA measurements were averaged for all tests. The average responses for each modality and each test at epochs with phasic changes in activity were used as temples to compare to the individual responses for each modality and test. In the squat test, the epoch corresponded with posture changes from stand to squat and squat to stand. For the cold pressor test, the epoch reflected the participant immersing their hand in the ice water, while the epoch in the diving reflex test corresponded with the cold mask being placed onto the face. The Valsalva maneuver had four epochs that reflect dynamic changes in both HR and BP. For deep breathing, the epoch was the entire first two minutes of the task.
Extracted features from individual responses
From all the physiological modalities we monitored, participants’ cardiovascular measures (HR, MAP, and RMSSD) registered average responses with a peak above 1σ of their baseline, thus deemed significant responses (as detailed in Methods), while the pupil dilation and EDA measures did not show a consistent, significant response and were discarded from further processing and feature extraction (Supplementary Fig. 2). Figure 5 shows examples of individual responses, the calculated features, and the corresponding transformation of the average signal based on those features. The blue traces in each panel represent an individual’s phasic response within a single epoch for an autonomic test, while the black trace in each panel is the template from the average response over all individuals and all sessions for that modality and that test. After applying the template matching algorithm, the calculated duration scaling (H), vertical scaling (Vs), and the vertical offset (Vo) can transform the average template to the red trace, which shows the estimated fit in each case. The features illustrated in these examples are representative across individuals, tests, and modalities.
Closer inspection of individual responses can provide insight into physiological meaning of the vertical and horizontal scaling, as well as the vertical offset. As evident in the examples plotted in Fig. 5a and b, the features extracted during the squat test reveal certain response properties, specifically in the transition from squatting to standing (at t = 0), for HR and MAP, respectively. In the first example (Fig. 5a), HR response was faster (H < 1) and larger (Vs > 1) than the average response during this squatting-standing transition. Meanwhile, in the second example (Fig. 5b), the MAP response was also faster (H < 1) and larger (Vs > 1) than the average response.
Similar interpretations can be made from the responses during the diving reflex test (Fig. 5c-d), where t = 0 indicates the time that the refrigerated mask was placed on the participant’s face; 5c shows an example of a HR response, and 5d shows an example of a response in MAP. In both cases, the individual’s responses were faster (H < 1) and larger (Vs > 1) than the average response.
For the deep breathing test, the vertical scaling feature of the HR response scales the average peak-to-peak HR over all individuals based on the maximum and minimum HRs in each breathing cycle within the first two minutes of the test. The response in the specific example of Fig. 5e was more than twice as large as the average response (Vs = 2.27). As mentioned in the methods, since deep breathing was visually guided, all peaks were entrained to the respiratory cycle, and there was no need for horizontal scaling feature extraction.
Similarly, the features and fit of the RMSSD measure of HRV during the first two-minute interval of deep breathing reveal the vertical and horizontal scaling using the template matching algorithm; the RMSSD response features, as calculated for the specific example in Fig. 5f, reveal a faster (H < 1) and larger (Vs > 1) than the average RMSSD response.
Sympathetic and parasympathetic features and their correlations to BMI and session time
The strength of this analysis is that all of the features from all of the tests could be interpreted as vertical and horizontal scaling measurements. The data features were grouped by vertical and horizontal scaling and by sympathetic and parasympathetic drive. Next, the features were averaged for each participant over all four testing sessions (see Table 1). These averages (and corresponding standard deviations) were correlated with physiological characteristics (BMI and age). Figure 6 represents sympathetic features (top row), and parasympathetic features (bottom row), left and right columns represent H (duration scale) and V (amplitude scale), respectively. Each data point is the average (±SD) calculated feature for a single participant. A regression demonstrated a significant increase in sympathetic amplitude scaling with increasing BMI (p = 0.000377). Furthermore, this result was validated by 10 repeats of 7-fold cross-validation; the data left out was used to reproduce BMI based on the trend line produced. The mean absolute percentage error was approximately 11.7% with a standard deviation of 0.579.
To test for an effect of circadian rhythm, the grouped features were compared across morning and afternoon sessions. The average feature for each individual was calculated over two AM sessions and two PM sessions; in Fig. 7a, each trace represents a single individual, where blue and orange traces represents higher values in AM or PM sessions, respectively. While the differences were not significant, these trends are summarized in Fig. 7b; most (14 and 12 participants, respectively) participants had longer (H > 1) sympathetic and parasympathetic responses in PM sessions. For vertical scaling features, larger (V > 1) sympathetic responses were observed in PM sessions for 14 participants; parasympathetic vertical scaling features were larger in AM sessions for a higher number of participants (11 participants). Two participants did not have both morning and afternoon sessions and were left out of this analysis.
Recent studies have been focused on developing quantitative standards based on biomarkers to aid with diagnosis, prognosis, and estimates of treatment efficacy (Lötsch & Ultsch, 2018). To this end, this work sought to develop a reproducible and sensitive metric to quantify significant changes in ANS function by determining features that represented duration and amplitude scales compared to the average response. Instead of a priori developed indices, a template matching method was used to estimate indices that characterized physiologically-relevant indices that represent slight deviations in sympathetic and parasympathetic individual responses. These data show that sympathetic amplitude responses significantly increased and parasympathetic responses decreased with increasing BMI, and as such are consistent with past work that actually demonstrated significant sympathetic overactivity and decreased parasympathetic activity in patients with increased BMI greater than 25 (da Silva et al., 2009; Guarino et al., 2017; Molfino et al., 2009). While higher BMI may yield increased energy expenditure necessary for body weight, it has occurred to others that the activation of the sympathetic ANS may be an initial or primary driving force in weight maintenance regulation (Molfino et al., 2009). From a therapeutic standpoint, the links between high BMI and baroreceptor dysfunction, hypertension, organ damage, and cardiovascular disease suggest consideration of sympathetic inhibition.
Our data also show trends in ANS responses in relation with circadian rhythm. Sympathetic responses were both longer in duration and larger in amplitude during afternoon testing. For the parasympathetic features, longer duration responses occurred in afternoon testing, while larger amplitude responses were observed in morning testing sessions. Most past studies relating diurnal variations of autonomic function compare sympathetic and parasympathetic markers between day and night, and have demonstrated, especially in hypertensive patients, sympathetic activity spikes present in early morning (6:00 AM to 8:00 AM), with little change for the rest of the day and nighttime (Kario et al., 2004; Marfella et al., 2003; Middlekauff & Sontz, 1995; Panza et al., 1991). Although we did not test night time conditions and the recorded differences are not significant, our data suggest systematic variability cycles within the typical 24-h period, and as such, may also offer therapeutic opportunities with novel temporal characteristics.
Past studies of the ANS report HR and BP changes in response to stressor tests (Bellenger et al., 2016; Borresen & Lambert, 2008; Castiglioni et al., 2011; Freeman, 2006; Gonzaga et al., 2017; Martinez et al., 2010; McGuire et al., 2005; Michael et al., 2017; Pagani & Lucini, 2001; Radtke et al., 2016; Rossi et al., 2015; Taylor et al., 2003; Thayer et al., 2010; Van de Borne et al., 1994; Ziegler et al., 1992), but there is a dearth of studies that measure many of the remaining signals simultaneously and during controlled autonomic perturbations. The approach to ANS measurement taken in this work shows that parallel expansion of the modalities of raw physiological signals measured broadens the analysis from simple vitals to numerous measures that include temporal measures of HRV, pupillometry, respiration, and EDA.
This initial application of signal processing and machine learning on a set of standard physiological measures presents some challenges. First, this work depended only on cardiovascular changes during all autonomic tests. Other recordings, especially pupil diameter and EDA, were measured and analyzed in preliminary studies not discussed, but the peak average response did not meet the rigorous threshold of 1σ over baseline. Additionally, other autonomic measurements like photoplethysmogram (PPG), seismocardiogram (SCG), respiratory effort, pre-ejection period (PEP) of the heart, and other time-domain and frequency-domain HRV measures were not recorded. The reliability of our HRV measure, therefore, suggest that recording and analyzing these other responses might increase the reliability in defining a normal ANS index.
The focus on a small healthy cohort between the ages of 18 and 45 and BMI less than 30 provided a baseline that reflects the significant differences and supports the plausibility that future recording efforts with larger and more diverse sample populations will be fruitful. Specifically, significant responses were not always measured for this healthy population in sympathetic or parasympathetic duration or amplitude scales, but there may be observable changes in patient populations or individuals outside of these age and BMI constraints. Relative to the average template from a healthy cohort, any patterns of significantly elevated, diminished, accelerated, or slowed responses to the same autonomic tests in a patient cohort could be used to diagnose specific dysautonomias. The limitations of our sample extend to the limitation of the recording modalities collected but not analyzed. For example, an impaired EDA has been linked to early stages of diabetic neuropathy (Khalfallah et al., 2012; Petrofsky & McLellan, 2009), and impaired HR, EDA, and pupil dilation have been linked with post-traumatic stress disorder (PTSD) (Mckinnon et al., 2020; Pole, 2007), responses that would not have been observed here. Nevertheless, the methods that focus on a healthy cohort and cardiovascular signals generated a reasonable set of healthy responses that may be used to construct a normal template. Such a range of normal responses might well be used to quantify autonomic variations in actual patient populations.
A third limitation is that this method relies on phasic response to a specific battery of autonomic tests, as opposed to continuous changes in recording modalities. Because the method is not ongoing, it cannot provide a continuously updating biomarker to estimate ANS function and balance. Additionally, this battery focuses on cardiovascular responses due to peripheral signaling in the body, and may explain the lack of robust and significant responses in EDA and pupil diameter. Future work may also include other tests, such as a cognitive aptitude assessment, social stress test or mental arithmetic (Bauerly et al., 2019; Duschek et al., 2009; Gurel et al., 2020a; Tornatzky & Miczek, 1993), that can induce more central nervous system mediated responses and, therefore, increase the utility of EDA, pupil diameter, and other recorded modalities. Recent work has also shown that biomarkers, like HR and PPG amplitude, can be used to predict responses to transcutaneous cervical vagus nerve stimulation (tcVNS) and model dynamic characteristics of an adapting ANS (Gazi et al., 2020). While it is unclear how this may scale for other conditions or interventions, modeling biomarker responses can be applied to continuously monitoring vital-sensing devices to calculate real-time risk scores and further comprehensive index values related to autoimmune health.
This method was developed toward creating data-driven approaches to comprehensively and objectively quantify the ANS. Modern methods of computational science, including machine learning and artificial intelligence techniques, have been used to decode complex clinical and experimental data by detecting patterns, classifying signals, and extracting information to inform diagnostic and treatment actions (Debnath et al., 2020; Norgeot et al., 2019; Wiens & Shenoy, 2018). Continuous data from many sensors, including those in this study and adding electroencephalography (EEG) or other neural recording devices, can be used to train such a model on recordings from healthy, able-bodied individuals to characterize ANS balance. With so many features extracted from each patient, advanced classification strategies can be built to explore non-linear relationships between features and calculate weighted indices to quantify the sympathetic and parasympathetic branches in patient populations. Since the battery of autonomic tests can be completed within 30 min and all sensors are non-invasively placed, clinical translation can be simplified. Future studies could focus on varying patient populations, as disturbances in autonomic regulation have been described in a variety of diseases, including those resulting from focal injury, such as spinal cord injuries (Krassioukov et al., 2012) and stroke (Dütsch et al., 2007), and diffuse disorders, such as sepsis and infection (Badke et al., 2018; Ferreira & Bissell, 2018), rheumatoid arthritis (Koopman et al., 2016a; Koopman et al., 2016b), Crohn’s disease (Engel et al., 2015), and diabetes mellitus (Serhiyenko & Serhiyenko, 2018; Verrotti et al., 2014). Additionally, dysautonomias have been described in numerous cardiovascular conditions (Broadstone et al., 1991; Carthy, 2013; Kishi, 2012; Shen & Zipes, 2014; Vinik et al., 2013) and central nervous system disorders, including Alzheimer’s disease (Femminella et al., 2014), Parkinson’s disease (Goldstein, 2014), Huntington’s disease (Diago et al., 2017; Kobal et al., 2010), and psychiatric conditions including depression, schizophrenia, and PTSD (Jung et al., 2019), among others. Another possible application could be evaluation of targeted neuromodulation; specifically, there is clinical interest in stimulating the vagus nerve, which is involved with responses in cardiovascular, pulmonary, gastrointestinal, renal, hepatic, and endocrine systems (Chavan et al., 2017; Pavlov et al., 2018). Vagus nerve stimulation (VNS) has been used in previous studies for multiple conditions, including refractory epilepsy (Rong et al., 2014; Stefan et al., 2012), depression (Kong et al., 2018; Rong et al., 2016), PTSD (Bremner et al., 2020), pre-diabetes (Huang et al., 2014), tinnitus (Shim et al., 2015), stroke (Redgrave et al., 2018), and others, including oromotor dysfunction, rheumatoid arthritis, and obesity (Guiraud et al., 2016). These studies have used a range of electrical stimulation settings and sites, and there is no optimal dose or set of parameters (Badran et al., 2018); the mechanism of VNS and responses are not well understood. While there have been studies that report mixed or no reported significant effects of VNS on HRV, pupil diameter, and evoked potentials, the preliminary effects on clinical populations are clear (Burger et al., 2020; Gurel et al., 2020b; Libbus et al., 2017). The taVNS protocol used in this study (pulse width of 300 μs at a continuous rate of 30 Hz for 5 min) did not produce any significant effects (Debnath et al., n.d.), but short-term taVNS (3.4 s ON, 26–27 s OFF) elicited robust pupil dilation and alpha oscillations compared to sham stimulation, showing that significant pupillary and EEG markers can be observed (Sharon et al., 2020). A set of biomarkers or calculated features to accurately and consistently quantify these changes and relate them to the ANS in a patient-centered approach can be extremely helpful in a number of clinical applications.
The sympathetic and parasympathetic parameters determined in this study will be valuable to diagnose autonomic function and underlying disorders, as well as predict responses to targeted modulated therapies. While the use of autonomic modulation has shown promise in treating cardiovascular, autoimmune, and nervous system disorders, the template matching method in this work can offer additional insight towards the effects of stimulation and medication in patient populations.
Availability of data and materials
The data that support the findings of this study are available from the corresponding author upon reasonable request.
autonomic nervous system
heart rate variability
mean arterial pressure
post-traumatic stress disorder
quantitative direct and indirect reflex test
quantitative sudomotor axon reflex test
root mean square of successive R-R interval differences
sympathetic skin response
transcutaneous auricular vagus nerve stimulation
transcutaneous cervical vagus nerve stimulation
thermoregulatory sweat testing
vagus nerve stimulation
Akselrod S, Gordon D, Ubel FA, Shannon DC, Barger AC, Cohen RJ. Power spectrum analysis of heart rate fluctuation: a quantitative probe of beat-to-beat cardiovascular control. Science. 1981;213(4504):220–2. https://doi.org/10.1126/science.6166045.
Allen MT, Shelley KS, Boquet AJ. A comparison of cardiovascular and autonomic adjustments to three types of cold stimulation tasks. Int J Psychophysiol. 1992;13(1):59–69. https://doi.org/10.1016/0167-8760(92)90021-3.
Badke CM, Marsillio LE, Weese-Mayer DE, Sanchez-Pinto LN. Autonomic nervous system dysfunction in pediatric Sepsis. Front Pediatr. 2018;6:280. https://doi.org/10.3389/fped.2018.00280.
Badran BW, Mithoefer OJ, Summer CE, LaBate NT, Glusman CE, Badran AW, et al. Short trains of transcutaneous auricular vagus nerve stimulation (taVNS) have parameter-specific effects on heart rate. Brain Stimul. 2018;11(4):699–708. https://doi.org/10.1016/j.brs.2018.04.004.
Baek HJ, Cho C-H, Cho J, Woo J-M. Reliability of ultra-short-term analysis as a surrogate of standard 5-min analysis of heart rate variability. Telemed J E Health. 2015;21(5):404–14. https://doi.org/10.1089/tmj.2014.0104.
Barman SM, Yates BJ. Deciphering the neural control of sympathetic nerve activity: status report and directions for future research. Front Neurosci. 2017;11:730. https://doi.org/10.3389/fnins.2017.00730.
Bauerly KR, Jones RM, Miller C. Effects of social stress on autonomic, behavioral, and acoustic parameters in adults who stutter. Journal of Speech, Language, and Hearing Research. 2019;62(7):2185–202. https://doi.org/10.1044/2019_JSLHR-S-18-0241.
Beckers F, Ramaekers D, Aubert AE. Approximate entropy of heart rate variability: validation of methods and application in heart failure. Cardiovasc Eng. 2001;1(4):177–82. https://doi.org/10.1023/A:1015212328405.
Bellenger CR, Fuller JT, Thomson RL, Davison K, Robertson EY, Buckley JD. Monitoring athletic training status through autonomic heart rate regulation: a systematic review and Meta-analysis. Sports Med. 2016;46(10):1461–86. https://doi.org/10.1007/s40279-016-0484-2.
Billman GE. The LF/HF ratio does not accurately measure cardiac sympatho-vagal balance. Front Physiol. 2013;4. https://doi.org/10.3389/fphys.2013.00026.
Borresen J, Lambert MI. Autonomic control of heart rate during and after exercise. Sports Med. 2008;38(8):633–46. https://doi.org/10.2165/00007256-200838080-00002.
Bremner JD, Gurel NZ, Wittbrodt MT, Shandhi MH, Rapaport MH, Nye JA, et al. Application of Noninvasive Vagal Nerve Stimulation to Stress-Related Psychiatric Disorders. J Pers Med. 2020;10(3). https://doi.org/10.3390/jpm10030119.
Broadstone VL, Roy T, Self M, Pfeifer MA. Cardiovascular autonomic dysfunction: diagnosis and prognosis. Diabet Med 8 Spec No: S88–93, 1991. https://doi.org/10.1111/j.1464-5491.1991.tb02165.x, S2, S93.
Burger AM, Van der Does W, Brosschot JF, Verkuil B. From ear to eye? No effect of transcutaneous vagus nerve stimulation on human pupil dilation: a report of three studies. Biol Psychol. 2020;152:107863. https://doi.org/10.1016/j.biopsycho.2020.107863.
Carthy ER. Autonomic dysfunction in essential hypertension: a systematic review. Ann Med Surg (Lond). 2013;3(1):2–7. https://doi.org/10.1016/j.amsu.2013.11.002.
Castiglioni P, Parati G, Rienzo MD, Carabalona R, Cividjian A, Quintin L. Scale exponents of blood pressure and heart rate during autonomic blockade as assessed by detrended fluctuation analysis. J Physiol. 2011;589:355–69. https://doi.org/10.1113/jphysiol.2010.196428.
Chavan SS, Pavlov VA, Tracey KJ. Mechanisms and therapeutic relevance of neuro-immune communication. Immunity. 2017;46(6):927–42. https://doi.org/10.1016/j.immuni.2017.06.008.
Coote JH, Chauhan RA. The sympathetic innervation of the heart: important new insights. Auton Neurosci. 2016;199:17–23. https://doi.org/10.1016/j.autneu.2016.08.014.
Cracchiolo M, Sacramento JF, Mazzoni A, Panarese A, Carpaneto J, Conde SV, et al. Decoding neural metabolic markers from the carotid sinus nerve in a type 2 diabetes model. IEEE Transactions on Neural Systems and Rehabilitation Engineering. 2019;27(10):2034–43. https://doi.org/10.1109/TNSRE.2019.2942398.
De Becker P, Dendale P, De Meirleir K, Campine I, Vandenborne K, Hagers Y. Autonomic testing in patients with chronic fatigue syndrome. Am J Med. 1998;105(3):22S–6S. https://doi.org/10.1016/S0002-9343(98)00168-5.
Debnath S, Barnaby DP, Coppa K, Makhnevich A, Kim EJ, Chatterjee S, et al. Machine learning to assist clinical decision-making during the COVID-19 pandemic. Bioelectron Med. 2020;6:14. https://doi.org/10.1186/s42234-020-00050-8.
Debnath S, Levy TJ, Zanos S, Zanos TP. Noninvasive, multimodal assessment of physiological responses to transcutaneous auricular vagus nerve stimulation. 2021 10th International IEEE/EMBS Conference on Neural Engineering. 2021;1109–12. https://doi.org/10.1109/NER49283.2021.9441419.
Devor M, Jänig W, Michaelis M. Modulation of activity in dorsal root ganglion neurons by sympathetic activation in nerve-injured rats. J Neurophysiol. 1994;71(1):38–47. https://doi.org/10.1152/jn.19188.8.131.52.
Diago E, Pérez J, Lasaosa S, Alebesque A, Horta S, Kulisevsky J, et al. Circadian rhythm and autonomic dysfunction in presymptomatic and early Huntington’s disease. Parkinsonism Relat Disord. 2017;44:95–100. https://doi.org/10.1016/j.parkreldis.2017.09.013.
Doytchinova A, Hassel JL, Yuan Y, Lin H, Yin D, Adams D, et al. Simultaneous noninvasive recording of skin sympathetic nerve activity and electrocardiogram. Heart Rhythm. 2017;14(1):25–33. https://doi.org/10.1016/j.hrthm.2016.09.019.
Droguett VSL, Santos A da C, Medeiros CE de, Marques DP, Nascimento LS do, Brasileiro-Santos M do S. Cardiac autonomic modulation in healthy elderly after different intensities of dynamic exercise. Clinical Interventions in Aging 10 Dove Press. 2015:203–8.
Du N, Bai S, Oguri K, Kato Y, Matsumoto I, Kawase H, et al. Heart rate recovery after exercise and neural regulation of heart rate variability in 30-40 year old female Marathon runners. J Sports Sci Med. 2005;4(1):9–17.
Ducla-Soares JL, Santos-Bento M, Laranjo S, Andrade A, Ducla-Soares E, Boto JP, et al. Wavelet analysis of autonomic outflow of normal subjects on head-up tilt, cold pressor test. Valsalva manoeuvre and deep breathing Experimental Physiology. 2007;92(4):677–86. https://doi.org/10.1113/expphysiol.2007.038026.
Duschek S, Muckenthaler M, Werner N, Reyes del Paso GA. Relationships between features of autonomic cardiovascular control and cognitive performance. Biol Psychol. 2009;81:110–7. https://doi.org/10.1016/j.biopsycho.2009.03.003, 2.
Dütsch M, Burger M, Dörfler C, Schwab S, Hilz MJ. Cardiovascular autonomic function in poststroke patients. Neurology. 2007;69(24):2249–55. https://doi.org/10.1212/01.wnl.0000286946.06639.a7.
Engel T, Ben-Horin S, Beer-Gabel M. Autonomic dysfunction correlates with clinical and inflammatory activity in patients with Crohn’s disease. Inflamm Bowel Dis. 2015;21:2320–6. https://doi.org/10.1097/MIB.0000000000000508.
Ernst G. Heart-rate variability-more than heart beats? Front Public Health. 2017;5:240. https://doi.org/10.3389/fpubh.2017.00240.
Femminella GD, Rengo G, Komici K, Iacotucci P, Petraglia L, Pagano G, et al. Autonomic dysfunction in Alzheimer’s disease: tools for assessment and review of the literature. J Alzheimers Dis. 2014;42(2):369–77. https://doi.org/10.3233/JAD-140513.
Ferreira JA, Bissell BD. Misdirected sympathy: the role of Sympatholysis in Sepsis and septic shock. J Intensive Care Med. 2018;33(2):74–86. https://doi.org/10.1177/0885066616689548.
Freeman R. Assessment of cardiovascular autonomic function. Clin Neurophysiol. 2006;117(4):716–30. https://doi.org/10.1016/j.clinph.2005.09.027.
Gazi AH, Gurel NZ, Richardson KLS, Wittbrodt MT, Shah AJ, Vaccarino V, et al. Digital cardiovascular biomarker responses to transcutaneous cervical Vagus nerve stimulation: state-space modeling, prediction, and simulation. JMIR mHealth and uHealth. 2020;8(9):e20488. https://doi.org/10.2196/20488.
Gibbons CH, Cheshire WP, Fife TD. Model Coverage Policy: Autonomic Testing [Online]. https://www.aan.com/siteassets/home-page/tools-and-resources/practicing-neurologist%2D%2Dadministrators/billing-and-coding/model-coverage-policies/14autonomicmodel_tr.pdf.
Goldberger JJ, Arora R, Buckley U, Shivkumar K. Autonomic nervous system dysfunction: JACC focus seminar. J Am Coll Cardiol. 2019;73(10):1189–206. https://doi.org/10.1016/j.jacc.2018.12.064.
Goldstein DS. Dysautonomia in Parkinson’s disease: neurocardiological abnormalities. Compr Physiol. 2014;4:805–26. https://doi.org/10.1002/cphy.c130026.
Goldstein DS, Cheshire WP. Beat-to-beat blood pressure and heart rate responses to the Valsalva maneuver. Clin Auton Res. 2017;27(6):361–7. https://doi.org/10.1007/s10286-017-0474-y.
Gonzaga LA, Vanderlei LCM, Gomes RL, Valenti VE. Caffeine affects autonomic control of heart rate and blood pressure recovery after aerobic exercise in young adults: a crossover study. Sci Rep. 2017;7(1):14091. https://doi.org/10.1038/s41598-017-14540-4.
Guarino D, Nannipieri M, Iervasi G, Taddei S, Bruno RM. The role of the autonomic nervous system in the pathophysiology of obesity. Front Physiol. 2017;8. https://doi.org/10.3389/fphys.2017.00665.
Guiraud D, Andreu D, Bonnet S, Carrault G, Couderc P, Hagège A, et al. Vagus nerve stimulation: state of the art of stimulation and recording strategies to address autonomic function neuromodulation. J Neural Eng. 2016;13(4):041002. https://doi.org/10.1088/1741-2560/13/4/041002.
Gurel NZ, Huang M, Wittbrodt MT, Jung H, Ladd SL, MdMH S, et al. Quantifying acute physiological biomarkers of transcutaneous cervical vagal nerve stimulation in the context of psychological stress. Brain Stimulation. 2020a;13(1):47–59. https://doi.org/10.1016/j.brs.2019.08.002.
Gurel NZ, Wittbrodt MT, Jung H, Shandhi MH, Driggers EG, Ladd SL, et al. Transcutaneous cervical vagal nerve stimulation reduces sympathetic responses to stress in posttraumatic stress disorder: A double-blind, randomized, sham controlled trial. Neurobiology of Stress. 2020b;13:100264. https://doi.org/10.1016/j.ynstr.2020.100264.
Heath ME, Downey JA. The cold face test (diving reflex) in clinical autonomic assessment: methodological considerations and repeatability of responses. Clin Sci (Lond). 1990;78(2):139–47. https://doi.org/10.1042/cs0780139.
Hilz MJ, Dütsch M. Quantitative studies of autonomic function. Muscle Nerve. 2006;33(1):6–20. https://doi.org/10.1002/mus.20365.
Hilz MJ, Stemper B, Sauer P, Haertl U, Singer W, Axelrod FB. Cold face test demonstrates parasympathetic cardiac dysfunction in familial dysautonomia. American journal of physiology-regulatory. Integrative and Comparative Physiology. 1999;276(6):R1833–9. https://doi.org/10.1152/ajpregu.1999.276.6.R1833.
Huang F, Dong J, Kong J, Wang H, Meng H, Spaeth RB, et al. Effect of transcutaneous auricular vagus nerve stimulation on impaired glucose tolerance: a pilot randomized study. BMC Complement Altern Med. 2014;14(1):203. https://doi.org/10.1186/1472-6882-14-203.
Illigens BMW, Gibbons CH. Sweat testing to evaluate autonomic function. Clin Auton Res. 2008;19(2):79–87. https://doi.org/10.1007/s10286-008-0506-8.
Jung W, Jang K-I, Lee S-H. Heart and brain interaction of psychiatric illness: a review focused on heart rate variability, cognitive function, and Quantitative electroencephalography. Clin Psychopharmacol Neurosci. 2019;17(4):459–74. https://doi.org/10.9758/cpn.2019.17.4.459.
Kang JH, Kim JK, Hong SH, Lee CH, Choi BY. Heart rate variability for quantification of autonomic dysfunction in fibromyalgia. Ann Rehabil Med. 2016;40(2):301–9. https://doi.org/10.5535/arm.2016.40.2.301.
Kario K, Pickering TG, Hoshide S, Eguchi K, Ishikawa J, Morinari M, et al. Morning blood pressure surge and hypertensive cerebrovascular disease: role of the alpha adrenergic sympathetic nervous system. Am J Hypertens. 2004;17(8):668–75. https://doi.org/10.1016/j.amjhyper.2004.04.001.
Khalfallah K, Ayoub H, Calvet JH, Neveu X, Brunswick P, Griveau S, et al. Noninvasive galvanic skin sensor for early diagnosis of Sudomotor dysfunction: application to diabetes. IEEE Sensors J. 2012;12(3):456–63. https://doi.org/10.1109/JSEN.2010.2103308.
Kinoshita T, Nagata S, Baba R, Kohmoto T, Iwagaki S. Cold-water face immersion per se elicits cardiac parasympathetic activity. Circ J. 2006;70(6):773–6. https://doi.org/10.1253/circj.70.773.
Kishi T. Heart failure as an autonomic nervous system dysfunction. J Cardiol. 2012;59(2):117–22. https://doi.org/10.1016/j.jjcc.2011.12.006.
Kobal J, Melik Z, Cankar K, Bajrovic FF, Meglic B, Peterlin B, et al. Autonomic dysfunction in presymptomatic and early symptomatic Huntington’s disease. Acta Neurol Scand. 2010;121(6):392–9. https://doi.org/10.1111/j.1600-0404.2009.01251.x.
Kong J, Fang J, Park J, Li S, Rong P. Treating depression with transcutaneous auricular Vagus nerve stimulation: state of the art and future perspectives. Front Psychiatry. 2018;9. https://doi.org/10.3389/fpsyt.2018.00020.
Koopman FA, Chavan SS, Miljko S, Grazio S, Sokolovic S, Schuurman PR, et al. Vagus nerve stimulation inhibits cytokine production and attenuates disease severity in rheumatoid arthritis. Proc Natl Acad Sci U S A. 2016a;113(29):8284–9. https://doi.org/10.1073/pnas.1605635113.
Koopman FA, Tang MW, Vermeij J, de Hair MJ, Choi IY, Vervoordeldonk MJ, et al. Autonomic dysfunction precedes development of rheumatoid arthritis: a prospective cohort study. EBioMedicine. 2016b;6:231–7. https://doi.org/10.1016/j.ebiom.2016.02.029.
Krassioukov A, Biering-Sørensen F, Donovan W, Kennelly M, Kirshblum S, Krogh K, et al. International standards to document remaining autonomic function after spinal cord injury. J Spinal Cord Med. 2012;35(4):202–11. https://doi.org/10.1179/1079026812Z.00000000053.
Kyriakou K, Resch B, Sagl G, Petutschnig A, Werner C, Niederseer D, et al. Detecting moments of stress from measurements of wearable physiological sensors. Sensors (Basel). 2019;19(17). https://doi.org/10.3390/s19173805.
Libbus I, Nearing BD, Amurthur B, KenKnight BH, Verrier RL. Quantitative evaluation of heartbeat interval time series using Poincaré analysis reveals distinct patterns of heart rate dynamics during cycles of vagus nerve stimulation in patients with heart failure. J Electrocardiol. 2017;50(6):898–903. https://doi.org/10.1016/j.jelectrocard.2017.06.007.
Lötsch J, Ultsch A. Machine learning in pain research. Pain. 2018;159(4):623–30. https://doi.org/10.1097/j.pain.0000000000001118.
Low PA. Composite autonomic scoring scale for laboratory quantification of generalized autonomic failure. Mayo Clin Proc. 1993;68(8):748–52. https://doi.org/10.1016/S0025-6196(12)60631-4.
Low VA, Sandroni P, Fealey RD, Low PA. Detection of small-fiber neuropathy by sudomotor testing. Muscle Nerve. 2006;34(1):57–61. https://doi.org/10.1002/mus.20551.
Mainardi LT. On the quantification of heart rate variability spectral parameters using time-frequency and time-varying methods. Philosophical transactions of the Royal Society a: mathematical. Phys Eng Sci. 2009;367(1887):255–75. https://doi.org/10.1098/rsta.2008.0188.
Marfella R, Giugliano D, di Maro G, Acampora R, Giunta R, D’Onofrio F. The squatting test: a useful tool to assess both parasympathetic and sympathetic involvement of the cardiovascular autonomic neuropathy in diabetes. Diabetes. 1994;43:607–12. https://doi.org/10.2337/diab.43.4.607.
Marfella R, Gualdiero P, Siniscalchi M, Carusone C. Verza, Marzano S, Esposito K, Giugliano D. morning blood pressure peak, QT intervals, and sympathetic activity in hypertensive patients. Hypertension. 2003;41(2):237–43. https://doi.org/10.1161/01.HYP.0000050651.96345.0E.
Martinez JM, Garakani A, Kaufmann H, Aaronson CJ, Gorman JM. Heart rate and blood pressure changes during autonomic nervous system challenge in panic disorder patients. Psychosom Med. 2010;72(5):442–9. https://doi.org/10.1097/PSY.0b013e3181d972c2.
Masi EB, Levy T, Tsaava T, Bouton CE, Tracey KJ, Chavan SS, et al. Identification of hypoglycemia-specific neural signals by decoding murine vagus nerve activity. Bioelectronic Medicine. 2019;5(1):9. https://doi.org/10.1186/s42234-019-0025-z.
McGuire KMB, Greenberg MA, Gevirtz R. Autonomic effects of expressive writing in individuals with elevated blood pressure. J Health Psychol. 2005;10(2):197–209. https://doi.org/10.1177/1359105305049767.
Mckinnon AI, Gray NS, Snowden RJ. Enhanced emotional response to both negative and positive images in post-traumatic stress disorder: evidence from pupillometry. Biol Psychol. 2020;154:107922. https://doi.org/10.1016/j.biopsycho.2020.107922.
Merfeld DM. Signal detection theory and vestibular thresholds: I. Basic theory and practical considerations Exp Brain Res. 2011;210:389–405. https://doi.org/10.1007/s00221-011-2557-7, 3-4.
Michael S, Graham KS, Davis GMO. Cardiac autonomic responses during exercise and post-exercise recovery using heart rate variability and systolic time intervals—a review. Front Physiol. 2017;8. https://doi.org/10.3389/fphys.2017.00301.
Middlekauff HR, Sontz EM. Morning sympathetic nerve activity is not increased in humans. Implications for mechanisms underlying the circadian pattern of cardiac risk. Circulation. 1995;91(10):2549–55. https://doi.org/10.1161/01.cir.91.10.2549.
Molfino A, Fiorentini A, Tubani L, Martuscelli M, Fanelli FR, Laviano A. Body mass index is related to autonomic nervous system activity as measured by heart rate variability. Eur J Clin Nutr. 2009;63(10):1263–5. https://doi.org/10.1038/ejcn.2009.35.
Mourot L, Bouhaddi M, Regnard J. Effects of the cold pressor test on cardiac autonomic control in normal subjects. Physiol Res. 2009;58(1):83–91. https://doi.org/10.33549/physiolres.931360.
Munoz ML, van Roon A, Riese H, Thio C, Oostenbroek E, Westrik I, et al. Validity of (ultra-)short recordings for heart rate variability measurements. PLoS One. 2015;10(9):e0138921. https://doi.org/10.1371/journal.pone.0138921.
Norgeot B, Glicksberg BS, Butte AJ. A call for deep-learning healthcare. Nat Med. 2019;25(1):14–5. https://doi.org/10.1038/s41591-018-0320-3.
Novak P. Quantitative Autonomic Testing. JoVE (Journal of Visualized Experiments). 2011:e2502. https://doi.org/10.3791/2502.
Pagani M, Lucini D. Autonomic dysregulation in essential hypertension: insight from heart rate and arterial pressure variability. Auton Neurosci. 2001;90(1-2):76–82. https://doi.org/10.1016/S1566-0702(01)00270-3.
Pagani M, Lombardi F, Guzzetti S, Sandrone G, Rimoldi O, Malfatto G, et al. Power spectral density of heart rate variability as an index of sympatho-vagal interaction in normal and hypertensive subjects. J Hypertens Suppl. 1984;2:S383–5.
Panza JA, Epstein SE, Quyyumi AA. Circadian variation in vascular tone and its relation to alpha-sympathetic vasoconstrictor activity. N Engl J Med. 1991;325(14):986–90. https://doi.org/10.1056/NEJM199110033251402.
Pavlov VA, Chavan SS, Tracey KJ. Molecular and functional neuroscience in immunity. Annu Rev Immunol. 2018;36(1):783–812. https://doi.org/10.1146/annurev-immunol-042617-053158.
Perini R, Veicsteinas A. Heart rate variability and autonomic activity at rest and during exercise in various physiological conditions. Eur J Appl Physiol. 2003;90(3-4):317–25. https://doi.org/10.1007/s00421-003-0953-9.
Petrofsky JS, McLellan K. Galvanic skin resistance—a marker for endothelial damage in diabetes. Diabetes Technol Ther. 2009;11(7):461–7. https://doi.org/10.1089/dia.2008.0096.
Pittenger GL, Mehrabyan A, Simmons K, Null A, Dublin C, Barlow P, et al. Small fiber neuropathy is associated with the metabolic syndrome. Metab Syndr Relat Disord. 2005;3(2):113–21. https://doi.org/10.1089/met.2005.3.113.
Pole N. The psychophysiology of posttraumatic stress disorder: a meta-analysis. Psychol Bull. 2007;133(5):725–46. https://doi.org/10.1037/0033-2909.133.5.725.
Porta A, Tobaldini E, Guzzetti S, Furlan R, Montano N, Gnecchi-Ruscone T. Assessment of cardiac autonomic modulation during graded head-up tilt by symbolic analysis of heart rate variability. Am J Phys Heart Circ Phys. 2007;293(1):H702–8. https://doi.org/10.1152/ajpheart.00006.2007.
Radtke T, Poerschke D, Wilhelm M, Trachsel LD, Tschanz H, Matter F, et al. Acute effects of Finnish sauna and cold-water immersion on haemodynamic variables and autonomic nervous system activity in patients with heart failure. Eur J Prev Cardiol. 2016;23(6):593–601. https://doi.org/10.1177/2047487315594506.
Redgrave JN, Moore L, Oyekunle T, Ebrahim M, Falidas K, Snowdon N, et al. Transcutaneous auricular Vagus nerve stimulation with concurrent upper limb repetitive task practice for Poststroke motor recovery: a pilot study. J Stroke Cerebrovasc Dis. 2018;27(7):1998–2005. https://doi.org/10.1016/j.jstrokecerebrovasdis.2018.02.056.
Rong P, Liu A, Zhang J, Wang Y, He W, Yang A, et al. Transcutaneous vagus nerve stimulation for refractory epilepsy: a randomized controlled trial. Clin Sci. 2014:CS20130518. https://doi.org/10.1042/CS20130518.
Rong P, Liu J, Wang L, Liu R, Fang J, Zhao J, et al. Effect of transcutaneous auricular vagus nerve stimulation on major depressive disorder: a nonrandomized controlled pilot study. J Affect Disord. 2016;195:172–9. https://doi.org/10.1016/j.jad.2016.02.031.
Rossi RC, Vanderlei LCM, Gonçalves ACCR, Vanderlei FM, Bernardo AFB, Yamada KMH, et al. Impact of obesity on autonomic modulation, heart rate and blood pressure in obese young people. Auton Neurosci. 2015;193:138–41. https://doi.org/10.1016/j.autneu.2015.07.424.
Russo MA, Santarelli DM, O’Rourke D. The physiological effects of slow breathing in the healthy human. Breathe (Sheff). 2017;13(4):298–309. https://doi.org/10.1183/20734735.009817.
Salahuddin L, Cho J, Jeong MG, Kim D. Ultra short term analysis of heart rate variability for monitoring mental stress in mobile settings. Annu Int Conf IEEE Eng Med Biol Soc. 2007;2007:4656–9. https://doi.org/10.1109/IEMBS.2007.4353378.
Sandroni P, Benarroch EE, Low PA. Pharmacological dissection of components of the Valsalva maneuver in adrenergic failure. J Appl Physiol. 1991;71(4):1563–7. https://doi.org/10.1152/jappl.19184.108.40.2063.
Scheen AJ, Philips J-C. Squatting test: a dynamic postural manoeuvre to study baroreflex sensitivity. Clin Auton Res. 2012;22(1):35–41. https://doi.org/10.1007/s10286-011-0140-8.
Serhiyenko VA, Serhiyenko AA. Cardiac autonomic neuropathy: risk factors, diagnosis and treatment. World J Diabetes. 2018;9(1):1–24. https://doi.org/10.4239/wjd.v9.i1.1.
Shaffer F, Ginsberg JP. An overview of heart rate variability metrics and norms. Front Public Health. 2017;5. https://doi.org/10.3389/fpubh.2017.00258.
Sharon O, Fahoum F, Nir Y. Transcutaneous vagus nerve stimulation in humans induces pupil dilation and attenuates alpha oscillations. J Neurosci NJ-RM-1361-20. 2020;41(2):320–30. https://doi.org/10.1523/jneurosci.1361-20.2020.
Shen MJ, Zipes DP. Role of the autonomic nervous system in modulating cardiac arrhythmias. Circ Res. 2014;114(6):1004–21. https://doi.org/10.1161/CIRCRESAHA.113.302549.
Shields RW. Heart rate variability with deep breathing as a clinical test of cardiovagal function. CCJM. 2009;76(4 suppl 2):S37–40. https://doi.org/10.3949/ccjm.76.s2.08.
Shim HJ, Kwak MY, An Y-H, Kim DH, Kim YJ, Kim HJ. Feasibility and safety of transcutaneous Vagus nerve stimulation paired with notched music therapy for the treatment of chronic tinnitus. J Audiol Otol. 2015;19(3):159–67. https://doi.org/10.7874/jao.2015.19.3.159.
da Silva A, do Carmo J, Dubinion J, Hall JE. Role of sympathetic nervous system in obesity related hypertension. Curr Hypertens Rep. 2009;11:206.
Stavrakis S, Kulkarni K, Singh JP, Katritsis DG, Armoundas AA. Autonomic modulation of cardiac arrhythmias: methods to assess treatment and outcomes. JACC: Clinical Electrophysiology. 2020;6(5):467–83. https://doi.org/10.1016/j.jacep.2020.02.014.
Stefan H, Kreiselmeyer G, Kerling F, Kurzbuch K, Rauch C, Heers M, et al. Transcutaneous vagus nerve stimulation (t-VNS) in pharmacoresistant epilepsies: a proof of concept trial. Medline Abbreviated Title: Epilepsia. 2012;53(7):e115–8. https://doi.org/10.1111/j.1528-1167.2012.03492.x.
Stein PK, Kleiger RE, Rottman JN. Differing effects of age on heart rate variability in men and women. Am J Cardiol. 1997;80(3):302–5. https://doi.org/10.1016/S0002-9149(97)00350-0.
Sumner CJ, Sheth S, Griffin JW, Cornblath DR, Polydefkis M. The spectrum of neuropathy in diabetes and impaired glucose tolerance. Neurology. 2003;60(1):108–11. https://doi.org/10.1212/wnl.60.1.108.
Taylor A, McCartney N, Kamath M, Wiley R. Isometric training lowers resting blood pressure and modulates autonomic control. Med Sci Sports Exerc. 2003;35(2):251–6. https://doi.org/10.1249/01.mss.0000048725.15026.b5.
Thayer JF, Yamamoto SS, Brosschot JF. The relationship of autonomic imbalance, heart rate variability and cardiovascular disease risk factors. Int J Cardiol. 2010;141(2):122–31. https://doi.org/10.1016/j.ijcard.2009.09.543.
Tornatzky W, Miczek KA. Long-term impairment of autonomic circadian rhythms after brief intermittent social stress. Physiol Behav. 1993;53(5):983–93. https://doi.org/10.1016/0031-9384(93)90278-N.
Umetani K, Singer DH, McCraty R, Atkinson M. Twenty-four hour time domain heart rate variability and heart rate: relations to age and gender over nine decades. J Am Coll Cardiol. 1998;31(3):593–601. https://doi.org/10.1016/S0735-1097(97)00554-8.
Van de Borne P, Nguyen H, Biston P, Linkowski P, Degaute JP. Effects of wake and sleep stages on the 24-h autonomic control of blood pressure and heart rate in recumbent men. Am J Phys Heart Circ Phys. 1994;266(2):H548–54. https://doi.org/10.1152/ajpheart.1994.266.2.H548.
Verrotti A, Prezioso G, Scattoni R, Chiarelli F. Autonomic neuropathy in diabetes mellitus. Front Endocrinol (Lausanne). 2014;5. https://doi.org/10.3389/fendo.2014.00205.
Victor RG, Leimbach WN, Seals DR, Wallin BG, Mark AL. Effects of the cold pressor test on muscle sympathetic nerve activity in humans. Hypertension. 1987;9(5):429–36. https://doi.org/10.1161/01.HYP.9.5.429.
Vinik AI, Erbas T, Casellini CM. Diabetic cardiac autonomic neuropathy, inflammation and cardiovascular disease. J Diabetes Investig. 2013;4(1):4–18. https://doi.org/10.1111/jdi.12042.
Vogel ER, Sandroni P, Low PA. Blood pressure recovery from CME Valsalva maneuver in patients with autonomic failure. Neurology. 2005;65(10):1533–26. https://doi.org/10.1212/01.wnl.0000184504.13173.ef.
Weimer LH. Autonomic testing: common techniques and clinical applications. Neurologist. 2010;16(4):215–22. https://doi.org/10.1097/NRL.0b013e3181cf86ab.
Wiens J, Shenoy ES. Machine learning for healthcare: on the verge of a major shift in healthcare epidemiology. Clin Infect Dis. 2018;66(1):149–53. https://doi.org/10.1093/cid/cix731.
Wirch JL, Wolfe LA, Weissgerber TL, Davies GAL. Cold pressor test protocol to evaluate cardiac autonomic function. Appl Physiol Nutr Metab. 31:235–43. https://doi.org/10.1139/h05-018.
Yamamoto Y, Hughson RL, Peterson JC. Autonomic control of heart rate during exercise studied by heart rate variability spectral analysis. J Appl Physiol. 1991;71(3):1136–42. https://doi.org/10.1152/jappl.19220.127.116.116.
Yeh C-P, Huang H-C, Chang Y, Chen M-D, Hsu M. The Reliability and Validity of a Modified Squat Test to Predict Cardiopulmonary Fitness in Healthy Older Men. BioMed Research International 2018 Hindawi: e4863454, 2018.
Yoshimura N, White G, Weight FF, de Groat WC. Patch-clamp recordings from subpopulations of autonomic and afferent neurons identified by axonal tracing techniques. J Auton Nerv Syst. 1994;49(1):85–92. https://doi.org/10.1016/0165-1838(94)90024-8.
Young HA, Benton D. Heart-rate variability: a biomarker to study the influence of nutrition on physiological and psychological health? Behav Pharmacol. 2018;29(2 and 3):140–51. https://doi.org/10.1097/FBP.0000000000000383.
Zanos TP. Recording and decoding of vagal neural signals related to changes in physiological parameters and biomarkers of disease. Cold Spring Harb Perspect Med. 2019;9(12):a034157. https://doi.org/10.1101/cshperspect.a034157.
Zanos TP, Silverman HA, Levy T, Tsaava T, Battinelli E, Lorraine PW, et al. Identification of cytokine-specific sensory neural signals by decoding murine vagus nerve activity. PNAS. 2018;115(21):E4843–52. https://doi.org/10.1073/pnas.1719083115.
Zhang J. Effect of age and sex on heart rate variability in healthy subjects. J Manip Physiol Ther. 2007;30(5):374–9. https://doi.org/10.1016/j.jmpt.2007.04.001.
Ziegler D, Laux G, Dannehl K, Spüler M, Mühlen H, Mayer P, et al. Assessment of Cardiovascular Autonomic Function: Age-related Normal Ranges and Reproducibility of Spectral Analysis, Vector Analysis, and Standard Tests of Heart Rate Variation and Blood Pressure Responses. Diabet Med. 1992;9:166–75. https://doi.org/10.1111/j.1464-5491.1992.tb01754.x.
The authors would like to thank all participants who took part in testing sessions.
The research was supported by internal funding from the Feinstein Institutes for Medical Research and Northwell Health.
Ethics approval and consent to participate
This study was approved by the Northwell Health Institutional Review Board, IRB #19–0461. All participants were given detailed information about the research study, and completed a Consent for Participation in a Research Study form.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1 Supplementary Fig. 1. Normalizing pupillometry data.
Gaze location was used to calibrate the pupil diameter data. The RGB values of the fixation point in each frame were converted to relative luminance (L = 0.2126R + 0.7152G + 0.0722B, based on the luminosity function) and then linearly fit to raw pupil diameter measurements. By dividing by the slope of the line, the effects of brightness were minimized. Additionally, all sessions were completed in a quiet room with lowered ambient light.
Additional file 2 Supplementary Fig. 2. Average electrodermal activity and pupil dilation during autonomic testing (76 sessions).
The individual calculated responses (gray lines) were accumulated and averaged (black line) to extract an average response for each modality during each test. Each column represents a different calculated signal (electrodermal activity and pupil dilation). The dotted black traces correspond to a 95% confidence interval. RMSSD was not calculated for the cold pressor test and Valsalva maneuver due to time constraints necessary to accurately convey heart rate variability. Squat Test: vertical lines reflect changes in posture from standing to squatting and then squatting to standing. Cold Pressor Test: the first vertical line reflects when the participant immersed their hand into the ice water. The second vertical line represents the maximum of three minutes. The average trace only represents the individual traces available at that time point, as participants removed their hand at their own discretion. All participants kept their hand in the ice water for at least 30 s. Deep Breathing: vertical lines reflect when the deep breathing rate (6 breaths/minute) began and ended. In the third column for the RMSSD, only the first two minutes were analyzed. Diving Reflex: vertical lines reflect when the refrigerated gel-mask was placed on and removed from the participant’s face. A five second removal period is designated before the one minute of recovery. Valsalva Maneuver: vertical lines reflect the phases of the effort, from baseline, five seconds designated for inhalation and preparation, 15 s of the Valsalva maneuver, 10 s at the end of maneuver, and a final minute of recovery.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Debnath, S., Levy, T.J., Bellehsen, M. et al. A method to quantify autonomic nervous system function in healthy, able-bodied individuals. Bioelectron Med 7, 13 (2021). https://doi.org/10.1186/s42234-021-00075-7