Article Text


Heart rate variability following youth concussion: how do autonomic regulation and concussion symptoms differ over time postinjury?
  1. Melissa Paniccia1,
  2. Lee Verweel1,
  3. Scott G Thomas2,
  4. Tim Taha2,
  5. Michelle Keightley3,4,5,
  6. Katherine E Wilson1,
  7. Nick Reed1,4,5
  1. 1 Concussion Centre, Bloorview Research Institute, Holland Bloorview Kids Rehabilitation Hospital, Toronto, Ontario, Canada
  2. 2 Faculty of Kinesiology and Physical Education, University of Toronto, Toronto, Ontario, Canada
  3. 3 Bloorview Research Institute, Holland Bloorview Kids Rehabilitation Hospital Toronto, Toronto, Ontario, Canada
  4. 4 Rehabilitation Sciences Institute, Faculty of Medicine, University of Toronto, Toronto, Ontario, Canada
  5. 5 Department of Occupational Science and Occupational Therapy, University of Toronto, Toronto, Ontario, Canada
  1. Correspondence to Dr Melissa Paniccia; melissa.paniccia{at}


Objectives Following youth concussion, objective physiological indicators are needed to corroborate changes in subjective clinical measures. The objectives of the current study were to: (1) explore the effect of concussion on heart rate variability (HRV) across days postinjury in youth athletes aged 13–18 years old, compared with healthy age-matched and sex-matched controls and (2) examine the relationship between postconcussion symptom domains (physical, cognitive, fatigue and emotional) and HRV.

Methods Prospective, longitudinal, case-control study (N=44). This study comprised 29 concussed athletes between the ages of 13 and 18 years old (21 females, 8 males) and 15 age-matched and sex-matched controls). All participants completed baseline testing, which included demographic information (age, sex, concussion history), self-reported concussion symptoms (Post-Concussion Symptom Inventory [PCSI]) and a 24-hour heart rate recording via the Polar RS800CX system. The PCSI and HRV were collected weekly while the participant was symptomatic and then 1, 3 and 6 months following symptom resolution. HRV variables included time and frequency domain measures. Data visualisations and mixed effects modelling were used to derive parsimonious models.

Results HRV increased across days postinjury. Concussion symptom domains (physical, cognitive, fatigue and emotional) all had a significant main effect on HRV; concussed participants who reported more symptoms had higher HRV compared with those who reported fewer symptoms. Visualisations of HRV depict the recovery trajectory as non-linear across time. No significant differences on HRV measures were found between concussed and control participants.

Conclusion These preliminary findings provide the foundation to understand the varied trajectory and relationship between objective physiological measures and subjective symptom reporting.

  • concussion
  • youth
  • heart rate variability
  • recovery
  • symptoms

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See:

Statistics from

What are the new findings?

  • This study contributes to the knowledge base in bridging traditional subjective measures with objective physiological outcomes in an understudied youth athlete population.

  • Concussed youth athletes demonstrated non-linear trends of change in heart rate variability (HRV) across the recovery trajectory.

  • HRV was significantly associated with postconcussion symptom domains (physical, cognitive, emotional and fatigue) along days postinjury.

  • This study demonstrated that changes in HRV were present irrespective of symptom status (symptomatic/asymptomatic) and injury status (concussed/control).

How might the findings impact clinical practice in near future?

  • Findings from this study are still in an early stage of interpretation and are not yet ready for clinical uptake.

  • Findings from this study do set the foundation in understanding that a potential stress response may be at play along the recovery trajectory following concussion in youth athletes.

  • Clinicians are encouraged to continue using an individualised approach to assessment and treatment, while being cognizant of the underlying stress response.


It is estimated that 4 million children and youth present to emergency departments worldwide with concussion.1–3 In the assessment of concussion, well-established protocols have included neuropsychological testing and subjective symptom reporting to gauge clinical recovery.2–5 However, their application to youth athletes in particular presents with challenges. Elucidating the aetiology of concussion symptoms is difficult because postconcussion symptoms are non-specific.6 7 They can be found across a variety of clinical populations, including depression and anxiety8 and psychosocial stressors.9–11 Previous studies have also shown that athletes may minimise their symptoms over concerns of delayed return to play.12 13 Thus, objective physiological measures are needed to elucidate a physiological recovery trajectory that may present differently from a clinical recovery trajectory.

HRV, a measure of beat-to-beat variation in heart rate, has emerged as a non-invasive tool for the exploration of autonomic nervous system function following concussion.14–16 The exploration of HRV and concussion in sport populations has not been extensive, with few studies investigating the university athlete,17–19 and no studies examining the youth athlete population. In a recent prospective matched control study of university athletes, some perturbations in HRV (LFnu, HFnu) were found to persist beyond symptom resolution, which illustrated the potential of prolonged autonomic nervous system (ANS) disturbance following concussion.19 Similarly, Gall et al 20 found that concussed athletes who reported recent symptom resolution still demonstrated abnormal HRV during exercise compared with healthy controls. This finding was replicated by a recent study, which demonstrated that psychological systems and sleep disturbance following concussion in athletes were resolved at the return-to-play time point; however, ANS disturbances were still present (ie, decreased HRV), beyond return-to-play.17 It has been cited that more substantial follow-up will be valuable in improving the understanding of the natural history of change in ANS function following concussion.14 17 Further, it is not clear if the same trends between subjective symptom reporting and objective HRV measures would manifest similarly in the youth population.

The current consensus statement on concussion in sport3 has identified the following major gaps in assessing the value of physiological measures: (1) lack of longitudinal designs and (2) careful correlation with clinical measures, with little known on how these trends manifest within a youth population. Thus, this study aims to fill these gaps with the following exploratory research questions: (1) what is the effect of concussion on HRV across days postinjury in youth athletes aged 13–18 years old? (2) What is the relationship between postconcussion symptom domains (physical, cognitive, fatigue and emotional) and HRV? (3) What is the effect of having a concussion on HRV compared with healthy controls? As such, the hypotheses are the following: (1) HRV will decline immediately post injury and increase along the recovery trajectory; (2) HRV will have a negative relationship with postconcussion symptoms such that higher symptom reporting will be reflected in lower HRV and (3) concussed participants will demonstrate a decrease in HRV following concussion whereas controls will not demonstrate change compared with their baseline.


All participants and their parents provided informed consent prior to their participation in this study.


A convenience sample of 553 youth athletes (ages 13–18 years old) was recruited from local sport community organisations. Participants were part of a larger preinjury baseline testing study and completed a multimodal assessment consisting of physical, neurocognitive, subjective symptom reporting and physiological (HRV) measures (Reed et al, 2014).21 This study solely focused on the relationship between subjective concussion symptom reporting and objective HRV measures. Inclusion criteria were the following: male or female; ages 13–18 years old. Participants with a diagnosis of developmental delay, neurological condition, cardiac disease, experiencing symptoms from a previous concussion or non-English speaking were excluded from the study.


Demographic collection form

Age, sex and concussion history (ie, number of previous concussions) were collected using a questionnaire. Type of primary sport (ie, sport played most often and at the highest level) and level of play (eg, representative/competitive, house-league or recreational) were also collected.

Acute Concussion Evaluation (ACE) form

The ACE form is a postinjury questionnaire on injury characteristics, symptoms, risk factors for prolonged recovery (eg, history of migraine) and red flags that would warrant going to an emergency department.22 This measure was only used for those participants who sustained a concussion.

Post-Concussion Symptom Inventory (PCSI)

The PCSI is a 22-item self-report questionnaire used to assess postconcussion symptoms with youth ages 13–18 years old. It captures severity of symptoms ranging from 0=‘not a problem’, 3=‘somewhat of a problem’ and 6=‘severe problem’. Summed domain scores were derived for the following four domains: physical domain (eg, headache, dizziness, blurred vision); cognitive domain (eg, trouble concentrating, confused, mentally foggy); emotional domain (eg, sad, irritable, nervous) and fatigue domain (eg, fatigue, drowsiness). The PCSI has been found to have good validity and reliability among this age group.23 Internal consistency has been found to be moderate to strong, with Cronbach’s alpha ranging from 0.79 to 0.93 for the symptom subscales.23

Heart rate variability (HRV)

Heart rate recording was obtained from the Polar RS800CX watch and chest strap, with a sampling rate of 1000 Hz (RS800CX; Polar Electro, Kemple, Finland). Table 1 provides a definition of the time and frequency domain measures used in this study.24 A 24-hour recording methodology was suitable to the research objectives of this study, as they are more ecologically valid.25

Table 1

HRV variable definitions


On receiving their informed consent, all participants completed a baseline/preinjury assessment21 that included the completion of the demographic collection form (age [years], sex, concussion history [yes/no], number of previous concussions, type of sport and level of competition in their primary sport), PCSI, and participants were fitted with a Polar RS800CX watch and chest strap (RS800CX; Polar Electro, Kemple, Finland). Due to the group-based testing environment of the larger baseline study, individual availability and school schedules, the start and stop time of the HRV data collection varied. All participants were instructed to carry out their usual daily activities, with the exception of removing the device (and putting it back on) if they went swimming, took a bath, or played a contact sport.

Youth athletes who sustained a concussion were instructed to receive a diagnosis from a physician. Concussed participants were then followed by research personnel with health science backgrounds (occupational therapy, kinesiology); (1) weekly follow-up assessment while symptomatic; (2) 1, 3 and 6 months postsymptom resolution (asymptomatic). The first follow-up included completion of the ACE form by the research personnel and the PCSI and HRV by the participant. Within this first follow-up, all concussed participants received concussion education on management strategies according to the evidence informed ‘Concussion & You’ programme.26 Here, the following domains of concussion management strategies were addressed: (1) energy conservation; (2) sleep hygiene; (3) hydration and nutrition and (4) gradual return to school and sport. All subsequent follow-up visits included the completion of these measures (except for the ACE form). Age (±6 months) and sex-matched control participants were recruited from the baseline cohort and the same measures were collected at each follow-up with the exception of the ACE form at the first follow-up visit.

Data analysis

HRV data cleaning and processing

The RHRV package within R: A Language and Environment for Statistical Computing was used to process and clean the HRV data file. The FilterNIHR function was used to filter ectopic heart beats by using an adaptive threshold for rejecting beats whose value exceeds the cumulative mean threshold. Second, standard values were obtained by identifying minimum (25 beats per minute) and maximum (200) heart beats per minute. Regarding frequency domain measures (ie, HF and HFnu), window frames were set to 300 s, 50% overlap; thus, the measures were calculated as the mean of successive 5 min epochs. The InterpolatedNIHR function (4 Hz) created an equidistant time series. Power spectrum density, via fast Fourier transform, was used to derive the frequency domain variables: HF (0.15–0.4 Hz); LF (0.04–0.15 Hz). Recordings that were <14 hours and not continuous were excluded from the analysis as the duration of the recording has been found to influence HRV indices.27 For the concussed group (n=29), length of recording was M=19.45 hours, SD=4.37. Within the control group (n=15), length of recording was M=22.90 hours, SD=3.26.

Statistical analysis

The R program was used for all statistical analyses. Generalised linear mixed models were employed to explore the effect of concussion on HRV, along the recovery trajectory. Parsimonious models generated for this study resulted in a random intercept model (fixed main effects). Participant ID was the random intercept, allowing each individual participant’s trend to vary along days postinjury which is appropriate given the different variance structures across groups.28 This approach accounts for the variability in the number of follow-up assessments per participant and the variable amount of time between those time points. In terms of number of follow-up assessments, concussed group: Mdn=3, range=1–9; control: Mdn=3, range=1–7. Formal goodness of fit statistics were tested for each outcome variable with Akaike Information Criterion (AIC), in which lower values represent better fit. Sex, age and history of concussion were accounted for in these models. However, inclusion of those covariates in the models resulted in higher AIC values (lower fit) and omitted to satisfy the parsimony criteria. Baseline HRV variables were also accounted for in the models as variation exists among healthy individuals.29 The statistical threshold for significance was set to p<0.05.


A total of 44 participants were included in the analysis (concussed=29; control=15). Concussed participants and matched controls were similar on all HRV variables as well as PCSI values at baseline. As well, there were no significant sex differences at baseline on HRV variables, with the exception that females had lower HFnu compared with males. Descriptives for HRV values for concussed and control participants at baseline can be found in table 2. In terms of PCSI baseline values, females reported significantly more physical symptoms, males reported significantly more fatigue symptoms and females reported more overall symptoms compared with males. Participant demographic information and injury characteristics can be found in table 3. Finally, given the variable trajectory in number of follow-ups and time between follow-ups, table 4 depicts the number of HRV measures across time since injury.

Table 2

Descriptive HRV values for concussed (n=29) and control (n=15) participants at baseline

Table 3

Participant demographic information and injury characteristics

Table 4

Number of assessments according to HRV measure, across time since injury


No significant differences in SDNN were found when considering days postinjury, postconcussion symptom domains (physical, cognitive, emotional and fatigue), symptom status (symptomatic/asymptomatic), and concussed participants were not significantly different from controls (table 5).

Table 5

Mixed effect model estimates


A decrease in RMSSD appeared to occur 15 days postinjury and decreased until day 30, followed by levelling off at approximately 50 days (B=0.0001, p=0.02; figure 1). Postconcussion symptom domains, symptom status and injury status were not significant injury variables in this model (table 5).

Figure 1

Relationship between RMSSD, concussion symptom domain and days postconcussion where day 0 marks the day of first follow-up. The green and purple trend line represents the lowess smooth line, which maps an average trajectory for all concussed individuals.


Mixed modelling revealed a main effect of days postinjury (B =0.0009, p=0.029), whereby with increasing days post injury, pNN50 increased. Here, the visual trajectory of this finding reflected an immediate decrease from day of injury, decreasing until 30 days postinjury. PNN50 values then increase until day 65, followed by a plateau (figure 2). A positive, significant relationship was found between the physical (B =0.016, p=0.019), cognitive (B =0.023, p=0.012) and fatigue (B =0.028, p=0.030) domains and pNN50, whereby concussed participants who reported more symptoms in each domain were found to have increased pNN50. There were no significant effects of the emotional symptom domain, symptom status and injury status (table 5).

Figure 2

Relationship between pNN50, concussion symptom domain and days postconcussion where day 0 marks the day of first follow-up. The green and purple trend line represents the lowess smooth line, which maps an average trajectory for all concussed individuals.


Mixed modelling revealed a main effect of days postinjury (B =0.001, p=0.005), whereby concussed participants were found to have increasing values of HF as days postinjury increased. All four postconcussion symptom domains, physical (B =0.064, p<0.001), cognitive (B =0.138, p<0.001), fatigue (B =0.131, p<0.001) and emotional (B =0.305, p<0.001) were also found to have significant main effect on HF, whereby those who reported more symptoms in each of those domains displayed increased HRV. A main effect of concussion status was also found, in which concussed participants were found to have significantly higher HF compared with controls (B =0.101, p<0.001). Within the concussed participants, symptomatic vs asymptomatic participants did not significantly differ in their values of HF (table 5).


Mixed modelling revealed a main effect of days postinjury (B =0.0008, p<0.001), whereby concussed participants were found to have increasing values of HFnu as days postinjury increased. All four postconcussion symptom domains were also found to have a positive and significant main effect on HFnu. With increased report of concussion symptoms in the physical (B =0.014, p<0.001), cognitive (B =0.033, p<0.001), emotional (B =0.069, p<0.001) and fatigue domain (B =0.022, p<0.001), significant increases in HFnu were found. Within the concussed participants, symptomatic vs asymptomatic participants did not significantly differ in their values of HFnu (table 5).


Findings from this study suggest a general trend of increasing HRV along the recovery trajectory. The relationship between concussion symptoms and HRV revealed differential trends along days postinjury, which reflects a non-linear physiological pattern in concussion recovery. Age and sex differences were not present in this study; however, the current study’s sample was predominantly female. History of concussion did not appear to play a role in the findings. With a larger sample size and adequate control, these potential demographic differences may be elucidated in future studies.

HRV along days postinjury

This study found a general increase over days postinjury in concussed youth athletes across all HRV variables except SDNN. However, in the initial 30–40 days of concussion, HRV is shown to decrease (14%–25% decrease was observed in this time period). Following this period, increases are observed until approximately day 75 (8%–26% increase across HRV measures) and returns to baseline thereafter. This trend should be interpreted with caution, as there are fewer data as days postinjury increases. These findings are somewhat in contrast to a prospective cohort study of varsity athletes, which revealed significant differences in the acute phase of concussion (days 1–7), but these were only found in RMSSD compared with other HRV measures.30 It is important to note that looking at HRV during an isolated period of time (acute injury) may depict changes in HRV as linear, that is, an expectation that HRV may be linearly decreasing or increasing. The strength of this repeated measures study is in the visualisation where the trajectory of time revealed trends of decreasing HRV, followed by trends of increasing HRV. While the mechanism of change is still unknown, it is hypothesised that an ‘uncoupling’ exists between the ANS and cardiovascular systems whereby dysregulated shifts occur between the parasympathetic and sympathetic branches of the ANS.31 Finally, no real differences were found when comparing concussed youth athletes to controls and further investigation is needed to explore this difference.

HRV and postconcussion symptom domains

A consistent result was found across some of the HRV measures (pNN50, HF, HFnu) whereby those who reported more concussion symptoms displayed increased HRV and those who reported fewer concussion symptoms had lower HRV. This finding differs significantly from research to date. In concussed adult samples, physiological markers of increased sympathetic output have been found (eg, decreased HRV).32–34 As well, increases in symptom reporting or psychological stress have been associated with decreased HRV.17 19 35 It is worth noting that these cited studies were conducted in more controlled and short-term data collection environments (eg, assessing change in HRV following positional changes); thus, their direct comparison to this study warrants caution. Nonetheless, a potential hypothesis for this opposite trend includes the role of consistent clinical education follow-ups within this study. All concussed participants received one-on-one concussion management strategies regarding graduated return to activity and symptom limited activity. It may be possible that this served as a clinically meaningful interaction for the concussed youth athletes. For example, those with increased symptoms may have followed the management strategy to rest and graduate their level of activity. This in turn may have served as modulatory feedback to the ANS and influenced homeostasis (increased HRV).18 The implication of this novel finding is that changes in HRV may have the potential to be altered based on support and consistent clinical education. While this was not the intent of the study protocol, it behoves the concussion community to further investigate external factors, which may alter the recovery trajectory for youth athletes.


Overall, there appears to be large within-group variability in examining HRV following concussion. It has been postulated that autonomic dysregulation may present as a subclinical phenomena for those who have sustained a concussion36 and a physiological threshold may exist for those differences to be detectable.25 It is still unclear if the findings in this study represent a typical response to a stressful event (having a concussion) or if autonomic dysregulation is occurring above and beyond this response. Within the context of a large preinjury baseline study, it was not feasible to collect other meaningful variables such as physical activity repertoires. This is important to consider as changes in exercise training (frequency, duration and intensity) have been shown to significantly alter vagal tone (increase HRV over time).27 37

This study employed a non-controlled, 24-hour recording protocol, in which participants were able to go about their daily lives. This is in contrast to the majority of HRV concussion studies, which employ resting state position protocols or observe change in HRV according to exercise exertion. Restrictions on physical activity in this study were not feasible or appropriate for long-term recordings, given the diverse activity repertoires of youth (school, sport, extracurricular activity). This is an important consideration as changes in exercise training have been shown to significantly alter vagal tone (increase HRV over time).27 37 In the present study, the physical activity levels of concussed and control participants were likely different due to concussed participants’ inability to engage in physical activity and due to the training phase of each youth (ie, competitive, training or off season). Specific qualifications to physical activity (ie, frequency, duration and intensity) across the various follow-up time points will be important to consider for future study and likely will contribute valuable information to assist in the interpretation of change in HRV.

Other factors related to cognitive load and sleep-wake cycles are also important to consider. The role of cognitive demand/activity in everyday life (ie, academic demands, extracurricular activities) was not considered in this study. This study did not collect information in the form of a cognitive diary or employ other forms of cognitive activity tracking, and it is important to consider these clinical measures to allow a fulsome interpretation of the autonomic nervous system. Further, a potentially confounding effect may have been the inability to control for sleep-wake cycles. Sleep factors can significantly impact how an individual carries out their daily activities and can serve as a feedback loop to the ANS to regulate stress responses.38 Thus, a standardised assessment (ie, sleep questionnaire) to quantify these sleep factors is needed.

According to the most recent consensus statement, a single ‘physiological time window’ for concussion recovery does not exist due to methodological differences in assessing physiological measures and study design.3 Multiple studies have suggested physiological dysfunction that outlasts clinical measures of recovery,33 34 39 but this has not been conclusively established. The challenge present in this study was the inability to directly compare results to other similar concussion studies given the difference in recording protocol. However, regarding the baseline values, normative studies in adolescents examining 24-hour recordings show similarity across the measures collected in this study. For example, Faulkner, Hathaway, and Tolley40 examined healthy adolescents (15 years old ±1.6) and reported similar baseline values and ranges on pNN50, RMSSD, SDNN. Nonetheless, the finding that baseline values between concussed and control groups were similar and that baseline values compare to analogous studies may indicate promise in using the 24-hour recording as an ecologically valid methodology to observe change in the ANS. This assessment protocol may be the initial first step in using a dynamic protocol over a sufficient length of time to gauge how the ANS fluctuates in response to concussion over time. Taken together, while this 24-hour protocol introduces variability and noise within the HRV signal, the ability of trends to be seen despite this, is indicative of a potentially salient signal, and one that may be more ecologically valid41 42 when gauging the physiological response to concussive injury.

Finally, while robust hierarchical modelling was employed in this study to address the variability in data points across time since injury, the data were heavily concentrated within acute periods following concussion and more sparse within the chronic period. The power of the analysis approach to detect reliable differences may be influenced and thus skewed to reflect earlier time points in recovery. However, the authors present the visualisation of these recovery trajectories in an effort to be transparent about the distribution of data.


The authors thank the Rehabilitations Sciences Institute at the University of Toronto for their support throughout the doctoral program. We would also like to acknowledge the efforts of the members of the CIHR ‘NeuroCare’ Team and the members of the Concussion Centre (Bloorview Research Institute), specifically Talia Dick, James Murphy and Katherine Mah. We are grateful to the youth and families for their participation in this research. Jordan Collins is also thanked for his work and expertise in statistical analysis.


  1. 1.
  2. 2.
  3. 3.
  4. 4.
  5. 5.
  6. 6.
  7. 7.
  8. 8.
  9. 9.
  10. 10.
  11. 11.
  12. 12.
  13. 13.
  14. 14.
  15. 15.
  16. 16.
  17. 17.
  18. 18.
  19. 19.
  20. 20.
  21. 21.
  22. 22.
  23. 23.
  24. 24.
  25. 25.
  26. 26.
  27. 27.
  28. 28.
  29. 29.
  30. 30.
  31. 31.
  32. 32.
  33. 33.
  34. 34.
  35. 35.
  36. 36.
  37. 37.
  38. 38.
  39. 39.
  40. 40.
  41. 41.
  42. 42.
View Abstract


  • Contributors MP conceptualised the study objectives, data collection and analysis and drafted manuscript. LV contributed to the conceptualisation of the study with particular expertise in methodology and data processing, contributed to literature review and manuscript draft. SGT participated in the methodology and data processing design and editing manuscript. TT participated in the conceptual design with expertise in data processing, editing manuscript. MK conceptualised study objectives and edited manuscript. KEW participated in data collection, literature review and editing manuscript. NR is the principal investigator on this manuscript with a broader role in supervising the proposal, data collection and analysis, methodology, diligent manuscript draft review and editing for submission.

  • Funding This is work was funded by the Canadian Institutes of Health Research (#127048), the Ontario Neurotrauma Foundation and the Ontario Brain Institute. The Ontario Brain Institute is an independent non-profit corporation, funded partially by the Ontario government.

  • Disclaimer This work was originally featured within the first author’s doctoral thesis, as part of a larger body of work on neurophysiological variation (Paniccia, 2018 (University of Toronto, published doctoral thesis). The opinions, results, and conclusions are those of the authors and no endorsement by the Ontario Brain Institute is intended or should be inferred.

  • Competing interests None declared.

  • Patient consent Not required.

  • Ethics approval Holland Bloorview Research Ethics Board.

  • Provenance and peer review Not commissioned; externally peer reviewed.

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.