Abstract
BACKGROUND:
Increasing evidence suggests autonomic nervous system (ANS) dysfunction may occur following mild traumatic brain injury (mTBI). Measures of heart rate, heart rate variability, blood pressure and baroreceptor sensitivity can be used to evaluate ANS dysfunction following mTBI.
OBJECTIVE:
Summarize the evidence for ANS dysfunction in adults following mTBI.
METHODS:
A search of Embase, MEDLINE, Cochrane Central Register, PsycINFO, CINAHL and SPORTDiscus databases was conducted. Search topics included: mTBI and ANS. Identified abstracts were independently reviewed by 2 reviewers followed by full text screening. Risk of bias was assessed using a modified SIGN checklist. A structured synthesis was performed.
RESULTS:
Thirty-nine studies (combined 1,467 participants diagnosed with mTBI) evaluating ANS function were included. ANS function was evaluated under various conditions including: rest, during exertion, cold pressor test, Valsalva maneuver, using face cooling and eyeball pressure paradigms. Short-term or ultra-short-term recordings were most common. The majority of studies (28/39) were rated as “unacceptable” for quality of evidence.
CONCLUSIONS:
Altered parameters of ANS function have been reported in multiple conditions following mTBI, both acutely and in the post-acute/chronic stages of recovery. However, due to methodological limitations, conclusions regarding the severity and timing of ANS dysfunction following mTBI cannot be drawn.
Keywords
Introduction
Mild traumatic brain injury (mTBI), often referred to as concussion, can occur following a biomechanical force to the head or body resultant in a neurometabolic cascade and energy crisis (Giza & Hodva, 2014). While diagnostic criteria differ based on a variety of definitions of mTBI and concussion, it is an injury often characterized by a loss of consciousness less than 30 minutes, post-traumatic amnesia less than 24 hours and Glasgow Coma Scale score of 13–15 (Mild Traumatic Brain Injury Committee of the Head Injury Interdisciplinary Special Interest Group, 1993). There is an annual estimated incidence of mTBI of 1,153/100,000 in Canada (Langer et al., 2020). Adults generally recover in the 10–14 day period following injury (McCrory et al., 2017), however up to 30% will go on to develop persistent post-concussive symptoms, including: headache (Lucas et al., 2013), dizziness (Theadom et al., 2016), fatigue (Stulemeijer et al., 2006) and exercise intolerance (Kozlowski et al., 2013). Currently, the pathophysiology underpinning mTBI recovery and symptom resolution is not fully understood; however, it is hypothesized several of the aforementioned post-concussive symptoms may be associated with autonomic dysfunction. Therefore, evaluation of autonomic function post-mTBI may provide insight into physiologic mechanisms underpinning mTBI recovery and symptom presentation.
The autonomic nervous system (ANS), composed of the parasympathetic and sympathetic branches, functions to maintain bodily homeostasis (Wehrwein et al., 2016). Sympathetic preganglionic neurons exit the ventral horns of the spinal cord, synapsing on sympathetic postganglionic neurons in the paravertebral sympathetic chain before terminating on effector organs, including the heart and lungs (Wehrwein et al., 2016). Activation of the sympathetic nervous system relaxes tracheal and bronchial smooth muscle and increases atrial muscle contractility, conduction velocity of the atrial ventricular node and heart rate (HR) (through innervation of the sinoatrial node) (Wehrwein et al., 2016). The parasympathetic branch of the ANS includes cranial nerves III, VII, IX and X (Wehrwein et al., 2016). Cranial nerve X, the vagus nerve, regulates several organ systems, including the heart and lungs. Increased vagal activation of the sinoatrial node further decreases its intrinsic firing rate, thus decreasing HR (Wehrwein et al., 2016).
Several measures can be used to evaluate ANS function. The ones most commonly employed are: heart rate variability (HRV) and cardiac baroreceptor sensitivity (BRS). HRV is normally monitored either at rest (i.e., supine, seated or upright quiet-stance) or during physiologic challenge (i.e., sit-stand protocol, isometric hand grip). HRV describes the variation in time intervals between consecutive heart beats and provides a measure of the ability of the ANS to adapt to external stimuli/challenge (Shaffer & Ginsberg, 2017). Higher time-domain HRV is generally considered to be an indicator of cardiovascular health, especially in terms of cardiovascular-related mortality (Task Force of the European Society of Cardiology and The North American Society of Pacing and Electrophysiology, 1996). Due to the exponential decay-like relationship between time-domain HRV parameters and HR (Monfredi et al., 2014), an exertional task is expected to produce an increase in HR and corresponding decrease in HRV. However, altered HR/HRV responses to physiologic challenge have been observed in individuals with mTBI compared to healthy control groups (Pertab et al., 2018). Measures of BRS provide another method of evaluating autonomic function. The baroreflex serves to regulate HR in response to blood pressure (BP) changes sensed by baroreceptors in the carotid arteries and aortic arch (Duschek et al., 2013). Research suggests BRS may be blunted post-mTBI (Pertab et al., 2018). Whether altered autonomic function post-mTBI is due to deconditioning (in those with persistent symptoms), altered central integration or another mechanism has not yet been well characterized.
Given the potential for metrics of autonomic function to both help understand the pathophysiology of post-concussive symptoms and provide an objective measure of recovery, a review of the literature including adults with all mechanisms of mTBI is warranted. While there are previous reviews of the literature in this field (Blake et al., 2016; Callaway & Kosofsky, 2019; Pertab et al., 2018; Purkayastha, 2019) they have several limitations. Namely, the prior reviews either were not systematic in nature or they did not sufficiently address/discuss the methodological limitations of the included studies. Further, as this is a rapidly advancing field of research, this review includes recently published studies that were not available to be included in prior reviews. Therefore, the objective of this review was to systematically characterize measures of autonomic function, including HRV, BP and BRS, in patients with mTBI compared to controls or in mTBI groups over time.
Materials and methods
Protocol and registration
The review was first registered in November 2019 (updated April 2020 and January 2021) with the National Institute for Health Research’s International Prospective Register of Systematic Reviews. The registered protocol (PROSPERO CRD42020157698) can be accessed at: https://www.crd.york.ac.uk/prospero/ display_record.php?RecordID = 157698. Institutional ethics approval was not required as no primary data collection occurred. The methods of this review followed Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines (Liberati et al., 2009).
Information sources
Embase (1974 –January 21, 2021), Ovid MEDLINE (1946 –January 22, 2021), Cochrane Central Register of Clinical Trials (December 2020), PsycINFO (1803 –January Week 3, 2021), CINAHL (1973 –January 2021) and SPORTDiscus (1930 –January 2021) were searched. In databases that allowed, results were restricted to human studies. Search results were limited to English language studies. Filters were applied to select for adult and young adult populations as appropriate within the specific database searches. The final search was completed on January 24, 2021. No study authors were contacted to identify investigations beyond published work. Authors were encouraged to share additional references meeting inclusion criteria.
Search
Development of the search strategy was done with assistance of a University of Calgary health sciences librarian. Search terms were related to 1) mTBI and 2) ANS. An example search strategy is provided for Embase in Table 1 and search strategies for additional databases are also available (Supplementary Tables 1-5). Records were exported directly from databases and imported into EndnoteX9 (EndNote Team, 2013). All results were exported from EndnoteX9 (EndNote Team, 2013) and imported into Covidence systematic review software where duplicates were removed.
Example search strategy
Example search strategy
Inclusion and exclusion criteria for study design, participants, outcome measures and comparators are provided in Table 2.
Eligibility criteria
Eligibility criteria
Following duplicate removal, all abstracts were screened by two independent reviewers based on eligibility criteria. Conflicts regarding inclusion/exclusion of studies were resolved by a third independent reviewer. Full text screening of included abstracts was completed by two independent reviewers and the same process was applied for conflict resolution.
Data extraction
Data extraction and risk of bias (RoB) were performed by two reviewers independently for all studies meeting the inclusion criteria following full text screening. The following data was collected: participant characteristics (age, sex/gender), mTBI diagnosis and classification, control group, assessment timepoint(s), time since injury, physiologic measure with units and duration, test condition(s) and summarized results.
Each study was assessed for RoB by two independent reviewers using modified checklists evaluated and adapted to meet Scottish Intercollegiate Guidelines Network (SIGN) criteria (Healthcare Improvement Scotland, n.d.). A modified version of the SIGN Cohort checklist was used for cross-sectional and cohort study designs. A modified version of the SIGN Case-control checklist was used for case-control study design. The injury (mTBI) was considered as the “exposure” or “case” (for case-control studies) and metrics of ANS dysfunction as the “outcome”. The SIGN checklists contain 13–15 questions with four response categories: “yes”, “no”, “not applicable” and “can’t say”. In the cohort study checklist: A) Item 1.5 regarding drop outs was used to assess whether or not drop outs were reported; B) Item 1.11 was divided into two items, 1.11 and 1.12, which respectively assessed, “the outcome measure(s) is/are reliable” and “the outcome measure(s) is/are valid”; C) For item 1.13, “exposure level of prognostic factor is assessed more than once”, symptomatology/presence of persistent symptoms was considered the “exposure” and for studies with a single timepoint symptomatology had to be assessed once for the item to be rated as “yes”. In the case-control study checklist: the items, “the outcome measure(s) is/are reliable” (1.10) and “the outcome measure(s) is/are valid” (1.11) were added. Evaluation is based on reliability and validity of subject selection, assessment, confounding factors and statistical analysis. The overall assessment of the study was based on the ability to minimize RoB. Discrepancies were discussed and if a consensus was not achieved by two reviewers, conflicts were resolved by a third independent reviewer.
The corresponding authors for three included studies were contacted. The first researcher was contacted to provide supplementary tables referenced in the manuscript (Ding et al., 2020) (replied), the second researcher was contacted to clarify results (Huang et al., 2019) (did not reply) and a third author was contacted regarding result discrepancies within the manuscript (Kozlowski et al., 2013) (replied).
Data synthesis
Due to heterogeneity of study methods, participant demographics, outcome measures and mTBI diagnostic criteria, data was not pooled for a meta-analysis. Demographics and results were summarized.
Data availability
Data sharing is not applicable to this article as no new data were created or analyzed in this study.
Results
Study selection
A total of 4,311 studies were identified, of which, 1,064 duplicates were removed. Two references were shared by a study author (Haider et al., 2020; Pyndiura et al., 2020). The title and abstract of 3,247 studies were screened with 3,092 not meeting inclusion criteria (see Table 2). Full texts of 155 studies were assessed for eligibility and 116 were excluded. The most common reason for exclusion was incorrect study design. Ultimately, 39 studies were included in the qualitative synthesis (see Fig. 1 for PRISMA flowchart).

PRISMA flowchart.
Study characteristics, participant demographics and RoB grade are summarized in Table 3.
Study characteristics
Study characteristics
Abbreviations: ApEN, approximate entropy; BDI, Beck Depression Inventory; BMI, body mass index; BP, blood pressure; BRS, baroreceptor sensitivity; CISG, Concussion in Sport Group; COMPASS-31 = Composite Autonomic Symptom Score 31 questionnaire; CPT, cold pressor test; CT, computerized tomography; DBP, diastolic blood pressure; EP, eyeball pressure stimulation; FC, face cooling; HF = high frequency; HR, heart rate; HRV, heart rate variability; HRSD, standard deviation of heart rate; hx = history; hx-, no mTBI history; hx3 +, history of≥3 mTBI; ICHD, International Classification of Headache Disorders; IHGC, isometric hand grip contraction; LF, low frequency; LOC, loss of consciousness; MAP, mean arterial pressure; MD, medical doctor; mTBI, mild traumatic brain injury; pNN50, percentage of successive RR intervals that differ by more than 50ms; PTA, post-traumatic amnesia; PTSD, Post-traumatic stress disorder; QVTI, QT interval variability index; RMSSD, root mean square of successive RR interval differences; RoB, risk of bias; SampEN, sample entropy; SBP, systolic blood pressure; SDNN, standard deviation of NN intervals; SysSlope, systolic slope; SRT, Sharpened Romberg Test (balance test); TBI, traumatic brain injury; TP, total power; VM, valsalva maneuver; 5SS, 5s squat, 5s stand (squat-stands); 10SS, 10s squat, 10s stand (squat-stands). †using progressive treadmill test. ‡exercise testing only completed if symptoms had returned to baseline; if participants were symptomatic follow-ups were conducted at 3d intervals until athlete was asymptomatic at rest/had returned to play. §data from same cohort. ¶systolic slope (SysSlope) of the arterial pressure wave = calculated from the rate of rise (change in pressure divided by change in time) of the systolic upstroke. ||abnormal SRT = loss of balance in first 30s of 60s balance test (across 4 trials) where participant stands with one foot in front of the other, toes touching heel, arms crossed across chest, eyes closed. #uncomplicated mTBI (meeting diagnostic criteria for mTBI with no intercranial lesions on conventional imaging). ††normative data was taken from the average of each 5min segment over 24hrs recorded from 240 controls; mTBI did not follow the same experimental procedure. ‡‡data from same cohort. aSung et al. Psychophysiol. bSung et al. Clin Neuropsychol.
Of the included studies, twenty-one were cross-sectional studies (Abaji et al., 2016; Bishop et al., 2017; Ding et al., 2020; Dobney et al., 2018; Hanna-Pladdy et al. 2001; Hilz et al., 2020; Howard et al., 2018; Huang et al., 2019; Johnson et al., 2020; Kozlowski et al., 2013; Mirow et al., 2016; Pyndiura et al., 2020; Russell et al., 2020; Solbakk et al., 2005; Tan et al., 2009; Truong & Ciuffredia, 2016), sixteen were cohort studies (Brandt et al., 2020; Clausen et al., 2016; Dobson et al., 2017; Gall et al., 2004; La Fountaine et al., 2018; Leddy et al., 2010; Liao et al., 2016; Purkayastha, 2019; Senthinathan et al., 2017; Sung, 2016; Sung, 2016; Wright et al. 2017, 2018), while two (Haider et al., 2020; Hutchison et al., 2017) were case-control studies. All studies, with the exception of two (Leddy et al., 2010; Mirow et al., 2016), compared mTBI participants to non-head injured control groups. Of these, five studies (Gall et al., 2004; Huang et al., 2019; Hutchison et al., 2017; Purkayastha, 2019; Senthinathan et al., 2017) included controls with a remote history of mTBI, three studies (Johnson et al., 2020; Russell et al., 2020) included controls with no mTBI in the last twelve months and eleven studies (Bishop et al., 2017; Clausen et al., 2016; Hilz et al., 2020; La Fountaine et al., 2018; Solbakk et al., 2005; Tan et al., 2009) did not report on mTBI history of the control group. Two studies (Wright et al., 2018) compared post-mTBI values to preseason values.
In twenty-two studies (Abaji et al., 2016; Bishop et al., 2017; Ding et al., 2020; Dobney et al., 2018; Haider et al., 2020; Hanna-Pladdy et al., 2001; Hilz et al., 2020; Howard et al., 2018; Huang et al., 2019; Johnson et al., 2020; Kozlowski et al., 2013; Mirow et al., 2016; Pyndiura et al., 2020; Russell et al., 2020; Solbakk et al., 2005; Tan et al., 2009; Truong & Ciuffreda, 2016) autonomic dysfunction was assessed at a single timepoint (within 72 hours to greater than 10 years post-injury). Multiple assessments were conducted in fifteen studies (Brandt et al., 2020; Dobson et al., 2017; Gall et al., 2004; Hutchison et al., 2017; La Fountaine et al., 2018; Liao et al., 2016; Purkayastha, 2019; Senthinathan et al., 2017; Sung, 2016; Sung, 2016; Wright et al., 2018). Timing of follow up assessments ranged from within 24 hours up to 18 months after the initial assessment. Outcomes were assessed pre and post sub-symptom threshold exercise interventions in two studies (Clausen et al., 2016; Leddy et al., 2010).
Demographics
A total of 1,467 mTBI participants were included in the thirty-nine studies. The number of mTBI participants across studies ranged from 3 –331 with a mean of 38±60. Eight studies included only male participants (Abaji et al., 2016; Bishop et al., 2017; Gall et al., 2004; La Fountaine et al., 2018; Mirow et al., 2016; Tan et al., 2009; Wright et al., 2017; Wright et al., 2018) and in most studies reporting sex/gender, males outnumbered female mTBI participants with a male:female ratio of 768:619. Two studies (Dobney et al., 2018; Truong & Ciuffreda, 2016) did not specify the sex/gender of mTBI participants; one study (Truong & Ciuffreda, 2016) did not include sex/gender of mTBI or control participants. Twenty studies (Brandt et al., 2020; Dobney et al., 2018; Dobson et al., 2017; Haider et al., 2020; Hanna-Pladdy et al., 2001; Hilz et al., 2017; Howard et al., 2018; Huang et al., 2019; Hutchison et al., 2017; Johnson et al., 2020; Kozlowski et al., 2013; La Fountaine et al., 2011; Pyndiura et al. 2020; Russell et al., 2020; Senthinathan et al., 2017; Solbakk et al., 2005) reported participant sex and eight studies (Hilz et al., 2020; La Fountaine et al., 2016; Leddy et al., 2010; Liao et al., 2016; Mirow et al., 2016; Sung, 2016; Tan et al., 2009) reported participant gender. Ten studies (Abaji et al., 2016; Bishop et al., 2017; Clausen et al., 2016; Ding et al., 2020; Gall et al., 2004; La Fountaine et al., 2018; Purkayastha, 2019; Sung, 2016; Wright et al., 2018) did not specify if reporting sex/gender or interchanged the two terms.
The mean age of mTBI participants across studies ranged from 18.1 to 41.2 years with an overall mean age of 26.9±8.4 years. Hutchison et al. (Hutchison et al., 2017) did not separately report the age of mTBI and control groups, stating that the average age was not significantly different.
Injury characteristics
Time since injury varied greatly between studies, with nine studies assessing participants within 72 hours of mTBI (Bishop et al., 2017; Dobson et al., 2017; Gall et al., 2004; La Fountaine et al. 2009, 2018; Russell et al., 2020; Wright et al., 2018), nine studies within the first 7 days (Huang et al., 2019; Hutchison et al., 2017; Johnson et al., 2020; La Fountaine et al., 2011; Liao et al., 2016; Purkayastha, 2019; Senthinathan et al., 2017; Sung, 2016; Sung, 2016), one within 10 days (Johnson et al., 2018), one within the first 2 weeks (Brandt et al., 2020) and the remaining seventeen studies (Abaji et al., 2016; Clausen et al., 2016; Ding et al., 2020; Haider et al., 2020; Hanna-Pladdy et al., 2001; Hilz et al., 2020; Howard et al., 2018; Kozlowski et al., 2013; Leddy et al., 2010; Mirow et al., 2016; Pyndiura et al., 2020; Solbakk et al., 2005; Truong & Ciuffreda, 2016) ranged from > 45 days to > 10 years post-mTBI. Time since injury was not reported in two studies (Dobney et al., 2018; Tan et al., 2009).
Seventeen studies (Abaji et al., 2016; Bishop et al., 2017; Dobson et al., 2017; Gall et al., 2004; Haider et al., 2020; Huang et al., 2019; Hutchison et al., 2017; Johnson et al., 2020; La Fountaine et al., 2018; Purkayastha, 2019; Senthinathan et al., 2017; Wright et al., 2018) included participants with sport-related mTBI, thirteen (Brandt et al., 2020; Ding et al., 2020; Hanna-Pladdy et al., 2001; Hilz et al., 2020; Howard et al., 2018; Kozlowski et al., 2013; Leddy et al., 2010; Liao et al., 2016; Solbakk et al., 2005; Sung, 2016; Sung, 2016; Truong & Ciuffreda, 2016) included multiple mechanisms of injury for mTBI (including motor vehicle collision, fall, sport-related, assault or other) and three studies with participants who sustained an mTBI in military combat (Mirow et al., 2016; Russell et al., 2020; Tan et al., 2009). Six studies (Clausen et al., 2016; Dobney et al., 2018; Hilz et al., 2017; Pyndiura et al., 2020) did not specify the mechanism of injury.
Fourteen studies (Bishop et al., 2017; Brandt et al., 2020; Clausen et al., 2016; Ding et al., 2020; Gall et al., 2004; Howard et al., 2018; Hung et al., 2019; Johnson et al., 2020; Kozlowski et al., 2013; La Fountaine et al., 2009; Mirow et al., 2016; Solbakk et al., 2005; Truong & Ciuffreda, 2016) included mTBI participants who were symptomatic. Participants (mTBI) were asymptomatic in six studies (Abaji et al., 2016; Dobney et al., 2018; Haider et al., 2020; Hilz et al., 2016; Pyndiura et al., 2020). One study (Hanna-Pladdy et al., 2001) compared symptomatic and asymptomatic mTBI participants and seven studies followed participants across timepoints (acutely symptomatic until they were asymptomatic) (Leddy et al., 2010; Hutchison et al., 2017; La Fountaine et al., 2011; Purkayastha, 2019; Senthinathan et al., 2017; Wright et al., 2018). Eleven studies did not report on symptoms at the time of assessment(s) (Dobson et al., 2017; Hilz et al., 2020; La Fountaine et al., 2018; Liao et al., 2016; Russell et al., 2020; Sung, 2016; Sung, 2016; Tan et al., 2009).
Assessment of autonomic dysfunction
Autonomic dysfunction was assessed at rest in thirty-four of the thirty-nine studies (Abaji et al., 2016; Bishop et al., 2017; Brandt et al., 2020; Ding et al., 2020; Dobney et al., 2018; Dobson et al., 2017; Gall et al., 2004; Haider et al., 2020; Hanna-Pladdy et al., 2001; Hilz et al., 2020; Huang et al., 2019; Hutchison et al., 2017; Johnson et al., 2020; Kozlowski et al., 2013; La Fountaine et al., 2018; Liao et al., 2016; Mirow et al., 2016; Purkayastha, 2019; Pyndiura et al., 2020; Russell et al., 2020; Senthinathan et al., 2017; Sung, 2016; Sung, 2016; Tan et al., 2009; Wright et al., 2018). Nine studies assessed autonomic dysfunction during a postural change (sit-to-stand, squat-stands) (Bishop et al., 2017; Ding et al., 2020; Dobson et al., 2017; Hilz et al., 2017; Pyndiura et al., 2020; Russell et al., 2020; Senthinathan et al., 2017; Wright et al., 2018), eight studies during exertion (treadmill, cycling or isometric hand grip contraction) (Abaji et al., 2016; Clausen et al., 2016; Dobney et al., 2018; Gall et al., 2004; Kozlowski et al., 2013; La Fountaine et al., 2016, Leddy et al., 2010), two studies during Valsalva maneuver (Dobson et al., 2017; Hilz et al., 2016), two studies during face cooling (Haider et al., 2020; Johnson et al., 2018) and two studies during a cognitive task (Hanna-Pladdy et al., 2001; Huang et al., 2019). Other conditions included eyeball pressure (Hilz et al., 2015), olfactory (Hilz et al., 2020) or visual (Solbakk et al., 2005) stimuli, cold pressor test (Johnson et al., 2020) and various light conditions (for assessment of pupillary light reflex) (Truong & Ciuffreda, 2016).
The most commonly assessed measures of autonomic function were: HR (assessed in 87% of studies (Abaji et al., 2016; Bishop et al., 2017; Clausen et al., 2016; Ding et al., 2020; Dobney et al., 2018; Dobson et al., 2017; Gall et al., 2004; Haider et al., 2020; Hanna-Pladdy et al., 2001; Hilz et al., 2020; Huang et al., 2019; Hutchison et al., 2017; Johnson et al., 2020; Kozlowski et al., 2013; La Fountaine et al., 2018; Leddy et al., 2010; Liao et al., 2016; Purkayastha, 2019; Pyndiura et al., 2020; Russell et al., 2020; Senthinathan et al., 2017; Solbakk et al., 2005; Sung, 2016; Sung, 2016; Wright et al., 2018), frequency-domain HRV metrics (62% of studies (Abaji et al., 2016; Bishop et al., 2017; Brandt et al., 2020; Ding et al., 2020; Dobney et al., 2018; Gall et al., 2004; Haider et al., 2020; Hilz et al., 2017; Huang et al., 2019; Hutchison et al., 2017; Johnson et al., 2018; La Fountaine et al., 2018; Liao et al., 2016; Mirow et al., 2016; Purkayastha, 2019; Pyndiura et al., 2020; Russell et al., 2020; Senthinathan et al., 2017; Sung, 2016; Sun, Lee et al., 2016)), BP (49% of studies (Bishop et al., 2017; Clausen et al., 2016; Ding et al., 2020; Dobney et al., 2018; Dobson et al., 2017; Hilz et al., 2020; Johnson et al., 2020; Kozlowski et al., 2013; La Fountaine et al., 2018; Leddy et al., 2010; Purkayastha, 2019; Wright et al., 2018)), time-domain HRV metrics (46% of studies (Abaji et al., 2016; Bishop et al., 2017; Dobney et al., 2018; Gall et al., 2004; Haider et al., 2020; Hilz et al., 2017; Huang et al., 2019; Hutchison et al., 2017; Johnson et al., 2018; La Fountaine et al., 2011; Liao et al., 2016; Purkayastha, 2019; Pyndiura et al., 2020; Sung, 2016; Sung, 2016; Tan et al., 2009)), BRS (15% of studies (Ding et al., 2020; Dobson et al., 2017; Hilz et al., 2017; La Fountaine et al., 2018)), non-linear HRV metrics (10% of studies (Abaji et al., 2016; Bishop et al., 2017; La Fountaine et al., 2009; Senthinathan et al., 2017)) and frequency-domain BP variability (10% of studies (Ding et al., 2020; Hilz et al., 2017)). Other metrics included measures of pupillary reflex (Truong & Ciuffreda, 2016) and BP and/or HR response to Valsalva maneuver (Dobson et al., 2017; Hilz et al., 2016). The Composite Autonomic Symptom Score 31 (COMPASS-31) questionnaire was evaluated in one study (Howard et al., 2018).
Risk of Bias (RoB)
Results from the RoB assessment are reported in Table 4. No studies were graded as “high-quality”. Eleven studies (Ding et al., 2020; Howard et al., 2018; Hutchison et al., 2017; Johnson et al., 2020; Kozlowski et al., 2013; Leddy et al., 2010; Purkayastha, 2019; Solbakk et al., 2005; Truong & Ciuffreda, 2016; Wright et al., 2018) were of “acceptable” quality, with the remaining 28 (Abaji et al., 2016; Bishop et al., 2017; Brandt et al., 2020; Clausen et al., 2016; Dobney et al., 2018; Dobson et al., 2017; Gall et al., 2004; Haider et al., 2020; Hanna-Pladdy et al., 2001; Hilz et al., 2020; Huang et al., 2019; Johnson et al., 2018; La Fountaine et al., 2018; Liao et al., 2016; Mirow et al., 2016; Pyndiura et al., 2020; Russell et al., 2020; Senthinathan et al., 2017; Sung, 2016; Sung, 2016; Tan et al., 2009) graded as “unacceptable”. The most common biases were a lack of reporting how many people were approached to participate, not specifying whether the assessment of outcomes was blinded to exposure status, unreliable measures of exposure (not specifying mTBI diagnostic criteria), unreliable or invalid outcome measure and not addressing potential confounders.
Risk of bias of included studies
Risk of bias of included studies
Reporting of study confounders
A summary of HR and HRV results is presented in Supplementary Table 6. BP, BRS and alternative autonomic measure results are presented in Supplementary Table 7. Correlations between metrics of autonomic function and injury characteristics are reported in Supplementary Table 7.
Discussion
In recent years, measurement of physiologic biomarkers of recovery has become a growing field of mTBI research. This includes assessment of autonomic parameters, such as those included in this review. These measures have the potential to provide an objective assessment of recovery and are therefore of great interest and clinically relevant. Evaluation of measures of autonomic function may prove especially important in the context of return-to-sport and allow for evaluation of physiologic changes that may persist beyond symptom resolution. The presented results offer evidence that a degree of ANS dysfunction may occur as a sequela of mTBI in adult cohorts. However, interpretation of the data is complicated by the wide variety of methods, reported variables, timepoints and patient populations evaluated in the literature. Additionally, significant covariates that may potentially confound or modify the relationship between mTBI and measures of autonomic function were not evaluated in many cases. While there is a paucity of literature with low RoB, the key findings of this review can be summarized as follows: It remains unclear if there are differences in HR, HRV and BP parameters at rest (in a range of supine, seated and upright quiet-stance postures) in individuals with mTBI compared to controls, although most studies suggest there is not a significant difference at rest. A wide range of (non-rest) test conditions have revealed differences in measures of autonomic function between mTBI and control groups. These include exertional tasks (i.e., isometric hand grip, squat-stands, submaximal treadmill test), Valsalva maneuver, face cooling, eyeball pressure, postural change (i.e., supine-to-stand, sit-to-stand), odorous and visual stimulations, cold pressor test, and measures of pupillary light reflex, among others. Findings suggest employing a physiologic challenge may be needed to reveal subtle differences in autonomic function between groups with and without a history of mTBI. A wide range of timepoints have been assessed in the literature from within 48 hours of injury (Dobson et al., 2017; La Fountaine et al., 2018) up to 15 years post-injury (Hanna-Pladdy et al., 2001). The majority of studies included were cross-sectional (one timepoint) (Abaji et al., 2016; Bishop et al., 2017; Ding et al., 2020; Dobney et al., 2018; Haider et al., 2020; Hanna-Pladdy et al., 2001; Hilz et al., 2020; Howard et al., 2018; Huang et al., 2019; Johnson et al., 2020; Kozlowski et al., 2013; Mirow et al., 2016; Pyndiura et al., 2020; Russell et al., 2020; Solbakk et al., 2005; Tan et al., 2009; Truong & Ciuffreda, 2016), whereas others assessed temporal changes (Brandt et al., 2020; Clausen et al., 2016; Dobson et al., 2017; Gall et al., 2004; Hutchison et al., 2017; La Fountaine et al., 2018; Leddy et al., 2010; Liao et al., 2016; Purkayastha, 2019; Senthinathan et al., 2017; Sung, 2016; Sung, 2016; Wright et al., 2018). While several studies found that between group differences noted acutely post-injury resolved at repeat assessments (Dobson et al., 2017; La Fountaine et al., 2011; Purkayastha,Williams et al., 2019), other studies evaluating one timepoint in the post-acute/chronic stage reported between group differences. The relationship between ANS dysfunction and post-concussive symptoms requires further investigation. None of the papers included in this review had a primary aim of characterizing this relationship. Especially in studies completing a single assessment in the post-acute/chronic stage of injury, clear symptom reporting is necessary to understand generalizability. Whether the participants (or a subset of the participants) are experiencing persistent symptoms should be described.
Proposed pathophysiology of ANS dysfunction following mTBI
The underlying pathophysiology of ANS dysfunction following mTBI remains unknown; however, pre-clinical and clinical work suggests this may result from either: a) diffuse axonal injury during the initial mTBI; b) reduced cerebral blood flow (CBF) post-mTBI or c) deconditioning following the injury. Higher brain centers, including the prefrontal cortex (anterior cingulate, insula, orbitofrontal and ventromedial cortex) and limbic structures (amygdala and hypothalamus) in addition to structures in the brainstem (periaqueductal gray matter, nucleus ambiguus, medulla) form the Central Autonomic Network (Benarroch, 1993). Therefore, depending on the mechanism of injury and whether there is diffuse axonal injury, centers responsible for control of HRV and other autonomic function may be affected.
CBF has been increasingly investigated following mTBI using a variety of methods including transcranial Doppler and arterial spin labeling (ASL) magnetic resonance imaging (MRI). Using ASL MRI, Meier et al. reported reduced CBF in the right dorsal midinsular cortex in collegiate athletes in the first week post-mTBI compared to controls (Meier et al., 2015). This difference did not persist at 1 month post-injury. Athletes who were slow to recover, characterized as > 14 days until clearance for return-to-play, had lower CBF in the right dorsal midinsular cortex at 1 month compared to those with a more rapid return-to-play (Meier et al., 2015). As noted previously, the insula is part of the Central Autonomic Network, playing an important role in cardiovascular ANS integration (Bernaroch, 1993). While preliminary studies suggest a link between ANS dysfunction and altered CBF, additional research is required.
Following mTBI, individuals may struggle with exercise intolerance (Kozlowski et al., 2013), characterized as the inability to return to one’s premorbid level of physical activity due to symptom exacerbation. This may result in reduced physical activity and subsequent deconditioning. It is understood that deconditioning (reduced exercise capacity) is associated with decreased parasympathetic (vagal) activity and thus an increased resting HR (Gourine & Ackland, 2019). For example, improvement in HR recovery, a measure of parasympathetic activity, observed following a 6 week training program is lost after only 4 weeks of detraining (Heffernan et al., 2007). Individuals with higher resting HR and lower percentage predicted peak oxygen-consumption have been shown to have slower HR recovery following exercise, a measure of an individual’s ability to recruit cardiac vagal tone (Machhada et al., 2017). Unfortunately, the link between physical activity/exercise intolerance and ANS function in individuals with mTBI has not been well characterized.
Methodological confounders
Previous reviews of the ANS literature (Blake et al., 2016; Callaway & Kosofsky, 2019; Pertab et al., 2018; Purkayastha, 2019) have failed to highlight the significant methodological confounders present in the included studies. Additionally, several of these were not systematic reviews, nor did they include participants will all mechanisms of injury. The current systematic review addressed these gaps in the literature by including a wide range of autonomic function outcomes, completing an RoB assessment, providing a comprehensive review of confounders for each study and outlining future directions for ANS research in the mTBI context. Advances in this field will be dependent on the methodological rigor of future research. Additional studies addressing confounders and with larger sample sizes will be necessary before these measures may be implemented clinically.
Providing a clear operational definition of the “exposure” (mTBI) is important for reproducibility and generalizability of the research. In the case of eighteen studies reliability and/or validity of exposure status could not be ascertained based on presented information. Future research should use accepted mTBI diagnostic criteria, such as the American Congress of Rehabilitation Medicine criteria (Mild Traumatic Brain Injury Committee of the Head Injury Interdisciplinary Special Interest Group, 1993) or the 5th Consensus statement on concussion in sport criteria (McCrory et al., 2017). Injury was self-reported with no additional detail or diagnostic criteria in three studies (Dobney et al., 2018; Hanna-Pladdy et al., 2001; Pyndiura et al., 2020). It is also important to consider the mTBI history of control groups. In eighteen studies (Bishop et al., 2017; Clausen et al., 2016; Gall et al., 2004; Hilz et al., 2020; Huang et al., 2019; Hutchison et al., 2017; Johnson et al., 2020; La Fountaine et al., 2018; Purkayastha, 2019; Russell et al., 2020; Senthinathan et al., 2017; Tan et al., 2009), controls either did have a reported lifetime history of mTBI or insufficient information on mTBI history was presented. Particularly when evaluating ANS function in the chronic stage of injury or in those with persistent symptoms, control groups should be free of lifetime mTBI history to confidently report between group differences.
Many studies failed to report confounders that could affect ANS function data, such as other medical conditions, medications, cardiorespiratory fitness and stage of menstrual cycle. Despite many of the cohorts being youth athletes, screening and recording of past medical history and medications affecting ANS function should be clearly reported. Unfortunately, this was not the case in many studies and limits the internal validity. This is especially important in older cohorts and those with persistent symptoms where polypharmacy may be more likely. Cardiorespiratory fitness should also be considered, especially in athletic populations. As mentioned previously, resting HR, and thus HRV (Monfredi et al., 2014), is affected by cardiorespiratory fitness (Machhada et al., 2017). Several studies (Clausen et al., 2016; Haider et al., 2020; Johnson et al., 2018) compared varsity athletes to recreational control athletes. In the absence of either a direct or indirect index of cardiorespiratory fitness completed by both groups, interpretation of study findings is limited. Those with higher cardiorespiratory fitness, and thus lower resting HR, would have higher time-domain HRV measures, compared to their peers with lower cardiorespiratory fitness, thus producing between group differences unrelated to mTBI history. Stage of the menstrual cycle is also known to affect HRV. A recent meta-analysis highlighted fluctuation in cardiac vagal activity across the menstrual cycle with a significant decrease from the follicular to luteal phase (Schmalenberger et al., 2019). Analysis comparing cycle sub-phases (i.e., menstrual, mid-to-late follicular, ovulatory, early-to-mid luteal, premenstrual) revealed a decrease in measures of cardiac vagal activity from menstrual to premenstrual sub-phases and mid-to-late follicular to premenstrual sub-phases, indicating that the change from follicular to luteal phases is likely driven by a decline in the premenstrual sub-phase (of the luteal phase) (Schmalenberger et al., 2019). The meta-analysis pooled multiple autonomic measures (HF, RMSSD, BRS) and separate analyses for each variable were not completed (Schmalenberger et al., 2019). Unfortunately, none of the included studies in our review reported if participants were eumenorrheic or menstrual cycle stage, which should be included in future studies.
Multiple studies did not specify if reporting sex/gender or interchanged the two terms and no studies explicitly reported both sex and gender. Sex refers to a set of biological attributes (primarily associated with physical and physiological features including chromosomes, gene expression, hormone levels and function, and reproductive/sexual anatomy), whereas gender refers to the socially constructed roles, behaviours, expressions and identities of girls, women, boys, men, and gender diverse people (Canadian Institutes of Health Research, 2020). Future research should consider reporting both sex and gender. There is a lack of evaluation of both sex and gender differences in autonomic function following mTBI.
Length of recording periods are generally divided as ultra-short-term (<3 minutes), short-term (5-20 minutes) or 24 hour recordings. It is inappropriate to compare HRV metrics when collected from epochs of different lengths, as longer recordings are associated with increased HRV (Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology, 1996). Unfortunately, one study (Tan et al., 2009) did compare data from short-term recordings (mTBI participants) to 24 hour control data. Ultra-short-term recordings for analysis of HRV remain controversial, yet these were used in the majority of included studies (Abaji et al., 2016; Brandt et al., 2020; Haider et al., 2020; Hilz et al., 2017; Huang et al., 2019; Johnson et al., 2018; La Fountaine 2009, 2011; Russell et al., 2020). Recent research suggests that recordings of at least 4 minutes should be used to approximate short-term data (Burma et al., 2021). Using 1 minute recordings, metrics can vary by ≥30% when compared to a 5 minute recording (Burma et al., 2021). Based on these findings, analysis using epochs < 4 minutes were rated as unreliable/invalid in the RoB assessment.
Recording methods should also be carefully considered. While higher-end HR monitors may obtain fairly close values to an electrocardiogram (ECG) at rest and during exertion (pNN50 being an exception) (Vanderlei et al., 2008), ECG remains the gold-standard for measurement of HRV. For specific variables, such as QTVI, a 12-lead ECG should be used in place of a 3-lead ECG (which was used in an included study (La Fountaine et al., 2011)) as electrode placement can cause variation in QT time of 25-30 ms (Salvi et al., 2012). Respiratory rate is also known to affect certain HRV variables, specifically HF. Respiratory rate was measured in nine studies (Hilz et al., 2020; Huang et al., 2019; La Fountaine et al., 2018; Purkayastha, 2019; Solbakk et al., 2005). HF is commonly referred to as the respiratory frequency as it includes frequencies between 0.15–0.4 Hz, corresponding to physiologic events occurring 9–24 times per minute. At rest, respirations occur ∼12–20 times per minute, thus falling into the HF range (Singh et al., 2018); however, altered breathing faster than the normal range, potentially due to test anxiety, would elevate this rate above the HF band and thus skew interpretation of results. Therefore, future studies investigating HRV following mTBI should ensure there is a steady and consistent rate and depth of breathing for all participants at all time points. Additional considerations should be taken for certain test conditions. For example, studies using face cooling (Haider et al., 2020; Johnson et al., 2018) or cold pressor test (Johnson et al., 2020) conditions should address whether participants (in the case of the included studies, athletes) use ice baths in their training regimen as this may affect response.
Small sample size is another limitation of the included research given there is large between-subject variability for many HRV variables (Nunan et al., 2010). Of the included studies, only twelve had a sample size ≥25 in both groups (mTBI, controls) (Brandt et al., 2020; Dobney et al., 2018; Hanna-Pladdy et al., 2011; Hilz et al., 2016; Howard et al., 2018; Hutchison et al., 2017; Liao et al., 2016; Pyndiura et al., 2020l; Russell et al., 2020; Sung, 2016; Sung, 2016; Truong & Ciuffreda, 2016). The n = 25 sample size as been used as a cut-off for small sample size papers in other systematic reviews of this population (Iverson et al., 2020). Notably, only three of the studies (Howard et al., 2018; Hutchison et al., 2017; Truong & Ciuffreda, 2016) with a larger sample (≥25 in both groups) were rated as “acceptable” on the risk of bisk assessment. Ability to detect between-group differences may be affected by these small sample sizes with potentially great inter-subject variability. A lack of evaluation for confounding/modification by factors that may have affected test results in addition to the use of unreliable/invalid autonomic function outcome measures (i.e., ultra-short-term recordings) are reflected in the RoB assessment where 72% (28 out of 39) of studies (Abaji et al., 2016; Bishop et al., 2017; Brandt et al., 2020; Clausen et al., 2016; Dobney et al., 2018; Dobson et al., 2017; Gall et al., 2004; Haider et al., 2020; Hanna-Pladdy et al., 2001; Hilz et al., 2020; Huang et al., 2019; Johnson et al., 2018; La Fountaine et al., 2018; Liao et al., 2016; Mirow et al., 2016; Pyndiura et al., 2020; Russell et al., 2020; Senthinathan et al., 2017; Sung, 2016; Sung, 2016; Tan et al., 2009) were rated as “unacceptable”.
Participant symptomatology
Whether a degree of ANS dysfunction persists in some or all cases following symptom resolution remains a question of great clinical interest. The relationship between specific post-concussive symptoms and autonomic dysfunction also requires further study. With return to sport protocols currently based on symptom exacerbation with increasing exertion/physical contact (McCrory et al., 2017), assessment of specific autonomic parameters could offer an objective measure of recovery. Preliminary data suggests certain autonomic measures may be altered in asymptomatic individuals in the post-acute stage of injury (Abaji et al., 2016; Haider et al., 2020), although these studies are limited by methodological confounders (see Table 5). In general, symptomatology is poorly characterized in the included literature, with sixteen studies (Abaji et al., 2016; Dobson et al., 2017; Gall et al., 2004; Hilz et al. 2011, 2017; Johnson et al., 2020; La Fountaine et al., 2018; Pyndiura et al., 2020; Senthinathan et al., 2017; Solbakk et al., 2005; Tan et al., 2009; Truong & Ciuffreda, 2016) not reporting results of a specific symptom questionnaire nor a general symptom battery. Unfortunately, several studies noted a symptom checklist had been completed, but results were not reported in the manuscript (Gall et al., 2004; Johnson et al., 2020; La Fountaine et al., 2009; Senthinathan et al., 2017; Tan et al., 2009). Additionally, the heterogeneity of questionnaires does not allow for easy comparison between studies and standardization should be considered. There is a great need for a study with the primary aim of evaluating change in multiple post-concussive symptoms and autonomic metrics across recovery timepoints in a large cohort. Finally, study of ANS dysfunction in cohorts with persistent post-concussive symptoms is limited. While several studies did complete assessments in the chronic stage of injury (>3 months post), symptoms were poorly characterized in many cases (Hilz et al., 2017; Pyndiura et al., 2020; Solbakk et al., 2005). Research in cohorts with persistent post-concussive symptoms, diagnosed using standardized criteria, is necessary to understand potential underlying physiologic mechanisms contributing to persistent symptoms, such as exercise intolerance, dizziness and headache.
Limitations
There are several limitations based on the methods of this review. Only published work was included in the review, thus the findings are subject to publication bias. Additionally, only English language papers were included. Our results and conclusions are therefore subject to language-of-publication bias. The majority of studies included in this review had a high RoB due to methodological shortcomings. As pre-clinical studies were not included, this review does not address translational research on this topic. It is understood that different critical appraisal tools uniquely assess the methodological quality of a publication and different RoB tools have been shown to lead to different conclusions (ratings of study quality) (Losilla et al., 2018). The SIGN checklist was chosen for this review as it is a comprehensive, easy to use RoB assessment with checklists tailored to different study designs (i.e., cohort, case-control studies). The heterogeneity of methods, outcomes and samples precluded a quantitative review or meta-analysis.
Future directions
There is great opportunity for the use of ANS function parameters, such as HRV, to understand recovery following mTBI and potentially guide return to play; however, the significant limitations of the presented literature limit its clinical applicability. We propose the following recommendations to build upon the current literature and ensure methodological rigor in the field. Additional research is needed to evaluate change in autonomic parameters across time points. While several studies have assessed metrics across time, due to a lack of control for confounders, small sample sizes and lack of strong methodological designs, future research is needed to understand the relationship between mTBI and ANS changes that may occur following injury. Completing both resting and exertional tasks at all timepoints should be considered. Symptomatology of study cohorts should be characterized using commonly used symptom batteries (i.e., symptom evolution score in SCAT5, or PCSS in ImPACT). Symptom data should be reported alongside all autonomic data. This is necessary to understand the link between specific post-concussive symptoms and autonomic function. Recording duration must be considered in the interpretation of HRV results. Ultra-short-term recordings (<3 minutes) remain controversial and recent data suggest that a recording of at least 4 minutes is required to approximate short-term recordings (>5 minutes) (Burma et al., 2021). Increasing methodological rigor in the field requires addressing confounders in the study design. These include reporting on past medical history and medications of participants, consumption of caffeine/tobacco and participation in exercise prior to data collection. There is a lack of investigation of sex/gender differences across recovery trajectory and test conditions. We know that females may be more likely to develop persistent post-concussive symptoms (Iverson et al., 2017), however this has not been investigated in the context of ANS dysfunction. Multiple studies included only male participants (Abaji et al., 2016; Bishop et al., 2017; Gall et al., 2004; La Fountaine et al., 2018; Mirow et al., 2016; Tan et al., 2009; Wright et al., 2018). Furthermore, none of the included studies reported on stage of the menstrual cycle, although this is known to influence HRV (Schmalenberger et al., 2020). Additionally, for certain test conditions, such as the isometric hand grip, there are known sexually dimorphic differences in force steadiness (Jakobi et al., 2018) that should be accounted for. This is also a reflection of the relatively small sample sizes of most studies that are not adequately powered to investigate these questions.
Conclusion
Measures of autonomic function, such as HRV and BRS, are objective, quantitative measures with the potential to be used to assess recovery trajectory following mTBI. This review summarized a broad scope of methods and variables to examine changes in autonomic function following mTBI. Currently, recovery is largely assessed by subjective metrics, such as symptom questionnaires. While objective measures of recovery continue to be investigated, there is significant opportunity for the use of these novel ANS provocation methods to evaluate both recovery across timepoints and the underlying pathophysiology of persistent symptoms. However, additional research using more rigorous methodology to enable evaluation of multiple covariates will be required before these tests can be used clinically. The potential for evaluation of these metrics across timepoints to inform return-to-sport, in addition to traditional symptom ratings, is of particular interest in the field.
Footnotes
Acknowledgments
The authors would like to acknowledge the help of Lorraine Toews, a health science librarian, for providing her expertise in the development of the search strategy for this review.
Conflict of interest
The authors have no conflict of interest to disclose.
The supplementary tables are available from
.
