Measuring Autonomic Symptoms with the Body Perception Questionnaire Short Form (BPQ-SF): Factor Analysis, Derivation of U.S. Adult Normative Values, and Association with Sensor-Based Physiological Measures ============================================================================================================================================================================================================= * Jacek Kolacz * Xiwei Chen * Evan J. Nix * Olivia K. Roath * Logan G. Holmes * Clarissa Tokash * Stephen W. Porges * Gregory F. Lewis ## Abstract **Objective** Autonomic regulation of organ and tissues may give rise to disruptions of typical functions. The Body Perception Questionnaire Short Form (BPQ-SF) includes items that were developed to assess autonomic symptoms in daily life. This pair of studies aimed to establish previously unexplored psychometric properties of the BPQ-SF, develop normative values for clinical and research use, and validate the self-reports with sensor-based measures. **Methods** Study 1 reports exploratory and confirmatory factor analysis (CFA) on BPQ-SF autonomic reactivity items from a large U.S. population-based online study (n = 2048). In study 2, BPQ-SF scores were examined for associations with heart period, respiratory sinus arrhythmia (RSA), and skin conductance (SC) during seated leg lifts in a community sample (n = 62). **Results** Study 1 results supported a 2-factor supra- and sub-diaphragmatic autonomic symptom solution (CFA: *RMSEA* = .040, *CFI* = .99, *TLI* = .99), though a 1-factor solution also fit the data well (*RMSEA* = .080, *CFI* = .99, *TLI* = .99). In study 2, flexible HP responses to lifts and rests were demonstrated at all autonomic symptom levels. However, low self-reported autonomic symptoms were associated with flexible dynamic RSA and SC, moderate symptoms with prolonged SC responses during rest periods, and high symptoms with little systematic changes in RSA and SC during leg lifts. **Conclusion** Results support the validity of self-reports of autonomic symptoms in research and clinical applications, with higher symptoms likely indicating impairment in autonomic flexibility. KEYWORDS * autonomic nervous system * autonomic symptoms * Body Perception Questionnaire * psychometrics * respiratory sinus arrhythmia * skin conductance ## Introduction The autonomic nervous system (ANS), which is responsive to emotional and metabolic challenges, is important for physical health and mental wellbeing. The pathways of the ANS and the brainstem regions with which it interfaces form a rapidly-acting brain-body interface that maintains functional regulation of organs and tissues, integrates information about the state of the body and environment, and controls responses to perceived and anticipated challenges and threats. (1–6) This system provides a physiological component of emotional arousal and regulation, (7–9) and is functionally integrated with brain regions that are involved in safety and threat appraisal (10) as well as the muscles of the face and head that are involved with social signaling of emotions. (11, 12) Its role in coordinating brain and body makes the ANS important for mental and physical health, with disrupted autonomic function being documented in PTSD, anxiety, autism, disorders of the gut, and chronic pain conditions such as fibromyalgia. (13–18) Furthermore, experimental evidence shows that alterations in caregiving environments can promote long-term changes in autonomic regulation (19–21), demonstrating how regulation of autonomic state is sensitive to social tuning and supporting the use of ANS state as an intervention target. The research and clinical importance of the ANS has led to a proliferation of sensor-based tools that quantify its function. (22) A complimentary and scalable source of ANS information can be acquired using self-reports, which document individuals’ experiences of the body’s function. Self-reports have valuable features as part of an assessment battery. First, self-reports are free or inexpensive to administer, and easy to implement at large scale for research, clinical monitoring, telehealth, and home use. In addition, measures of self-reported autonomic symptoms can contribute information about a range of ANS-innervated targets, compared to sensor-based measures which typically provide targeted information that quantifies performance of a single target organ. Although the function or dysfunction of any one organ can have a unique cause [e.g., tissue damage or a precise response to an environmental condition; (5)], combining information across multiple ANS-innervated targets provides a profile of underlying autonomic regulation. Despite these benefits, self-reports with established psychometric properties and validation are lacking. ### The Body Perception Questionnaire (BPQ-SF) The BPQ-SF (23; 24) is a self-report questionnaire developed to quantify experiences of body functions. The questionnaire includes autonomic reactivity subscales, which measure the disruption of organs innervated by the ANS. Items indexing autonomic reactivity are informed by the polyvagal theory (11,12,25), which posits that the brainstem structures regulating the ANS are integrated with the neural regulation of muscles of the face, neck, and head through the ventral vagal complex (VVC). Efferent pathways from the VVC innervate the heart, bronchi, esophagus, salivary, and lachrymal glands via visceromotor pathways, which are classically labeled as *autonomic* following Langley’s taxonomy. (1) However, outflow from the nucleus ambiguus of the VVC also innervates the laryngeal and pharyngeal muscles that are involved in swallowing and vocalization via somatomotor pathways, in addition to other muscles of the neck and face. (11) Informed by this neuroanatomical model of brainstem integration centers, the autonomic reactivity items of the BPQ include functions that are coordinated via visceromotor pathways (e.g., “I feel shortness of breath”) and those that include somatomotor pathways (e.g., “I have difficulty coordinating breathing and talking”). #### Dimensionality Data from 3 independent samples have supported two autonomic reactivity BPQ factors, corresponding to supra- and sub-diaphragmatic symptoms (above and below the diaphragm; Figure 1), with higher scores indicating greater autonomic disruption of typical organ function (24). This factor structure parallels the innervation patterns of the vagus nerve, which has outflow to structures above the diaphragm that emerge from the nucleus ambiguus (NA) and below the diaphragm from the dorsal motor nucleus of the vagus (DMNX; though both innervate the heart and bronchi). Subsequent studies with Italian and Chinese translations replicated this general factor structure. (26–28) ![Figure 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/04/28/2022.04.27.22274391/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2022/04/28/2022.04.27.22274391/F1) Figure 1. The Body Perception Questionnaire-Short Form (BPQ-SF) autonomic reactivity subscales are composed of items that measure experience of disruptions in circuits that affect function above and below the diaphragm. Because the BPQ has an ordered categorical response scale (5-point Likert-type scale ranging from “never” to “always”), analysis of item level data requires ordinal methods to avoid statistical distortions. (29) Sample size limitations of initial studies and subsequent factor analysis replications have resulted in the full ordinal outcomes becoming truncated to binary or 3-level responses to avoid overfitting the data influence of random noise. (24,26,27) These strategies have resulted in replicable, interpretable factor solutions but it is unclear whether the factor structure applies to full-item distributions, which are commonly used for subscale scoring. Thus, factor analysis on the English BPQ-SF with the full item distributions is needed. #### Norms Prior English-language studies have relied on convenience samples, limiting the generalizability of resulting score descriptive statistics. Normed values from a population-based sample are needed to inform score interpretation and clinical application. #### Validity The use of the BPQ to measure body awareness has been detailed in other studies, showing that BPQ scores are correlated with activity in brain regions associated with interoception (30) and research has suggested that the awareness subscale is a measure of interoceptive sensibility (31) or awareness of neutral and uncomfortable body sensations. (32) Associations of the body awareness subscale and other similar measures are explored in other studies. (24,27,32,33) However, to our knowledge, there has up to now been no systematic research to explore the association of the BPQ ANS reactivity subscales with sensor-based measures of physiological activity. ### The Present Studies Study 1 utilized a large national U.S. sample to examine the full-item response factor structure of the English-language Body Perception Questionnaire Short Form (BPQ-SF) and applied the results to develop population-based percentile and T scores. We hypothesized that factor analysis would converge on a structure similar to that which emerged from binary or three-level simplification of item response scales in past studies. Study 2 examined the associations of BPQ-SF autonomic reactivity subscales with sensor-based autonomic measures during seated leg raises. Electrocardiogram (ECG) and electrodermal activity (EDA) were used to assess: 1) changes in mean heart period, a broad metric of cardiac output, 2) respiratory sinus arrhythmia (RSA) amplitude, also known as high frequency heart rate variability, a measure of myelinated parasympathetic regulation, (34) 3) skin conductance reactivity frequency [an index of sympathetic regulation; (35)], and 4) vagal efficiency, a measure of the coordination of simultaneous changes in parasympathetic activity and heart period. (36, 37) According to principles of allostasis, (6, 38) flexible regulation of physiological activity to metabolic demands are key to maintaining adaptive responding, while less dynamic capability place an individual at risk of experiences of disrupted function. Given these principles, we hypothesized that higher BPQ autonomic reactivity scores would be related to reduced dynamic range of the ANS across sympathetic and parasympathetic metrics. We expected that high self-reported symptoms would be associated with less metabolic change in response to leg lifts, poorer ability to return to baseline during recovery periods (legs down), and lower vagal efficiency across the task period. Study 1. ## Method ### Procedure The study was approved by the Indiana University Institutional Review Board. All participants provided informed consent. Data were collected using an online survey consisting of the BPQ-SF, self-reported demographics, and several questionnaires not related to the present study. U.S. residents age 18 years and older were recruited through Qualtrics’ research panels, selected based on age, race, and gender to reflect the U.S. general population. Online panel data have been found to have similar internal reliability estimates and effect sizes between variables compared to conventional sampling methods (39) and U.S. residents recruited through Qualtrics are demographically similar to a national probability sample. (40) Recruitment and data quality check details are in Supplemental Materials. Factor analysis. Dimensionality was assessed by exploratory and confirmatory factor analysis (EFA and CFA) using the packages “psych” (41) and “lavaan” in R 4.1.1. (42) Factor analyses were applied to full item distributions using polychoric correlations and a weighted least squares estimator. (43) Exploratory factor retention was guided by loading structure, parallel analysis, and model fit. (44) Parallel analysis was conducted to examine the number of eigenvalues that deviated from the randomly simulated and re-sampled data, (45) with factors within +/-1 of the deviation from random data considered as viable solutions. (46) Goodness of fit was evaluated using the root mean squared error of approximation [RMSEA; (47, 48)], the Tucker-Lewis Index [TLI; (49)]; and the Comparative Fit Index [CFI; (50)]. Adequate fit was evidenced by an RMSEA value near .06 or lower as well as CFI and TLI values near .95 or greater. (51) EFA results were subject to geomin oblique rotation, (52; 53) which can reproduce correlated and uncorrelated factor structures. (44) The EFA solution was then applied to the second subsample using CFA and assessed using the fit indices described above. Calculation of percentiles and T scores. Percentiles were calculated from sum scores of subscales based on the factor analysis result. T-Scores were derived by normalizing ranks to a mean of 50 and standard deviation of 10. Unobserved scores were interpolated with a monotonic Hermite spline function. (54) ## Results The study 1 cohort was composed of 2048 adults. Descriptive statistics are in Table 1. The mean age was 46.34 years (SD = 17.19, range: 18-95) and 50.6% were female. The racial and ethnic distribution of the sample was similar to the general population. (55) Education level was diverse (14.9% graduate or professional degree, 66.9% at least a college degree, and 28.8% high school diploma). Responses were randomized into exploratory and confirmatory subsamples (each n = 1024). View this table: [Table 1.](http://medrxiv.org/content/early/2022/04/28/2022.04.27.22274391/T1) Table 1. Sample descriptive statistics for Study 1. Continuous variables are summarized with means and standard deviations. Categorical variables are summarized by counts and percentages. Exploratory factor analysis. Parallel analysis and fit indices suggested 2-5 factors. Examination of loadings showed that solutions with greater than 2 autonomic reactivity factors violated simple structure, with multiple factors having strong loadings on the same items and some items lacking any substantial loadings (Supplemental Materials). Thus, the 2-factor solution, corresponding to supra- and sub-diaphragmatic reactivity, was retained for CFA (Table 2). There was moderately high correlation between the two factors (r = .79). View this table: [Table 2.](http://medrxiv.org/content/early/2022/04/28/2022.04.27.22274391/T2) Table 2. Standardized exploratory and confirmatory factor analysis (EFA and CFA) loadings for the autonomic reactivity subscales of the Body Perception Questionnaire Short Form (BPQ-SF). Confirmatory factor analysis. The 2-factor solution fit well (*RMSEA* = .044 [90% CI: .039, .049], *CFI* = .99, *TLI* = .99). CFA loadings were similar to those found in the EFA (Table 2). Due to high correlation between factors (*r* = .81), we tested a 1-factor autonomic reactivity solution, which approached good fit (RMSEA = .081 [90% CI: .076, .085], CFI = .99, TLI = .99). Derivation of Percentiles and T Scores. Tables of percentile rank and T scores are in Supplemental Materials. Correlation between supra- and sub-diaphragmatic T scores was high (*r* = .73). Study 2. ## Method ### Procedure The study was approved by the Indiana University Institutional Review Board. Participants were recruited via flyers and screened for eligibility by phone and email. Exclusion criteria were: (1) medication use for heart conditions or blood pressure, (2) a pacemaker, (3) hospitalization for a head injury or major surgery in the past year, and (4) skin allergies. All participants provided informed consent and completed a demographic survey and the BPQ-SF on a laptop. They were then prepared for physiological data recording of electrocardiogram (ECG) with Lead II configuration and electrodermal activity (EDA) with electrodes applied thenar and hypothenar eminences of the right palm. Data were recorded at 1 kHz using AcqKnowledge Software 4.2 and the Biopac MP160 Data Acquisition System (BIOPAC Systems, Inc., Camino Goleta, CA, USA). For equipment and protocol details see Supplemental Materials. Seated straight-leg raises (SLR) were used to stimulate an autonomic response. This task was selected for its ease of execution, reliable demand on metabolic load, and low probability of leg movement disturbing sensors on the upper body and hands. First, participants sat quietly for a 2-minute baseline. Participants were then cued through 3 leg lifts during which they extended their legs for 30 seconds with shins parallel to the floor, followed by a 60-second rest with feet down (Figure 2, top panel). Research assistants marked timings at which legs were completely lifted and lowered. All participants were able to keep legs elevated for the full length of all 3 lifts. ![Figure 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/04/28/2022.04.27.22274391/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2022/04/28/2022.04.27.22274391/F2) Figure 2. Autonomic activity changes during seated leg lifts as a function of level of autonomic symptoms self-reported on the Body Perception Questionnaire Short Form (BPQ-SF). Top row: Mean heart period (HP), respiratory sinus arrhythmia (RSA), and skin conductance reactivity (SCR) rate were measured during seated baseline, leg lifts, and recovery (legs on floor). Middle row: Raw data and boxplots with between-group significance tests from repeated measures models of absolute values at baseline, leg lifts, and recovery periods. For within-group contrasts see Supplemental Materials. HP and RSA models are adjusted for age and BMI. Bottom row: Linear mixed model predicted values and 95% confidence intervals for HP, RSA, and SCR rate change scores from baseline to leg lifts, and from baseline to recovery periods as a function of self-reported autonomic symptoms. Dotted horizontal lines designate no change. ### Data cleaning and processing Physiological data cleaning and processing was conducted using AcqKnowledge 4.2 (BIOPAC Systems, Inc., Camino Goleta, CA, USA) and Matlab R2020a. (56) Inter-beat intervals (IBI) were measured by the timing between R-wave peaks from ECG data. IBIs were manually inspected for artifacts, arrhythmias, and missed beats, which may lead to bias and invalidation of beat-to-beat cardiac measures. (57) Missed beats were reconstructed, noise artifacts corrected, and ventricular arrhythmia fluctuations were removed to produce a clean, uninterrupted beat-to-beat signal using CardioEdit+ Software (Brain-Body Center for Psychophysiology and Engineering, University of North Carolina at Chapel Hill). Cleaned data were analyzed in 15-sec epochs. The first and last epochs of recovery periods were excluded from analysis to avoid transition periods and anticipatory autonomic activation immediately before and after leg lifts. #### Heart period (HP) HP was quantified as the mean time between successive heart beats in milliseconds, reflecting the sum of all influences on heart beat timing including sympathetic and parasympathetic effects. #### Respiratory sinus arrhythmia (RSA) amplitude RSA amplitude is an index of cardiac vagal tone, the effect of the outflow of the nucleus ambiguus via the myelinated vagus nerve pathway to the heart. (34) Calculations were conducted with the Porges-Bohrer method, (58, 59) using a 51-point moving polynomial filter and a second finite impulse response (FIR) type bandpass filter adapted for a 5-Hz resampling rate of the edited IBI time-series. Variability amplitude was quantified in the .12 - 1.0 Hz band, based on the respiratory frequency associated with spontaneous breathing in adults (.12 - .40 Hz) and expanded to include potentially faster respiration rates during the leg-up condition. (34) The result was natural log transformed to reduce skew. #### Vagal efficiency (VE) VE indexes the extent to which changes in cardiac vagal tone are associated with changes in cardiac output (36,37,60), calculated from within-individual slopes of synchronous HP and RSA values. The resulting metric describes the observed change in HP (ms) for a one-unit (ln(ms2)) change in RSA. #### Skin conductance response (SCR) frequency SCRs calculated from the raw EDA signal were used to index sympathetic activation. SCRs are caused by sympathetic activation that promotes the opening of palmar eccrine sweat pores, increasing skin surface conductivity. (35, 61) EDA data were downsampled to 7.812 Hz for analysis and visually inspected. Segments of data with signal noise or disruption were interpolated or excluded from analysis. Phasic data, which were separated from tonic data using smoothing baseline removal, were used to locate SCRs via trapezoidal numeric integration (minimum amplitude criterion = 03 µS) and summed over time to reflect firing rate per minute. ### Data analysis Statistical analyses were conducted using R 4.1 (42) and RStudio 1.4. (62) Physiological change scores were calculated by subtracting lift and recovery values from baseline. BPQ autonomic reactivity scores were converted to T scores and tertiles based on general population percentiles from Study 1. BMI was calculated from self-reported weight and height (kg/m2). Group differences were evaluated using Chi-Square and Welch’s unequal variance tests. (63) Change scores and the effects of seated leg lifts on HP and RSA were examined with linear mixed models. Standardized mean difference (Cohen’s d) was calculated for pairwise contrasts. (64) Due to skew and floor effects in SCR data, tobit regression and nonparametric repeated measures model were used. (65) Post hoc pairwise comparisons for non-parametric models were conducted using the two-tailed Mann-Whitney U test for between group comparisons and Wilcoxon signed-rank tests for within group comparisons. Rank-biserial correlations were calculated as an effect size for pairwise differences. (66) All statistical tests had alpha set at .05 and were two-sided. Based on their know associations, HP and RSA models were adjusted for age and BMI (75; associations between physiological variables and covariates are in Supplemental Materials). ## Results The sample was composed of 62 participants. Descriptive statistics are in Table 3. The mean age was 21.85 years (SD = 4.01, range: 18-36), 67.7% were female, and 48.4% had a bachelor’s degree or higher level of education. The mean BMI was 23.84 (SD = 4.53). View this table: [Table 3.](http://medrxiv.org/content/early/2022/04/28/2022.04.27.22274391/T3) Table 3. Sample descriptive statistics for Study 2. Continuous variables are summarized with means and standard deviations. Categorical variables are summarized by counts and percentages. Supra- and sub-diaphragmatic scores on the BPQ-SF were positively correlated (*r* = .59) and analysis results were comparable for the two subscales. Thus, results are reported on combined autonomic symptom scores (see Supplemental Materials for supra- and sub-diaphragmatic results reported separately). BPQ-SF autonomic symptom scores from the lab study were lower compared to the general population in Study 1 (low tertile n = 29, 46.8%, mid tertile n = 22, 35.5%, high tertile n = 11, 17.7%). BPQ autonomic symptom scores were not associated with age, gender, body mass index, or meeting overweight and obesity criteria (all p > .05). After data inspection, 9 EDA files were excluded due to signal drop out, poor quality, or equipment errors (final n = 53). ### Associations among physiological variables Correlations were calculated across the 7 postures (baseline, 3 sitting, 3 recovery). HP and RSA were moderately correlated (Pearson r = .49, p < .001). SCR rate was not significantly correlated with HP (Rho = -.17, p = .226) or RSA (Rho = -.06, p = .68). ### Associations between physiological variables and self-reported autonomic symptoms BPQ autonomic symptom groups did not differ significantly on HP, RSA, or SCR rate at baseline (all contrasts p > .10). BPQ scores were negatively associated with lower SCR rate during leg lifts (B = -0.24, SE = 0.083, z = -2.94, p = 0.003) and recovery periods (B = - 0.16, SE = 0.05, z = -2.99, p = 0.003), but no other linear associations were significant. All BPQ autonomic symptom tertile groups reduced mean heart period from baseline to leg lifts (Low = -57.21 [95% CI: -68.04, -46.38], d = 2.41; Mid = -64.53 [95% CI: -76.96, - 52.10], d = 2.72; High = -48.91 [95% CI: -66.49, -31.33], d = 2.06) and there were no significant differences among them during leg lifts (all contrasts p > .30). Only the mid autonomic symptom group maintained shorter mean heart period during recovery conditions compared to baseline (Figure 2 bottom left). This effect was significantly different from the low symptom group, in which the recovery periods were not distinguishable from baseline. Only the low and mid symptom tertile groups significantly decreased RSA in response to the leg lift task (Low = -.43 [95% CI: -.63, -.23], d = 1.12; Mid = -.53 [95% CI: -.76, -.31], d = 1.37). Both groups returned to baseline levels during recovery periods (Baseline-Recovery Change Low = -.10 [95% CI: -.29, .10], d = 0.25; Mid = -.07 [95% CI: -.29, .16], d = 0.17). The high symptom group did not show significant RSA change from baseline to leg lifts (B = -.15 [95% CI: -.47, .17], d = 0.39) or from baseline to recovery periods (B = -.14 [95% CI: -.18, .46], d = 0.37). SCR rate increased with leg lifts only in the low and mid groups (Low = 2.47 [95% CI: .82, 4.12], r = -0.60; Mid = 2.50 [95% CI: .65, 4.35], r = -0.71). During recovery periods, the SCR rate significantly declined from leg lifts in the low group (B = 1.70 [95% CI: .69, 2.71], r = 0.58) and was not significantly different from baseline. In contrast, the mid group maintained skin conductance activation above baseline values during recovery periods (Baseline-Recovery B = 2.19 [95% CI: .34, 4.03], r = -0.73). The high symptom group did not significantly differ from baseline during leg lifts or recovery periods. Overall, the high symptom group had a lower SCR rate during leg lifts and recovery compared to the low symptom group (leg lift p = .007; r = .64; recovery p = .002; r = .72), and the mid group was higher than the high group at recovery (p = .002; r = .77). VE was not associated with BPQ autonomic symptom levels (F (2, 53) = .83, p = .442). ## Discussion This is the first study to conduct factor analysis on full item responses of the English language Body Perception Questionnaire (BPQ-SF), derive normative values for research and clinical use, and assess the convergence of BPQ autonomic reactivity self-reports with sensor-based measures of parasympathetic and sympathetic function. The results demonstrate that 1) the autonomic reactivity subscales have a 1- or 2-factor structure when explored with full item response distributions in a large, general US population of adults; 2) that low levels of self-reported autonomic symptoms are associated with dynamic parasympathetic and sympathetic responses to a metabolic challenge, while high self-reported autonomic symptoms show less dynamically responsive patterns, and 3) that normative values from the general population study can be used as reference values for clinical studies and provide meaningful interpretation of autonomic symptoms. In study 1, the results of exploratory and confirmatory factor analysis in a large, general population sample converged with prior factor analysis that had been conducted using simplified item response ranges in studies with convenience and student samples. The results revealed two autonomic symptom factors that differentiated supra- and sub-diaphragmatic reactivity. Also consistent with prior studies, the supra- and sub-diaphragmatic subscales were strongly positively correlated (24, 26), suggesting that experiences of autonomic symptoms above and below the diagram are closely associated. In addition, supra- and sub-diaphragmatic reactivity subscale patterns showed similar associations with sensor-based physiological measures in lab data in Study 2. Prior studies have also found patterns of parallel associations between supra- and sub-diaphragmatic symptoms with trauma history, PTSD symptoms, and depression symptoms. (67, 68) This suggests that under some conditions supra- and sub-diaphragmatic reactivity may be scored as a single subscale, though distinctions between subscales should continue to be explored. It also underscores that processes that regulate organs above and below the diagram are functionally integrated and likely impacted in tandem. Consistent with previous studies, we found that the urge to vomit is associated with both supra- and sub-diaphragmatic reactivity factors (24, 26), likely reflecting the neuro-anatomical circuits that induce vomiting, which involve structures both above and below the diaphragm. (45, 69) Vomiting requires a coordinated outflow from both the dorsal motor nucleus of the vagus (which mainly innervates structures below the diaphragm) and the nucleus ambiguus (regulating structures above the diaphragm) which synchronize swallowing, respiration, cardiovascular activity, and gastrointestinal function. (69) Disruption at any of these pathways may result in sensitization towards an urge to vomit. In study two, participants in all autonomic symptom level groups responded metabolically to leg lifts with shortened average heart period. However, parasympathetic and sympathetic metrics showed different patterns of dynamic autonomic regulation associated with self-reported autonomic symptoms. Those with low autonomic symptoms had the most flexible autonomic regulation, withdrawing PNS and activating SNS during leg lifts, then increasing PNS while reducing SNS rapidly between challenges. This pattern reflects maximum flexibility in responding to momentary metabolic needs, a hallmark of efficient allostasis. (6, 38) This dynamic flexibility may protect individuals from extended or uncontrolled activation that lasts beyond immediate metabolic needs which could lead to experiencing autonomic symptoms. Similarly, moderate self-reported autonomic symptom levels were also associated with efficient withdrawal and re-engagement of the parasympathetic vagal brake. However, while this group exhibited sympathetic activation during leg lifts compared to baseline, this sympathetic activation was maintained even during rest periods. This was also reflected in average heart period change, where HP was shortened even during rest periods. Lack of dynamic flexibility in reducing sympathetic responses may cause some disruption in organ function over time, resulting in more experiences of autonomic symptoms. The highest levels of self-reported autonomic symptoms were associated with no consistent parasympathetic or sympathetic change. Individuals at these levels may be least organized in their autonomic response, not predictably engaging or disengaging the ventral vagal parasympathetic brake or the SNS, placing them most at risk of long-term disruptions in the function of organs and tissues innervated by the autonomic nervous system. However, this group nonetheless mounted a change in overall heart period to leg lifts and a return toward baseline during rest periods, a pattern that requires some neural regulation of the heart. The polyvagal theory posits that metabolic regulation via the unmyelinated cardioinhibitory vagal pathway to the heart via the dorsal vagal motor nucleus of the vagus (DMNX) would be associated with the highest levels of autonomic symptoms due to its being the evolutionarily oldest pathway and its intimate integration with the gut. (25) Electrophysiological studies have documented that direct stimulation of the DMNX will produce heart rate slowing (70–72) and that vagal activation influences sympathetic outflow (e.g., (73)). However, there is no known non-invasive marker of the cardioinhibitory influence of the DMNX and the data in this study cannot confirm whether it is involved in cardiac regulation during leg lifts in the high symptom group. VE was not associated with autonomic symptoms in this study. Prior studies have quantified VE during posture shifts that engage the baroreflex maintaining blood pressure to the brain (37, 60) or strenuous exercise on a stationary bicycle that requires strenuous metabolic muscle activation. (74) Compared to these studies, leg lifts may be a relatively weak ANS challenge, not able to provide enough metabolic engagement to permit robust VE measurement. Future studies that use paradigms that provide stronger challenge are needed. Viewed in sum, these studies support the use of self-reported autonomic symptoms via the BPQ as an index of autonomic disruption of typical organ and tissue function. However, there are important distinctions between sensor-based measures, which are typically used to index autonomic activity, and self-reports. First, self-reports and sensor-based measures differ on their ultimate goals. Whereas sensors typically focus on the dynamic function of a single organ or tissue, self-reports pool information over many innervation targets. Also, self-reported symptoms are likely most noticeable when their function is disrupted and are likely less sensitive to slight fluctuations in dynamic neural regulation. In addition, self-reports concatenate experience over long periods of time to provide enough opportunity to observe potential disruptions. Finally, compared to sensor-based measures, questionnaires have measurement limitations that include reliance on memory, reporting biases, attention, and other cognitive factors. However, the benefits of self-report measures in providing an overall summary of autonomic disruption at a low burden for participants makes them a potentially useful part of an assessment battery. Though the BPQ may index disruptions in regulation of the autonomic nervous system, the established factor structure is not specific to a particular autonomic branch. While some organs and tissues are innervated by a single pathway, such as the sympathetic regulation of sweat glands in the palm, most autonomic innervation targets are regulated by multiple branches. (5, 11) These multiple innervation paths, along with the co-regulation of the various branches of the ANS by higher level structures such as the hypothalamus, (5) make distinctions between specific autonomic pathways challenging for self-report measurement. BPQ-SF autonomic symptoms consistently cluster together across samples without specificity to one branch or another. The findings reported here, however, do support the use of the BPQ-SF autonomic symptom items as a measure of allostatic flexibility or autonomic resilience. High symptoms are associated with weaker flexibility to mount responses to challenge and recover after metabolic resources are not needed. Clinical and research use of self-report symptoms can be supported by the derivation of normative values and generation of cut scores. This paper describes the first derivation of normative percentile and T scores from a large general US sample. While study 2 used a prior division into tertiles, showing non-linear distinctions between low, moderate, and high symptom groupings, larger samples are needed to optimize cut scores to capture autonomic disruption with optimized sensitivity and specificity. ### Limitations Study 2 was based on a community convenience sample. We found that the sample had restricted age and BMI distribution, compared to a general population. Both factors have a known association with patterns of autonomic function. (75) This sample also had a larger proportion of low autonomic symptom scores and lower proportion of those at the highest levels compared to the online sample that was more representative of the range of BMI and age in a general population. Future studies need to examine the association with autonomic symptoms with sensor-based measures in more broad population-based samples. ## Conclusion These studies support the use of the BPQ-SF autonomic reactivity subscales and their validity to index autonomic disruption of organ functions. Assessing experiences of autonomic symptoms can be utilized in clinical monitoring of patient symptoms and in basic research into mechanisms that affect mental wellbeing and physical health. ## Supporting information Supplemental Material [[supplements/274391_file04.docx]](pending:yes) ## Data Availability All data produced in the present study are available upon reasonable request to the authors ## Conflicts of Interest and Source of Funding The authors have no financial conflicts of interest to declare. This work was supported by the Dillon Fund, the United States Association for Body Psychotherapy (USABP), and the Chaja Foundation. The funders had no role in the design, data collection, analysis, or manuscript preparation. ## Acknowledgements We wish to thank Juliet Dahlgren for assistance with physiological data processing and R. Mike Winters on an early version of physiological data analysis. ## Footnotes * **Author Contributions** The authors confirm contribution to the paper as follows. Study conceptualization and design: Jacek Kolacz and Greg Lewis. Data collection: Clarissa Tokash, Logan Holmes, Olivia Roath, Gregory Lewis, Jacek Kolacz. Data processing and analysis: Xiwei Chen, Evan Nix, Logan Holmes, Clarissa Tokash, Jacek Kolacz, and Greg Lewis. Data interpretation: Jacek Kolacz, Greg Lewis, and Stephen Porges. First draft of manuscript: Jacek Kolacz. Additional drafts of manuscript: All authors. All authors have reviewed the results and approved the final version. ACRONYMS : ANS : autonomic nervous system BMI : body mass index BPQ : Body Perception Questionnaire BPQ-SF : Body Perception Questionnaire Short Form CFA : confirmatory factor analysis CFI : Comparative Fit Index DMNX : dorsal motor nucleus of the vagus ECG : electrocardiogram EDA : electrodermal activity EFA : exploratory factor analysis FIR : finite impulse response HF-HRV : high frequency heart rate variability IBI : inter-beat intervals HP : heart period NA : nucleus ambiguus PNS : parasympathetic nervous system PTSD : post-traumatic stress disorder RMSEA : root mean squared error of approximation RSA : respiratory sinus arrythmia SCR : skin conductance reactivity SLR : seated straight-leg raise SNS : sympathetic nervous system TLI : Tucker-Lewis Index VE : vagal efficiency VVC : ventral vagal complex * Received April 27, 2022. * Revision received April 27, 2022. * Accepted April 28, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Langley JN. The autonomic nervous system. W. Heffer; 1921. 2. 2. W.B. Cannon. Organization for physiological homeostasis. Physiological Reviews. 9:399– 431. 3. 3.Hess WR. The Organisation of the Autonomic Nervous System. Basel, Switzerland: Benno Schwabe & Co.; 1948. 4. 4.Porges SW. Cardiac vagal tone: a physiological index of stress. Neuroscience & Biobehavioral Reviews. 1995;19:225–33. 5. 5.Jänig W. Integrative Action of the Autonomic Nervous System: Neurobiology of Homeostasis. Cambridge University Press; 2008. 6. 6.Sterling P. Allostasis: a model of predictive regulation. Physiology & behavior. 2012;106:5–15. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.physbeh.2011.06.004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21684297&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000301910000002&link_type=ISI) 7. 7.Porges SW, Doussard-Roosevelt JA, Maiti AK. Vagal tone and the physiological regulation of emotion. Monographs of the society for research in child development. 1994;167–86. 8. 8.Thayer JF, Lane RD. A model of neurovisceral integration in emotion regulation and dysregulation. Journal of Affective Disorders. 2000;61:201–16. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0165-0327(00)00338-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11163422&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000166745500006&link_type=ISI) 9. 9.MacCormack JK, Gaudier-Diaz MM, Armstrong-Carter EL, Arevalo JM, Meltzer-Brody S, Sloan EK, Cole SW, Muscatell KA. Beta-adrenergic blockade blunts inflammatory and antiviral/antibody gene expression responses to acute psychosocial stress. Neuropsychopharmacology. 2021;46:756–62. 10. 10.Thayer JF, Åhs F, Fredrikson M, Sollers III JJ, Wager TD. A meta-analysis of heart rate variability and neuroimaging studies: implications for heart rate variability as a marker of stress and health. Neuroscience & Biobehavioral Reviews. 2012;36:747–56. 11. 11.Porges SW. Orienting in a defensive world: Mammalian modifications of our evolutionary heritage. A Polyvagal Theory. Psychophysiology. 1995;32:301–18. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1469-8986.1995.tb01213.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7652107&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995RJ25200001&link_type=ISI) 12. 12.Porges SW. Cardiac vagal tone: a neurophysiological mechanism that evolved in mammals to dampen threat reactions and promote sociality. World Psychiatry. 2021;20:296. 13. 13.Kemp AH, Quintana DS, Gray MA, Felmingham KL, Brown K, Gatt JM. Impact of Depression and Antidepressant Treatment on Heart Rate Variability: A Review and Meta-Analysis. Biological Psychiatry. 2010;67:1067–74. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2009.12.012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20138254&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000277918000009&link_type=ISI) 14. 14.Chalmers JA, Quintana DS, Abbott MJ-A, Kemp AH. Anxiety Disorders are Associated with Reduced Heart Rate Variability: A Meta-Analysis. Frontiers in Psychiatry. 2014;5. 15. 15.Koenig J, Kemp AH, Beauchaine TP, Thayer JF, Kaess M. Depression and resting state heart rate variability in children and adolescents—A systematic review and meta-analysis. Clinical psychology review. 2016;46:136–50. 16. 16.Kolacz J, Porges SW. Chronic Diffuse Pain and Functional Gastrointestinal Disorders After Traumatic Stress: Pathophysiology Through a Polyvagal Perspective. Frontiers in Medicine. 2018;5. 17. 17.Kolacz J, Kovacic KK, Porges SW. Traumatic stress and the autonomic brain-gut connection in development: Polyvagal Theory as an integrative framework for psychosocial and gastrointestinal pathology. Developmental Psychobiology. 2019;61:796–809. 18. 18.Schneider M, Schwerdtfeger A. Autonomic dysfunction in posttraumatic stress disorder indexed by heart rate variability: a meta-analysis. Psychological medicine. 2020;50:1937– 48. 19. 19.McLaughlin KA, Sheridan MA, Tibu F, Fox NA, Zeanah CH, Nelson CA. Causal effects of the early caregiving environment on development of stress response systems in children. Proceedings of the National Academy of Sciences. 2015;112:5637–42. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTEyLzE4LzU2MzciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wNC8yOC8yMDIyLjA0LjI3LjIyMjc0MzkxLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 20. 20.Hastings PD, Kahle S, Fleming C, Lohr MJ, Katz LF, Oxford ML. An intervention that increases parental sensitivity in families referred to Child Protective Services also changes toddlers’ parasympathetic regulation. Developmental science. 2019;22:e12725. 21. 21.Porges SW, Davila MI, Lewis GF, Kolacz J, Okonmah Obazee S, Hane AA, Kwon KY, Ludwig RJ, Myers MM, Welch MG. Autonomic regulation of preterm infants is enhanced by Family Nurture Intervention. Developmental Psychobiology. 2019;61:942–52. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) 22. 22.Shaffer F, Ginsberg JP. An overview of heart rate variability metrics and norms. Frontiers in public health. 2017;258. 23. 23.Porges SW. Body Perception Questionnaire. Laboratory of Developmental Assessment, University of Maryland. 1993. 24. 24.Cabrera A, Kolacz J, Pailhez G, Bulbena Cabre A, Bulbena A, Porges SW. Assessing body awareness and autonomic reactivity: Factor structure and psychometric properties of the Body Perception Questionnaire-Short Form (BPQ-SF). International Journal of Methods in Psychiatric Research. 2018;27:e1596. 25. 25.Porges SW. The polyvagal theory: Neurophysiological foundations of emotions, attachment, communication, and self-regulation WW Norton & Company; 2011. 26. 26.Cerritelli F, Consorti G, D’Alessandro G, Galli M, Kolacz J, Porges SW. Cross-cultural adaptation and validity of the Italian version of the Body Perception Questionnaire. Under Review. 27. 27.Poli A, Maremmani AGI, Chiorri C, Mazzoni G-P, Orrù G, Kolacz J, Porges SW, Conversano C, Gemignani A, Miccoli M. Item Reduction, Psychometric and Biometric Properties of the Italian Version of the Body Perception Questionnaire—Short Form (BPQ-SF): The BPQ-22. International Journal of Environmental Research and Public Health. 2021;18:3835. 28. 28.Wang N, Ren F, Zhou X. Factor Structure and Psychometric Properties of the Body Perception Questionnaire–Short Form (BPQ-SF) Among Chinese College Students. Frontiers in Psychology. 2020;11:1355. 29. 29.Liddell TM, Kruschke JK. Analyzing ordinal data with metric models: What could possibly go wrong? Journal of Experimental Social Psychology. 2018;79:328–48. 30. 30.Critchley HD, Wiens S, Rotshtein P, Öhman A, Dolan RJ. Neural systems supporting interoceptive awareness. Nature neuroscience. 2004;7:189–95. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nn1176&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14730305&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000188468500018&link_type=ISI) 31. 31.Garfinkel SN, Seth AK, Barrett AB, Suzuki K, Critchley HD. Knowing your own heart: distinguishing interoceptive accuracy from interoceptive awareness. Biological psychology. 2015;104:65–74. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsycho.2014.11.004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25451381&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) 32. 32.Desmedt O, Heeren A, Corneille O, Luminet O. What Do Measures of Self-report Interoception Measure? Insights from A Systematic Review, Latent Factor Analysis, and Network Approach. 2021. 33. 33.Kobayashi R, Honda T, Machizawa M, Ichikawa N, Nakao T. Factor structure, reliability and validation of the Japanese version of the Body Perception Questionnaire-Body Awareness Very Short Form (BPQ-BAVSF-J). Japanese Journal of Research on Emotions. 2021;28:38–48. 34. 34.Bernston GG, Bigger Jr JT, Eckberg DL, Grossman P, Kaufman PG, Malik M, Nagaraja HN, Porges SW, Saul JP, Stone PH, van der Molen MW. Heart rate variability: origins, methods, and interpretive caveats. Psychophysiology. 1997;34:623–48. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1469-8986.1997.tb02140.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9401419&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997YH88700001&link_type=ISI) 35. 35.Dawson ME, Schell AM, Filion DL. The electrodermal system. In: Handbook of Physiology. 3rd ed. New York: Cambridge University Press; 2007. p. 159–81. 36. 36.Porges SW, Doussard-Roosevelt JA, Stifter CA, McCLENNY BD, Riniolo TC. Sleep state and vagal regulation of heart period patterns in the human newborn: An extension of the polyvagal theory. Psychophysiology. 1999;36:14–21. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S004857729997035X&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10098376&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000078095900002&link_type=ISI) 37. 37.Kolacz J, Kovacic K, Lewis GF, Sood MR, Aziz Q, Roath OR, Porges SW. Cardiac autonomic regulation and joint hypermobility in adolescents with functional abdominal pain disorders. Neurogastroenterology & Motility. 2021;33:e14165. 38. 38.Sterling P. Allostasis: a new paradigm to explain arousal pathology. Handbook of life stress, cognition and health. 1988. 39. 39.Walter SL, Seibert SE, Goering D, O’Boyle EH. A Tale of Two Sample Sources: Do Results from Online Panel Data and Conventional Data Converge? Journal of Business and Psychology. 2019;34:425–52. 40. 40.Boas TC, Christenson DP, Glick DM. Recruiting large online samples in the United States and India: Facebook, Mechanical Turk, and Qualtrics. Political Science Research and Methods. 2020;8:232–50. 41. 41.Revelle W. psych: Procedures for Personality and Psychological Research, Northwestern University. Version. 2020;1:12. 42. 42.R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2021. Available from: [https://www.R-project.org/](https://www.R-project.org/) 43. 43.Barendse MT, Oort FJ, Timmerman ME. Using exploratory factor analysis to determine the dimensionality of discrete responses. Structural Equation Modeling: A Multidisciplinary Journal. 2015;22:87–101. 44. 44.Fabrigar L, Wegener D t, MacCullum R. Strahan, e. J.(1999). evaluating the use of factor analysis in psychological research. Psychological Methods. 1999;4:272–99. 45. 45.Horn JL. A rationale and test for the number of components from the matrix of partial correlations. Psychometrika. 1965;30:179–85. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/BF02289447&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14306381&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1965CBT8600005&link_type=ISI) 46. 46.Lim S, Jahng S. Determining the number of factors using parallel analysis and its recent variants. Psychological Methods. 2019;24:452. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) 47. 47.Steiger JH, Lind JC. Statistically-based tests for the number of common factors. Annual meeting of the Psychometric Society; 1980 May; Iowa City, IA. 48. 48.Steiger JH. Structural model evaluation and modification: An interval estimation approach. Multivariate behavioral research. 1990;25:173–80. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1207/s15327906mbr2502_4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26794479&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990DZ95200004&link_type=ISI) 49. 49.Tucker LR, Lewis C. A reliability coefficient for maximum likelihood factor analysis. Psychometrika. 1973;38:1–10. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/BF02291170&link_type=DOI) 50. 50.Bentler PM. Comparative fit indexes in structural models. Psychological Bulletin. 1990;107:238. 51. 51.Hu L, Bentler PM. Cutoff criteria for fit indexes in covariance structure analysis: Conventional criteria versus new alternatives. Structural Equation Modeling: A Multidisciplinary Journal. 1999;6:1–55. 52. 52.Yates A. Multivariate exploratory data analysis: A perspective on exploratory factor analysis. Suny Press; 1987. 53. 53.Sass DA, Schmitt TA. A comparative investigation of rotation criteria within exploratory factor analysis. Multivariate Behavioral Research. 2010;45:73–103. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/00273170903504810&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26789085&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000276098300003&link_type=ISI) 54. 54.Dougherty RL, Edelman AS, Hyman JM. Nonnegativity-, monotonicity-, or convexity-preserving cubic and quintic Hermite interpolation. mathematics of computation. 1989;52:471–94. 55. 55.U.S. Census Bureau. 2020 Census. 2020. Available from: [https://www.census.gov/programs-surveys/decennial-census/decade/2020/2020-census-main.html](https://www.census.gov/programs-surveys/decennial-census/decade/2020/2020-census-main.html) 56. 56.Matlab R2020a. Natick, Massachusetts: The MathWorks Inc. 2020. 57. 57.Porges SW, Byrne EA. Research methods for measurement of heart rate and respiration. Biological Psychology. 1992;34:93–130. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0301-0511(92)90012-J&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1467397&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1992KA29300003&link_type=ISI) 58. 58.Porges SW. Method and apparatus for evaluating rhythmic oscillations in aperiodic physiological response systems. US4510944A, 1985 [cited 2020 Aug 4]. 59. 59.Lewis GF, Furman SA, McCool MF, Porges SW. Statistical strategies to quantify respiratory sinus arrhythmia: Are commonly used metrics equivalent? Biological Psychology. 2012;89:349–64. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsycho.2011.11.009&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22138367&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) 60. 60.Kovacic K, Kolacz J, Lewis GF, Porges SW. Impaired vagal efficiency predicts auricular neurostimulation response in adolescent functional abdominal pain disorders. Official journal of the American College of Gastroenterology| ACG. 2020;115:1534–38. 61. 61.Asahina M, Poudel A, Hirano S. Sweating on the palm and sole: physiological and clinical relevance. Clinical Autonomic Research. 2015;25:153–59. 62. 62.RStudio Team. RStudio: Integrated Development Environment for R. RStudio. PBC, Boston, MA; 2021. Available from: [http://www.rstudio.com/](http://www.rstudio.com/). 63. 63.Ruxton GD. The unequal variance t-test is an underused alternative to Student’s t-test and the Mann–Whitney U test. Behavioral Ecology. 2006;17:688–90. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/beheco/ark016&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000238450300026&link_type=ISI) 64. 64.Cohen J. Statistical power analysis for the behavioral sciences, 2nd edn. Á/L. Erbaum Press, Hillsdale, NJ, USA; 1988. 65. 65.Noguchi K, Gel YR, Brunner E, Konietschke F. nparLD: an R software package for the nonparametric analysis of longitudinal data in factorial experiments. Journal of Statistical software. 2012;50:1–23. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/jss.v050.i13&link_type=DOI) 66. 66.Cureton EE. Rank-biserial correlation. Psychometrika. 1956;21:287–90. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/BF02289138&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1956CFY3800007&link_type=ISI) 67. 67.Kolacz J, Dale LP, Nix EJ, Roath OK, Lewis GF, Porges SW. Adversity History Predicts Self-Reported Autonomic Reactivity and Mental Health in US Residents During the COVID-19 Pandemic. Frontiers in Psychiatry. 2020;11:1119. 68. 68.Kolacz J, Hu Y, Gesselman AN, Garcia JR, Lewis GF, Porges SW. Sexual function in adults with a history of childhood maltreatment: Mediating effects of self-reported autonomic reactivity. Psychological Trauma: Theory, Research, Practice, and Policy. 2020;12:281–90. 69. 69.Babic T, Browning KN. The role of vagal neurocircuits in the regulation of nausea and vomiting. European Journal of Pharmacology. 2014;722:38–47. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ejphar.2013.08.047&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24184670&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) 70. 70.Jordan D, Khalid MEM, Schneiderman N, Spyer KM. The location and properties of preganglionic vagal cardiomotor neurones in the rabbit. Pflügers Archiv. 1982;395:244–50. 71. 71.Perez MG, Jordan D. Effect of stimulating non-myelinated vagal axons on atrio-ventricular conduction and left ventricular function in anaesthetized rabbits. Autonomic Neuroscience. 2001;86:183–91. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1566-0702(00)00252-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11270096&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) 72. 72.Coote JH. Myths and realities of the cardiac vagus. The Journal of physiology. 2013;591:4073–85. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1113/jphysiol.2013.257758&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23878363&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) 73. 73.Brack KE, Coote JH, Ng GA. Interaction between direct sympathetic and vagus nerve stimulation on heart rate in the isolated rabbit heart. Experimental physiology. 2004;89:128–39. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1113/expphysiol.2003.002654&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15109218&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000188082300014&link_type=ISI) 74. 74.Dale LP. Maltreatment History and Physiological Regulation in Response to Physical and Emotional Stressors. 25th Annual International Trauma Conference; 2018 May; Boston, MA. 75. 75.Byrne EA, Fleg JL, Vaitkevicius PV, Wright J, Porges SW. Role of aerobic capacity and body mass index in the age-associated decline in heart rate variability. Journal of Applied Physiology. 1996;81:743–50. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8872642&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F28%2F2022.04.27.22274391.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996VD22800030&link_type=ISI)