INTRODUCTION
============

Direct or indirect exposure to actual or threatened death (e.g., via road traffic accidents, natural disasters, witnessing the death of others) or severe traumatic events (e.g., sexual violence/abuse, combat) may lead to posttraumatic stress disorder (PTSD). PTSD affects approximately 6.8--12% of the US adult population at some point in their lifetime \[[@b1-pi-2019-0112]\]. Symptoms of PTSD (e.g. intrusive memories, re-experience and avoidance) have negative impacts on both physical (e.g., digestive, respiratory and cardiovascular diseases) \[[@b2-pi-2019-0112],[@b3-pi-2019-0112]\] and psychological (e.g., negative mood, recklessness, substance abuse and comorbidity with other mental illnesses) \[[@b4-pi-2019-0112]\] characteristics. PTSD not only affects quality of life \[[@b5-pi-2019-0112]\] but also increases the risk of suicide \[[@b6-pi-2019-0112]\]. Thus, PTSD is a serious public health problem that places heavy burdens on both individuals and society. However, PTSD is often not diagnosed and treated in a timely manner due to the overlapping symptoms of other psychiatric disorders and the tendency of clinicians to inadequately examine for symptoms of traumatic stress during clinical practice \[[@b7-pi-2019-0112]\].

The autonomic nervous system (ANS), which includes the sympathetic nervous system (SNS) and the parasympathetic nervous system (PNS), plays an important role in the onset and progression of PTSD \[[@b8-pi-2019-0112]\]. The SNS is considered to offer the most immediate response to a stressor via the reaction of the sympatho-adrenomedullary arm reacting \[[@b9-pi-2019-0112]\]. Then, the activation of the ANS declines quickly results in short-lived response, as the PNS is activated and acetylcholine is released into the system after the immediate stress events have passed. Thus, when the acute stressors experienced by individuals are converted to chronic or recurrent stress, the SNS system may be continuously activated without the confrontation of the PNS. Sustained hyperactivity of the SNS and abnormal of the PNS is accompanied by increasing catecholamine (epinephrine and norepinephrine) levels and decreasing acetylcholine in individuals with PTSD \[[@b10-pi-2019-0112]\]. The dysregulation of the ANS may link the physiological and psychological pathologies of PTSD \[[@b11-pi-2019-0112]\]. Heart rate variability (HRV) is the most commonly used measure of ANS activity because it is noninvasive as well as easy and inexpensive to monitor \[[@b12-pi-2019-0112]\].

Metrics of HRV are widely used in the field of behavioral and biomedical sciences. In 1996, the European Society for Cardiology and the North American Society of Pacing and Electrophysiology formally defined HRV as "heart rate oscillations within the beat-to-beat variability." \[[@b13-pi-2019-0112]\] Since then, an increasing number of studies have focused on short-term (minutes to 1 h) or long-term (24 h) HRV \[[@b14-pi-2019-0112]\]. HRV measures the state of the ANS based on the notion of balance between the ANS and SNS. Variability is considered a proxy for the health of a system, and decreased variability being a sign of a pathological condition \[[@b15-pi-2019-0112]\]. And HRV indices can be used to better understand the peripheral physiological correlates of brain and behavioral processes (e.g., emotion and cognition) \[[@b16-pi-2019-0112]\], possibly by reflecting the functionality of the prefrontal cortex \[[@b17-pi-2019-0112]\]. HRV indices include the following commonly used indices: low-frequency power (LF, 0.01--0.15 Hz) and high-frequency power (HF, 0.15--0.40 Hz) in both ms^2^ and normalized units, the ratio of LF to HF (LF/HF), the standard deviation of normal-to-normal (NN) intervals (SDNN) and the root mean square of successive differences between NN intervals (RMSSD). LF and HF are frequency domain indices, and SDNN and RMSSD are time domain (TD) indices. HF-HRV is usually considered a proxy for the PNS state and is similar to respiratory sinus arrhythmia (RSA) \[[@b18-pi-2019-0112]\]. LF-HRV is regulated by the SNS and PNS, and high LF reflects increased SNS activity \[[@b19-pi-2019-0112]\]. The ratio of LF/HF might reflect the vagal balance. However, recent studies have suggested that the LF/HF ratio does not have a mathematical basis and is not an appropriate index of sympathovagal balance \[[@b20-pi-2019-0112]\]. Among time domain metrics, the RMSSD has been thought to be the most appropriate method for quantifying HRV, and it closely represents PNS \[[@b21-pi-2019-0112],[@b22-pi-2019-0112]\].

Studies show that PTSD is associated with HRV. For example, the risk of cardiovascular disease is increased in PTSD \[[@b23-pi-2019-0112]\]. HRV is a reliable index of this risk disease \[[@b24-pi-2019-0112]\]. Dysfunction of emotion is associated with HRV \[[@b25-pi-2019-0112]-[@b27-pi-2019-0112]\]. Moreover, some brain regions (e.g., anterior cingulate cortex and medial prefrontal cortex) with abnormalities are associated with HRV \[[@b28-pi-2019-0112]\]. However, findings regarding the alterations in HRV in patients with PTSD are inconsistent. Regarding TD, research has revealed that PTSD may be related to decreased RMSSD and SDNN \[[@b29-pi-2019-0112],[@b30-pi-2019-0112]\], although not consistently so \[[@b31-pi-2019-0112]\]. There is some evidence suggesting that HF is lower among individuals with PTSD \[[@b32-pi-2019-0112],[@b33-pi-2019-0112]\]; however, some studies have not reported this effect \[[@b34-pi-2019-0112]\], and others have only reported TD and not frequency domain indices \[[@b31-pi-2019-0112]\]. Some studies have reported elevated LF in individuals with PTSD \[[@b30-pi-2019-0112],[@b35-pi-2019-0112]\], although inconsistent results have been reported \[[@b36-pi-2019-0112],[@b37-pi-2019-0112]\]. The earliest review of the relationship between PTSD and HRV was published in 2007 by Friedman \[[@b38-pi-2019-0112]\]. A subsequent 2014 meta-analysis of HRV alterations with anxiety found that HF-HRV was reduced in PTSD patients relative to healthy controls, whereas LF-HRV did not differ significantly \[[@b39-pi-2019-0112]\]. A limitation of the above studies is that they did not consider potential confounding factors of HRV in PTSD patients. Given the importance of increasing statistical power and reducing false positive rates, an up-to-date meta-analysis of studies of HRV and PTSD is needed.

We aimed to summarize the current studies comparing HRV between PTSD patients and healthy controls. We concentrated on quantifying the frequency domain and TD HRV indices. Furthermore, we investigated the source of heterogeneity (e.g. clinical, methodological), and analyzed potential variables via meta-regression.

METHODS
=======

Search strategy
---------------

The meta-analysis was conducted according to the PRISMA guidelines and registered in the PROSPERO database (registration number CRD42018099682). At the time of registration, no similar reviews were listed in the PROSPERO database. Four databases (PubMed, EMBASE, Medline, and Web of Science) were searched for the terms "heart rate variability," "HRV," "vagal," "automatic nervous system," "ANS," "posttraumatic stress disorder," "PTSD," and "post-traumatic stress disorder." Articles were also obtained by reviewing cited references. We retrieved on those articles published after 1996 as the guidelines for HRV testing and calculations were published in 1996. The article search was completed in August 2018. The article inclusion criteria were as follows: 1) a definitive diagnosis of PTSD as defined by the Diagnostic and Statistical Manual of Mental Disorders (DSM-III), DSM-IV, DSM-5 or International Classification of Diseases (ICD-10) and comparisons with a healthy control group without traumatic events and psychiatric disorder; 2) HRV index recording and calculation in accordance with the guidelines of the Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology; 3) resting state HRV reported; and 4) study published in English. Two authors independently reviewed all potentially relevant manuscripts, and cases of disagreement were resolved via consultation with the third investigator.

Data extraction
---------------

Using only the frequency domain or time domain HRV to assess the patients with PTSD has limitations \[[@b40-pi-2019-0112]\], so frequency domain and time domain measures were included in this meta-analysis. If more than one TD measure was reported in a study, RMSSD was prioritized, followed by SDNN. We chose to report only RMSSD in this review because RMSSD and SDNN offer identical information about the HRV and the former is more stable \[[@b41-pi-2019-0112]\]. Regarding the frequency domain, HF-HRV and LF-HRV measures are the most commonly reported measures in the PTSD literature. Thus, the frequency domain measures extracted in this study were HF-HRV and LF-HRV. If the data were presented only in graphs in an article, we obtained the original data via email with the corresponding author. Where longitudinal data were reported, we extracted only the baseline data to reduce the influences of confounding factors. If HRV parameters were reported separately for day and night in a study, we extracted only those for night. Nighttime HRV is more stable than daytime HRV in reflecting the activity of the vagal nerve due to less body movement and fewer confounding factors. For the meta-regression and subgroup analyses, the following data were also extracted: year of publication, country, mean participant age, quality score of the study, the score of the depression scale, proportion of males, source of participants, recording duration and recording values.

Quality assessment
------------------

In 2009, Tak et al. \[[@b42-pi-2019-0112]\] developed a scale to assess HRV in methodological studies. The instrument evaluates study quality based on the following three dimensions: appropriate selection of participants, appropriate quantification of HRV and appropriate control for confounders. Each dimension includes three questions. The total score of this scale ranges from 0 to 18 points, with a higher score indicating higher quality.

Meta-analysis statistics
------------------------

### Effect size estimation

Meta-analytic computations were performed using Review Manager 5.1 <https://community.cochrane.org/help/tools-and-software/revman-5>) \[[@b43-pi-2019-0112]\]. A random effects (RE) model was used in this study. To investigate the magnitude of the effect size of each variable across the different studies, we calculated the standardized mean difference (SMD, excluding the influence of the unit). The significance of Hedges' g was determined using 95% confidence intervals (95% CIs). When the 95% CI did not overlap the null hypothesis value of zero, the group difference was considered significant at the α=0.05 level. In one study, only the maximum and minimum values were reported \[[@b44-pi-2019-0112]\]. We calculated the standard deviation according to the Cochrane guide \[[@b43-pi-2019-0112]\]. Heterogeneity was assessed using I^2^ values, with I^2^≥25%, I^2^≥50%, and I^2^≥75% interpreted as low, medium and high heterogeneity, respectively \[[@b45-pi-2019-0112]\].

### Meta-regression and subgroup analyses

Regression analysis and subgroup analysis were performed using Stata 15.0 (StataCrop LP, College Station, TX, USA) and Review Manager 5.1 (<https://community.cochrane.org/help/tools-and-software/revman-5>), respectively. When significant heterogeneity was observed (I^2^≥50%), we use mixed mode regression analysis and the restricted maximum-likelihood estimation (RMLE) method to explore the source of heterogeneity among studies. The potential variables investigated were year of publication, mean participant age, study quality score, the score of depression scale and the proportion of male participants. Most studies use the Beck Depression Scale (BDI) to assess depression, Ray's \[[@b46-pi-2019-0112]\] study use the Patient Health Questionnaire (PHQ-9), and Cohen's \[[@b33-pi-2019-0112]\] study used the Hamilton Depression Scale (HAMD). Some of the variables of interest were limited to a small number of studies (n\<10); in these cases, subgroup analyses were conducted. These variables included the country of study (Asia, Europe and America), source of participants (clinical, community and clinical & community), recording duration (short term and long term) and recording values (raw and log-transformed). To explore whether the effects of the HRV indices in PTSD patients were epiphenomenal to medication use (e.g., use of antidepressants, benzodiazepines, or antipsychotics), we analyzed studies that excluded current medicine users.

RESULTS
=======

Inclusion of studies
--------------------

According to the search terms, 499 articles were compiled from the four databases. Among these, 206 duplicates, 3 case reports, and 12 reviews were excluded. After reading the titles and abstracts, the remaining 47 articles were potentially eligible for inclusion. After reading the full texts of these 47 articles, 19 articles were finally included ([Figure 1](#f1-pi-2019-0112){ref-type="fig"}).

Study and sample characteristics
--------------------------------

The characteristics of the 19 studies that compared resting state HF, LF, and TD HRV measures between participants who met the inclusion criteria are summarized in [Table 1](#t1-pi-2019-0112){ref-type="table"}. Of the 19 studies included, 15 reported HF, 12 reported LF, and 10 studies reported RMSSD. 16 studies compared HRV between 2 groups (PTSD vs health controls without traumatic events) and 4 studies \[[@b35-pi-2019-0112],[@b36-pi-2019-0112],[@b47-pi-2019-0112],[@b48-pi-2019-0112]\] compared HRV between 3 groups (PTSD vs. healthy controls without traumatic events and healthy controls with traumatic events).

Main effects
------------

The meta-analysis revealed a significant effect of PTSD status on HF-HRV (Z=4.18, p\<0.0001), with individuals with PTSD (n=435) having lower HF-HRV than healthy controls \[n=940; g=-1.58 (-2.32, -0.84); I^2^=96%; k=14\]. [Figure 2](#f2-pi-2019-0112){ref-type="fig"} shows the HF-HRV results. Visual inspection of the funnel plot demonstrates only slight asymmetry, suggesting that nonsignificant studies were lacking ([Figure 3](#f3-pi-2019-0112){ref-type="fig"}). Next, to assess possible bias, we calculated the main effect after excluding extreme SMD values (i.e., \>2) and still found a significant effect \[Z=2.20, p=0.03; g=-0.31 (-0.59, -0.33); I^2^=67%; k=10\]. We also calculated the effect omitting the study that only reported RSA. Again, a significant effect was found \[Z=4.05, p\<0.0001; g=-1.65 (-2.45, -0.85); I^2^=96%; k=13\]. To explore whether the reduction in HRV in individuals with PTSD relative to controls was affected by medication use (including antidepressants, benzodiazepines, or antipsychotics), we repeated the calculation excluding medication users. There was a significant decrease in HF-HRV among unmedicated PTSD patients relative to unmedicated healthy controls \[Z=2.59, p=0.010; g=-1.60 (-2.82, -0.39); I^2^=96%; k=5\].

The meta-analysis with LF-HRV as the dependent variable did not reveal a significant effect of PTSD status (Z=1.42, p=0.15), indicating that there was no significant difference in LF between individuals with PTSD (n=403) and healthy controls \[n=915; g=-0.63 (-1.49, 0.24); I^2^=97%; k=12\]. We then calculated the main effect after excluding extreme SMD values (i.e., \>2) and again found no significant effect \[Z=1.97, p=0.05; g=-0.28 (-0.56, -0.00); I^2^=64%; k=8\]. Five of thirteen studies excluded medication users. Among these unmedicated PTSD patients (n=161) and healthy controls (n=329), there was no significant effect of PTSD status on LF-HRV \[Z=0.28, p=0.78; g=0.32 (-1.97, 2.61); I^2^=98%; k=4\].

When RMSSD-HRV was included as the dependent variable, a significant effect of PTSD status was observed (Z=4.80, p\<0.00001), with PTSD patients (n=304) having lower RMSSD than healthy controls \[n=325; g=-1.96 (-2.76, -1.16); I^2^=94%; k=9\]. The details of each study and SMD are shown in [Figure 4](#f4-pi-2019-0112){ref-type="fig"}. Visual observation of the funnel plot ([Figure 5](#f5-pi-2019-0112){ref-type="fig"}) reveals a lack of data in the lower right of the plot. We then calculated the effect of PTSD status after excluding studies with extreme SMD values (i.e., \>2); a significant effect was found \[Z=5.51, p\<0.00001; g=-0.89 (-1.20, -0.57); I^2^=42%; k=5\]. After excluding studies that did not report standard deviations, a significant effect was again found, suggesting that the declined vagal activity in PTSD was robust \[Z=4.30, p\<0.0001; g=-1.80 (-2.62, -0.98); I^2^=94%; k=8\]. Four of ten studies excluded medication users. Among the unmedicated PTSD patients (n=77) and healthy controls (n=82), a significant effect of PTSD status on RMSSD-HRV was observed \[Z=3.51, p=0.0005; g=-1.44 (-2.24, -0.64); I^2^=74%; k=3\].

The results of the meta-regression are presented in [Table 2](#t2-pi-2019-0112){ref-type="table"}. Year of publication, study quality score, mean participant age and proportion of males did not show significant heterogeneity (p\>0.05). The results of the subgroup analyses are shown in [Table 3](#t3-pi-2019-0112){ref-type="table"}. When we evaluated the effects of PTSD status for different countries separately, and varied results were obtained. A significant effect of PTSD status on RMSSD-HRV was found for three country groups; however, the effect on HF-HRV was not significant for studies in America. When evaluating the effects for different sources of participants (clinical, community and clinical & community), a significant effect of PTSD status on RMSSD-HRV was observed for all three groups; however, the effect in the clinical populations was large, whereas that in the community populations was medium. The effect of PTSD status on HF-HRV was significant for the clinical populations but not the community and clinical & community populations. The calculation of methodological heterogeneity demonstrated that all but one HRV indices subgroup demonstrated a significant effect. The PTSD participants showed lower HRV indices than did healthy controls regardless of the applied methods except in the case of HF-HRV, where no significant effect of PTSD status among studies applying long-time recording was found.

DISCUSSION
==========

This meta-analysis analyzed existing studies on HRV in PTSD. When we investigated the HRV indices, we found reductions in HF-HRV (but not LF-HRV) and RMSSD. This finding suggests that HRV alterations in PTSD are nuanced and that some indices of HRV have implications for the diagnosis and treatment of PTSD. To evaluate the robustness of our findings, we calculated differences in HRV in different subgroups, including the source of participants, methodology and medication-use subgroups. In the meta-regression analyses, there was no significant influence of the variables of interest on the resting HRV (HF-HRV and RMSSD) difference between the PTSD participants and healthy controls. The subgroup analyses revealed subgroup differences for HF-HRV: no significant reduction in HF-HRV in PTSD participants relative to controls was found in the American subgroup, in the community and clinical & community subgroups, and in the long-term recording subgroup. Only one of the subgroup analyses revealed variation in subgroups in the difference for RMSSD: no significant decline in RMSSD in individuals with PTSD relative to healthy controls was observed in the community subgroup. This finding suggests that RMSSD has greater robustness than HF as an HRV measure. However, given the small number of studies evaluated, future research is needed to confirm the reliability of these findings. Among the studies excluding medicated individuals, individuals with PTSD had lower HF-HRV and RMSSD than healthy controls. Thus, it is unlikely that the reductions in HF and RMSSD in PTSD individuals can be fully attributed to the intake of antipsychotic or benzodiazepine medications.

Effects of PTSD on the frequency and time domains
-------------------------------------------------

We found reductions in HF-HRV and RMSSD in PTSD, with a more robust reduction in the latter being more robust. Generally, high HRV reflects sympathovagal balance or well-conditioned adaptability of the ANS, whereas low HRV is a sign of deficiency and abnormal regulation of the ANS \[[@b49-pi-2019-0112]\]. PTSD is associated with the withdrawal of parasympathetic cardiac control, which is typically observed in the HF domain \[[@b29-pi-2019-0112]\]. The current results are consistent with those of the Hauschildt study \[[@b29-pi-2019-0112]\], which showed that a lower HF-HRV is significantly associated with PTSD. In TD analyses, RMSSD is used a measure of cardiac change and reflects the overall degree of ANS activity \[[@b50-pi-2019-0112]\]. The significantly reduced RMSSD in individuals with PTSD relative to that in healthy controls might also result from the decline in PNS function due to trauma events in PTSD individuals. The results are in accordance with the theory of neurovisceral integration, in which low PNS function is an indicator of impaired central-peripheral neural feedback. According to this theory, when an individual with low PNS function encounters environmental changes, he or she lacks sufficient psychophysiological resources to deal with them \[[@b27-pi-2019-0112]\]. Evidence suggests that individuals with high HRV at rest have higher scores on trait resilience psychometric scales \[[@b51-pi-2019-0112]\]. Individual with PTSD may continuously perceive threats and experience fear due to disinhibition of the threat circuits \[[@b52-pi-2019-0112],[@b53-pi-2019-0112]\].

The PNS can release acetylcholine to modulate the immune system and stimulate the cholinergic anti-inflammatory pathway to inhibit cytokine production \[[@b54-pi-2019-0112],[@b55-pi-2019-0112]\]. Accordingly, increased inflammation in patients with PTSD is observed. A recent meta-analysis reported that IL-1β, IL-6, TNF-α and IFN-γ levels were increased in patients with PTSD relative to healthy controls \[[@b56-pi-2019-0112]\]. In addition, neuroimaging studies have found that inflammation can alter brain structure and function. For example, IL-6 can modulate the connection between the amygdala and the anterior cingulate cortex (ACC) \[[@b23-pi-2019-0112]\]. Individuals with increased inflammation showed a stronger connection between the amygdala and the dorsomedial prefrontal cortex (PFC) than did individuals with normal levels of inflammation \[[@b57-pi-2019-0112]\]. These regions, as well as the ACC and media prefrontal cortex have been found to be associated with HRV \[[@b28-pi-2019-0112]\]. Thus, a reduction of HRV can be regarded as an endophenotype for PTSD, characterized by increased inflammation and alterations in brain structure and function. HRV as an endophenotype could play an important role in the diagnosis and treatment of PTSD and serve as a valuable biomarker in research.

Medicine use
------------

The present study suggests that the HRV decrease observed in PTSD is not completely attributable to medication use. In the studies of unmedicated individuals, we found that individuals with PTSD had lower HF-HRV and RMSSD than did healthy controls. In all of these studies, except one study that investigated drug-naive PTSD \[[@b36-pi-2019-0112]\], the participants had received antidepressant treatment at some time prior to the assessment. The washout period in these studies was approximately 4 weeks, and it is possible that the previous medication continued to have influence on the ANS at the time of the study. Few studies reporting whether HRV indices in patients with PTSD are influenced by medication, and their results are inconsistent. Cohen reported that HRV indices in patients with PTSD returned to normal after treatment with selective serotonin reuptake inhibitors (SSRIs) \[[@b58-pi-2019-0112]\]. However, other studies found that SSRIs do not improve HRV reductions in individuals with depressive symptoms \[[@b59-pi-2019-0112]\]. A review reported high dropout rates in SSRI trials and in trials with patients with PTSD \[[@b60-pi-2019-0112]\]. Participants in Wahbeh's study were consistently medicated for at least 4 weeks before assessment, and a lower peak frequency HF-HRV was observed in individuals with PTSD than in healthy controls \[[@b32-pi-2019-0112]\]. There is evidence that HRV reduction predates the development of PTSD. A large cohort study revealed that decreased HRV at baseline predicted PTSD at follow up \[[@b61-pi-2019-0112]\], indicating that HRV alterations might predate the development of PTSD. Therefore, it is possible that disturbance of the ANS may be a biomarker of the susceptibility to PTSD patients, not results.

Clinical implications
---------------------

Our results have a wide range of important diagnostic and treatment implications for PTSD. Decreased HRV might be considered as an endophenotype related to increased inflammation, alterations in brain structure and function, impaired executive dysfunction and emotion dysregulation in PTSD. Furthermore, HRV could serve as a valuable biomarker to monitor PTSD and treatment outcomes. An intervention that increases cardiac vagal tone may mitigate the imbalances in PTSD. A few intervention studies have aimed to increase HRV in PTSD patients via different approaches, with promising results obtained from relaxation training, mindfulness training \[[@b62-pi-2019-0112],[@b63-pi-2019-0112]\], and pharmacological treatment \[[@b58-pi-2019-0112]\]. In addition, individuals who show higher HRV during resting conditions seem to score higher on resilience questionnaires, to recover more efficiently from acute psychological stress and to be less vulnerable to the development of PTSD \[[@b65-pi-2019-0112]\]. Moreover, HRV reduction may predate PTSD. Therefore, we can intervene in susceptible populations by increasing HRV in advance to decrease the prevalence of PTSD.

Limitations
-----------

This study has several limitations. First, some studies had to be excluded because the results were given in the form of charts and we received no response from the contacted authors. However, all of these excluded studies revealed lower HRV for PTSD than for healthy controls. Thus, their exclusion would not change the results. Second, there was not a sufficient number of studies for a meta regression, and possible sources of bias need to be investigated in the future. Subgroup analysis was conducted across trials and not within trials; thus, we can conclude that the studied variables had no influence only at the study level. Third, because most of the articles did not report the trauma type or involved participants who had experienced different traumas, we were unable to explore HRV by trauma type.

We found a reduction in HRV in patients with PTSD relative to HRV in healthy controls, which cannot be completely attributed to medication use. Considering that HRV is potentially associated with central-peripheral integration, changes in inflammatory factors, PTSD-related brain regions and executive function, this appears to be a promising endophenotype for PTSD research.

This study was financially supported by the Special Project on Natural Chronic Non-infectious Diseases (grant nos. 2016YFC1307201); the National Natural Science Foundation of China (Grant Nos. 81701328, 81871061, and 81371484); the National Key Research & Development Program of China (grant nos. 2016YFC1307200); the China Postdoctoral Science Foundation (grant nos. 2017M612972) and the Postdoctoral Foundation of Sichuan University (grant nos. 2018SCU12042) to Dr. H.Z; and the Department of Science & Technology of Sichuan Province (grant nos. 2018 SZ0131).

The authors have no potential conflicts of interest to disclose.

**Author Contributions**

Conceptualization: Fenfen Ge. Data curation: Fenfen Ge, Ying Li. Formal analysis: Fenfen Ge. Funding acquisition: Wei Zhang. Investigation: Minlan Yuan. Methodology: Fenfen Ge. Project administration: Fenfen Ge. Resource: Wei Zhang. Software: Fenfen Ge. Supervision: Wei Zhang. Validation: Wei Zhang. Visualization: Fenfen Ge. Writing-original draft: Fenfen Ge. Writing-review & editing: Minlan Yuan, Wei Zhang

![Systematic search of the literature: PRISMA flowchart. PTSD: posttraumatic stress disorder, HRV: heart rate variability.](pi-2019-0112f1){#f1-pi-2019-0112}

![SMD in HF-HRV between individuals with PTSD and healthy controls. SMD: standard mean difference, HF-HRV: high-frequency heart rate variability, PTSD: posttraumatic stress disorder.](pi-2019-0112f2){#f2-pi-2019-0112}

![Meta-analysis funnel-plot comparing the high-frequency heart rate variability in individuals with posttraumatic stress disorder and healthy controls. SMD: standard mean difference.](pi-2019-0112f3){#f3-pi-2019-0112}

![SMD in RMSSD-HRV between individuals with PTSD and healthy controls. PTSD: posttraumatic stress disorder, SMD: standard mean difference, RMSSD-HRV: root mean square of successive R-R interval differences-heart rate variability.](pi-2019-0112f4){#f4-pi-2019-0112}

![Meta-analysis funnel-plot comparing the root mean square of successive R-R interval differences-heart rate variability in individuals with posttraumatic stress disorder and healthy controls. SMD: standard mean difference.](pi-2019-0112f5){#f5-pi-2019-0112}

###### 

Characteristics of included studies

  Study                   Age                                             Sample size (PTSD/controls)                    Quality   Trauma type   Medicine                                                                                                                                                         Exclusion criteria                                                                                                                       Length of recording                                                                                                                                                                                                       HRV indices   Conclusions    
  ----------------------- ----------------------------------------------- ---------------------------------------------- --------- ------------- ---------------------------------------------------------------------------------------------------------------------------------------------------------------- ---------------------------------------------------------------------------------------------------------------------------------------- ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- ------------- -------------- --------------------------------------------------------------------------------------------------
  Clinical                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
                          Cohen et al. \[[@b64-pi-2019-0112]\]            Matched                                        9/9       8             Combat, acts of terrorism, traffic accident                                                                                                                      Not reported                                                                                                                             Not reported                                                                                                                                                                                                              20 min ECG    HF (raw)       PTSD patients had significantly lower HRV compared to controls.
                          Cohen et al. \[[@b33-pi-2019-0112]\]            Matched                                        14/25     8             Combat, acts of terrorism, traffic accident                                                                                                                      Not reported                                                                                                                             Not reported                                                                                                                                                                                                              15 min ECG    HF (log)       PTSD has lower HRV values than controls.
  LF (log)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
                          Lakusic et al. \[[@b30-pi-2019-0112]\]          Matched                                        34/34     9             Combat                                                                                                                                                           Antidepressant (N=26)                                                                                                                    All: atrial fibrillation; sick sinus syndrome; second or third degree atrioventricular block, DM, CVD                                                                                                                     24 h ECG      HF (raw)       PTSD has lower HF, LF and RMSSD values than controls.
  LF (raw)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
  RMSSD (raw)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
                          Blechert et al. \[[@b66-pi-2019-0112]\]         Matched                                        23/32     11            Traffic or work-related accident, physical or sexual violence, war-related trauma, natural disasters, other trauma                                               Not use of benzodiazepines, blockers, sympathomimetic drugs, antipsychotics and tricyclic antidepressants                                History of psycho sis, drug abuse, a medical history of conditions that might affect the physiological systems under study and the                                                                                        5 min ECG     RSA(log)       PTSD has lower RSA than healthy controls.
                          Keary et al. \[[@b67-pi-2019-0112]\]            Matched                                        20/20     12            Not mentioned                                                                                                                                                    Not reported                                                                                                                             DM, CVD, BMI\>35                                                                                                                                                                                                          20 min ECG    HF (log)       LF and HF have no difference between the two groups.
  LF (log)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
                          Tan et al. \[[@b68-pi-2019-0112]\]              Matched                                        20/10     9             Combat                                                                                                                                                           Not reported                                                                                                                             Severe psychopathology, previous participation in another study involving HRV                                                                                                                                             15 min ECG    SDNN (raw)     PTSD has lower SDNN than healthy controls.
                          Shaikh al arab et al. \[[@b47-pi-2019-0112]\]   Matched                                        11/10     9             Traffic accident                                                                                                                                                 No psychiatric medication, benzo diazepinesanticonvulsants, antihypertensive drugs                                                       TBI, dementia, CVD, DM and history of substance abuse                                                                                                                                                                     24 h ECG      RMSSD (raw)    PTSD has lower RMSSD and SDNN than healthy controls.
  SDNN (raw)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
                          Slewa-Younan et al. \[[@b69-pi-2019-0112]\]     Matched                                        12/23     8             25 resettled Iraqi refugees                                                                                                                                      Cardiovascular and diabetes related drugs                                                                                                CVD, DM, neurological disorders, TBI and psychological disorders                                                                                                                                                          5 min ECG     HF (raw)       LF and HF have no difference between the two groups.
  LF (raw)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
                          Lee and Ttheus \[[@b48-pi-2019-0112]\]          Not reported                                   37/88     6             Combat                                                                                                                                                           Not reported                                                                                                                             All: male, no-veterans, less than four consecutive normal sinus rhythm R-R intervals                                                                                                                                      24 h ECG      RMSSD (raw)    PTSD has lower RMSSD and SDNN than healthy controls.
  SDNN (raw)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
                          Moon et al. \[[@b37-pi-2019-0112]\]             Not matched                                    34/27     8             Not mentioned                                                                                                                                                    Not reported                                                                                                                             Neurological disorder, substance abuse, mental retardation and TBI                                                                                                                                                        5 min ECG     HF (raw)       PTSD has lower HF than healthy controls. LF and RMSSD have no difference between the two groups.
  LF (raw)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
  RMSSD (raw)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
                          Chang et al. \[[@b36-pi-2019-0112]\]            Matched                                        32/192    12            Traffic accidents, natural disaster, human-made disaster, serious physical injury due to non-motor vehicle related accident; witnessd someone seriouly injured   Subjects who did not use any medication that have been reported to affect the ANS functioning for at least two weeks before evaluation   Pregnancy, smoking, DM, cancer, neuropathy, CVD, engaging in regular physical training exceeding 10 hours a week, history of substance dependence, organic brain disease or any concomitant major psychiatric disorders   5 min ECG     HF (raw)       PTSD has lower HF and LF than healthy controls.
  LF (raw)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
                          Thome et al. \[[@b70-pi-2019-0112]\]            Matched                                        57/41     9             Not mentioned                                                                                                                                                    No medication for mental illness and cardiovascular disease was used for one month prior to the study                                    TBI, neurological disorders, developmental disorders, pregnancy, mental illness, substance abuse                                                                                                                          6 min ECG     HF (log)       PTSD has lower HF, LF and RMSSD than healthy controls.
  LF (log)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
  RMSSD (log)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
                          Park et al. \[[@b71-pi-2019-0112]\]             Matched                                        45/43     12            Combat                                                                                                                                                           32 DM med                                                                                                                                All: mental disorder; current substance-use; uncontrolled hyperlipidemia, cerebrovascular disease and diabetes mellitus                                                                                                   15 min ECG    RMSSD (raw)    PTSD has lower RMSSD and HF than healthy controls.
  64 HD med               SDNN (raw)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              
  14 CVD med              HF (log)                                        LF has no difference between the two groups.                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
  LF (log)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
  Community                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       
                          Tucker et al. \[[@b35-pi-2019-0112]\]           Matched                                        13/32     10            Not mentioned                                                                                                                                                    No medication use                                                                                                                        Physical illness or, not meeting exposure criteria declining participation, and not making appointments                                                                                                                   5 min ECG     HF (raw)       PTSD has lower HF and LF than healthy controls.
  LF (raw)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
                          Agorastos et al.l \[[@b31-pi-2019-0112]\]       Matched                                        7/8       10            Combat                                                                                                                                                           Only one PTSD patients reported intake of antidepressants in the past                                                                    All: presence or history of any physical and mental comorbidities, previous or current abuse of/dependency on alcohol or other substances, smoking, BMI\>30, BMI\<18                                                      24 h ECG      RMSSD (raw)    RMSSD has no difference between the two groups.
  SDNN (raw)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      
                          Shah et al. \[[@b72-pi-2019-0112]\]             Not matched                                    14/25     11            Combat                                                                                                                                                           Antidepressant (51.6% PTSD)                                                                                                                                                                                                                                                                                                                                        24 h ECG      HF (log)       PTSD has lower LF and HF than healthy controls.
  LF (log)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        
  Clinical & Communtity                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           
                          Wahbeh and Oken \[[@b32-pi-2019-0112]\]         Matched                                        52/29     12            Combat                                                                                                                                                           Participants were required to be in good general health and on stable dosed of medications for at least 4 weeks prior to study onset     All: over 65 years old; mental disorder; substance dependence disorder within 3 months or current substance; two drinks or less of alcohol per day; sexual assault as primary PTSD events; TBI                            5 min ECG     HF (raw)       PTSD has lower LF than healthy controls.
  LF (raw)                HF has no difference between the two groups.                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
                          Meyer et al. \[[@b73-pi-2019-0112]\]            Not matched                                    18/23     10            Multiple traumata                                                                                                                                                No psychotropic medication                                                                                                               ALL:neurological disorders; current alcohol/drug abuse, medical illness, cardiac dysregulation                                                                                                                            5 min ECG     HF (raw)       PTSD has lower LF and RMSSD than healthy controls.
  LF (raw)                HF has no difference between the two groups.                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
  RMSSD (raw)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     
                          Ray et al. \[[@b46-pi-2019-0112]\]              Not matched                                    41/41     12            Combat                                                                                                                                                           No benzodiazepines and \|3-blockers                                                                                                      Veterans who were taking daytime benzodiazepines or \|3-blockers during initial screening were excluded                                                                                                                   5 min ECG     HF (log+raw)   PTSD has lower HF and RMSSD than healthy controls.
  RMSSD (raw)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     

DM: diabetes mellitus, CVD: cerebrovascular disease, MI: myocardial infarction, HP: hypertension, TBI: traumatic brain injury

###### 

Meta-regression analyses of the effects of potential variables

  Variable                          HF-HRV   RMSSD                                                                                           
  --------------------------------- -------- ------------------------- ------- ------- ------- ---- ------------------------ ------- ------- -------
  Year                              15       -0.109 (-0.272, -0.055)   0.084   -1.30   0.194   10   -0.281 (-0.697, 0.135)   0.212   -1.32   0.186
  Quality                           15       0.692 (-0.067, 1.451)     0.387   1.79    0.074   10   0.374 (-0.515, 1.263)    0.453   0.83    0.409
  Age                               13       -0.080 (-0.198, 0.038)    0.060   -1.33   0.184                                                 
  BDI                               10       0.016 (-0.186, 0.218)     0.103   0.16    0.874                                                 
  Proportion of male participants   13       -0.005 (-0.037, 0.028)    164.4   1.30    0.195   10   -0.019 (-0.054, 0.016)   0.018   -1.07   0.187

Stata 15.0 was used for meta-regression. HF-HRV: high-frequency heart rate variability, RMSSD: root mean square of successive R-R interval diﬀerences, BDI: Beck Depression Scale

###### 

Subgroup analyses of the effects of potential variables

                               HF-HRV   RMSSD-HRV                                                                        
  ---------------------------- -------- ----------- ---------------------- ----------- --- ------ ---------------------- ---------
  Country                                                                                                                
   Asia                        7        3.27        -2.17 (-3.47, -0.87)   0.0001      3   2.27   -2.62 (-4.88, -0.36)   0.0200
   Europe                      3        2.08        -1.37 (-2.66, -0.08)   0.0400      3   3.02   -2.50 (-4.11, -0.88)   0.0030
   America                     5        1.74        1.30 (-0.16, 2.76)     0.0800      4   4.00   -1.00 (-1.50, -0.51)   0.0001
  Source of participants                                                                                                 
   Clinical                    9        3.44        -1.92 (-3.02, -0.83)   \<0.00001   7   4.08   -2.39 (-3.53, -1.24)   0.00001
   Community                   6        1.19        0.76 (-0.49, 2.00)     0.2300      1   1.73   -0.96 (-2.05, 0.13)    0.0800
   Clinical & Communtiy        3        1.13        1.80 (-1.31, 4.90)     0.2600      2   5.33   -1.03 (-1.41, -0.65)   0.00001
  Recording duration                                                                                                     
   Short-time                  13       3.85        -1.81 (-2.74, -0.89)   0.0001      6   4.47   -2.36 (-3.39, -1.32)   0.00001
   Long-time Reported values   2        1.53        -0.27 (-0.61, 0.07)    0.1300      4   2.78   -1.15 (-1.96, -0.34)   0.0050
  Reported values                                                                                                        
   Raw                         8        2.73        -1.37 (-2.35, -0.39)   0.0060                                        
   Log transformed             7        2.99        -1.85 (-3.06, -0.64)   0.0030                                        

Review manager 5.1 was used for subgroup analysis. HF-HRV: high frequency-heart rate variability, RMSSD: root mean square of successive R-R interval diﬀerences, CI: confidence interval
