Study selection
Literature search identified a total of twenty-eight (28) studies that seemed relevant and were sought for full-text. Eight of these were excluded due to various reasons; Nakao et al. (2012)(13) used comorbid than healthy control; Goldstein et al. (2005)(14) did not assess our outcome of interest; Yoshiuchi et al. (2004)(15) used POTS participants comorbid with chronic fatigue syndrome; Singer et al. (2003)(16) intervened the control group with isoproterenol infusion similarly Freitas et al. (2000)(17) who intervened with cardio-selective beta-blocker and/or fludrocortisone and Stewart et al. (2007)(18) who employed hand-grip maneuver than HUTT. Furthermore, Bongiovanni et al. (2013)(19) and Aoki et al. (2008)(20) were excluded for not accessible full-text and use of Japanese language in the full-text retrieved, respectively. A total of twenty (20) studies fulfilled eligible criteria for inclusion. Figure 1, summarizes search results, screening and selection process.
Study characteristics
Table 1 summarizes the study characteristics of our twenty (20) studies that were eligible for inclusion in our study. Total number of participants reported was 1358, of these, 717 POTS and 641 were healthy participants. Regarding participants demographics, while other studies recruited both gender equally(21), other only recruited one gender participants (22), and other studies randomly involved both gender (23). While other studies matched the groups by age (14, 22), other studies did not (24). Furthermore, majority of studies reported participants’ ages central tendencies by mean, two studies utilized median instead (25, 26). While other studies used larger sample size (27), other used smaller sample sizes (28).
All twenty studies were case-control observational studies and none was interventional. These were conducted in different settings from a diverse number of countries all around the world. Eleven studies were done in USA, two in Australia and other were conducted in Israel (29), UK (23), Portugal (30), Japan (31), Germany (22), Korea republic (24) and The Netherlands (21), each contributing one study. This was thought to increase the external validity of this study.
Despite the fact that the search was not confined to any specified range of dates, none of the included study was found to have been published before the year 2000. Fifteen studies (75%) were published in the last decade.
Different studies reported different outcomes, but all aligned with our study questions. Eighteen studies compared heart rate, three studies compared RR-Interval, three studies compared rMSSD, six studies compared LF, eight studies compared HF, six compared LF/HF-Ratio and four studies compared LF(n.u) and HF(n.u) each. Studies comparing similar outcomes were analyzed together in the same forest-plot.
Table 1: Study characteristics
Study, Year
|
Study size
(POTS, Healthy)
|
Mean Aged
(POTS, Healthy)
|
Matched case-control or not?
|
Duration of HR/HRV parameter measurement (Angle of tilt)
|
Orthostasis induction method
|
Country of study
|
Outcome Recorded
|
Jacob 2019(29)
|
12,10
|
30±1.8, 32±3
|
Unmatched
|
30 Minutes (750)
|
HUTT
|
Israel
|
HF, HR
|
Owens 2018 (23)
|
20,20
|
36±10.84, 35±7.56
|
Unmatched
|
10 Minutes (600)
|
HUTT
|
UK
|
HF, LF
|
Goff 2017(32)
|
9,20
|
NA, NA
|
Unmatched
|
24 Hours
|
Daily life activity
|
Australia
|
rMSSD
|
Moon 2016(24)
|
46,67
|
28.9±1.9, 49.4±2.1
|
Unmatched
|
20 Minutes
|
Active standing
|
Korea Republic
|
HR
|
Freitas 2015(30)
|
10,12
|
29.4±8.5, 33.8±5.9
|
Matched
|
40 minutes (700)
|
HUTT
|
Portugal
|
HR, HF
|
Yoshida 2014(31)
|
70,38
|
13.7±0.1, 13.5±0.1
|
Unmatched
|
7 minutes (900)
|
Active standing
|
Japan
|
HR, LF/HF-Ratio
|
Medow 2014(25)
|
12,19
|
Median:
20.8, 21.4
|
Unmatched
|
10 Minutes (700)
|
HUTT
|
USA
|
HR, LF(n.u), LF(n.u)
|
Mallien 2014(22)
|
38,31
|
25.3±7, 26.2±6.3
|
Matched
|
Overnight
|
HUTT
|
Germany
|
HR, LF, HF, LF/HF-Ratio
|
Plash 2013(33)
|
15,15
|
36±3, 33±2
|
Unmatched
|
30 Minutes
|
Active standing
|
USA
|
HR
|
Ocon 2012(34)
|
16,20
|
21±1, 23±1
|
|
10 Minutes (750)
|
HUTT
|
USA
|
HR
|
Brewster 2012(35)
|
54,26
|
35±2, 27±1
|
Unmatched
|
5 Minutes
|
Active standing
|
USA
|
HR
|
Galbreath 2011(36)
|
17,17
|
27±9, 31±10
|
Unmatched
|
5 Minutes (600)
|
HUTT
|
USA
|
HR, HF, LF, rMSSD, RR-Interval, HF (n.u), LF (n.u)
|
Baumert 2011(37)
|
13,12
|
(32±13, 23±2),
|
Unmatched
|
10 Minutes (400)
|
HUTT
|
Australia
|
HR
|
Fu 2010(38)
|
27,16
|
26(21-33), 28(23,35)
|
Unmatched
|
45 Minutes (60%) *grip
|
HUTT
|
USA
|
HR
|
Ocon 2009(28)
|
9,7
|
NA, NA
|
Unmatched
|
10 Minutes (700)
|
HUTT
|
USA
|
RR-Interval, LF, HF, LF/HF-Ratio, LF(n.u), HF(n.u), HR
|
Garland 2007(39)
|
150,63
|
34.5±10.7, 30.2±9.3
|
Unmatched
|
5 Minutes
|
Active standing
|
USA
|
HR
|
Stewart 2006(40)
|
20,10
|
17±2, 17±1
|
Matched
|
10 Minutes (700)
|
HUTT
|
USA
|
HR, LF/HR-Ratio, HF, LF
|
Meier 2006(21)
|
21,39
|
15.5±2.2, 11.7±2.7
|
Unmatched
|
12 Minutes (600)
|
HUTT
|
The Netherlands
|
HR
|
Garland 2005(27)
|
136,191
|
29.1±8.0, 32.2 ±9.9
|
Unmatched
|
30 Minutes (600)
|
HUTT
|
USA
|
HR
|
Stewart 2000(26)
|
22,10
|
Median: 15.2, 15.8
|
Unmatched
|
30 Minutes (700)
|
HUTT
|
USA
|
HR, HF(n.u), LF(n.u), LF/HF-Ratio, HF, LF, rMSSD, RR-Interval
|
POTS - Postural orthostatic tachycardia; rMSSD - square root of mean of squares of successive R-R interval; LF - Low frequency power; HF - High frequency power; LF(n.u) - Low frequency power -normalized units; HF(n.u) - High frequency power -normalized units; HR – Heart Rate; NA – Data not accessed
Sources of bias
All 20 eligible articles included in this study were assessed for risk of bias in two levels; at study level and at review level. Regarding study level bias assessment; different studies involved different number sample sizes. Other studies included large number of participants (27) while other used low (28). It follows that, large sample sizes are more representative of general population as compared to small sample sized studies. Furthermore, none of these 20 studies reported to have had calculated the required sample size prior to their conduction.
Despite the fact that all studies were similar in that, they compared POTS versus healthy participants, some studies matched the comparison groups to reduce confounding factors while other studies did not (32). This might have introduced confounding factors to our study as factors such as female gender, BMI, physical fitness and race, each has been reported to independently alter heart rate variability (41).
Different studies utilized different methods to induce orthostasis, with other using Head-up Tilt test and others applying active standing (31).. Whether HUTT or active standing was used to induce orthostasis, , different durations, ranging from 5-45 minutes, were applied depending on the participant’s tolerance to orthostasis. Furthermore, different angles of tilt were set, ranging from 400 (37) to 750 in other studies. While majority involved awake patients, other studies (22) utilized sleeping participants. While other studies used ECG to measure heart rate variability in a short session (38), others used Holter device to record mean heart rate variability per 24 hours while a participants are carrying on with their daily activities(32). These different conditions were thought to increase heterogeneity hence influence our results.
At review level on the other hand, a number of loopholes for biases were also identified. Although, other studies had our data of interest, readily available to extract from tables in text, from one (23) study data had to be extracted by estimations and extrapolation from a graphical figure. This led to conducting sensitivity analysis excluding this study. Furthermore, overall mean ages of POTS and/or healthy group could not be calculated because data could not be accessed in other studies(28, 32), because median was utilized than mean (25, 26).
Heart rate
Figure 2 illustrates eighteen of twenty eligible studies that compared heart rate outcome between POTS versus Healthy participants. Overall mean difference between the two groups was 19.88 (15.24, 24.52) signifying higher heart rate in POTS group. This difference reached statistical significance (P-value<0.05). Random-effect model was used since heterogeneity, I2, was 99% (i.e. I2 >50%).
RR- interval
Figure 3a illustrates three of twenty eligible studies that compared time domain measure outcome between POTS versus Healthy participants in terms of mean RR intervals. Overall mean difference between the two groups was -162.89 (-172.65, -153.12) signifying lower Heart rate variability in terms of RR-interval in POTS group. This difference reached statistical significance (P-value<0.05). Fixed-effect model was used since heterogeneity, I2, was 0% (i.e. I2 <50%).
Root of mean of squares of successive R-R interval differences (rMSSD)
Figure 3b illustrates three of twenty eligible studies that compared time domain measure outcome between POTS versus Healthy participants in terms of rMSSD. Overall mean difference between the two groups was -15.16 (-18.28, -12.03) signifying lower Heart rate variability in terms of rMSSD in POTS group. The difference reached statistical significance (P-value<0.05). Fixed-effect model was used since heterogeneity, I2, was 2% (i.e. I2 <50%).
Low frequency power (LF)
Figure 4a illustrates five of twenty eligible studies that compared frequency domain measure outcome between POTS versus Healthy participants in terms of LF. Overall mean difference between the two groups was -80.89 (-211.37, 49.58) milliseconds2 signifying lower heart rate variability in terms of LF in POTS group. The difference, however did not reach statistical significance (P-value>0.05). Random-effect model was used since heterogeneity, I2, was 96% (i.e. I2 >50%).
High frequency power (HF)
Figure 4b illustrates seven of twenty eligible studies that compared frequency domain measure outcome between POTS versus Healthy participants in terms of HF. Overall mean difference between the two groups was -113.20 (-275.52, 49.13) milliseconds2 signifying lower heart rate variability in terms of HF in POTS group. The difference, did not reach statistical significance (P-value>0.05). Random-effect model was used since heterogeneity, I2, was 84% (i.e. I2 >50%).
Low frequency power /High frequency power ratio (LF/HF- ratio)
Figure 4c illustrates five of twenty eligible studies that compared frequency domain measure outcome between POTS versus Healthy participants in terms of LF/HF- ratio. Overall mean difference between the two groups was 0.29 (-0.25, 0.83) signifying higher heart rate variability in terms of LF/HF- ratio in POTS group. The difference, did not reach statistical significance (P-value>0.05). Random-effect model was used since heterogeneity, I2, was 95% (i.e. I2 >50%).
Low frequency power-normalized unit
Figure 4d illustrates four of twenty eligible studies that compared frequency domain measure outcome between POTS versus Healthy participants in terms of LF (n.u). Overall mean difference between the two groups was 0.05 (-0.04, 0.13) signifying higher heart rate variability in terms of LF (n.u.) in POTS group. The difference, however did not reach statistical significance (P-value>0.05). Random-effect model was used since heterogeneity, I2, was 96% (i.e. I2 >50%).
High frequency power-normalized unit
Figure 4e illustrates four of twenty eligible studies that compared frequency domain measure outcome between POTS versus Healthy participants in terms of HF (n.u). Overall mean difference between the two groups was -0.03 (-0.11, 0.04) signifying lower Heart variability in terms of HF (n.u.) in POTS group. The difference, however did not reach statistical significance (P-value>0.05). Random-effect model was used since heterogeneity, I2, was 91% (i.e. I2 >50%).
Sensitivity analysis
Eliminating three studies; one study for utilizing 24 hours parameters recording (32); another study for including a period of parameters measurements during sleeping (22); and one study (23), in which data were collected by estimates and extrapolations from a graphical figure, none of the outcome results changed statistical significance. The newly, obtained results were as follows; mean difference rMSSD =-15.41(-18.63,-12.2), p-Value<0.00001, I2=38%; mean difference HF =-156(-344.05,31.17), p-Value=0.1, I2=84%; mean difference LF/HF =0.39(-0.24,1.03), p-Value=0.23, I2=96%; and lastly, mean HR =20.98(16.27,25.69), p-Value<0.00001, I2=99%.
Publication bias
Figure 5 illustrates a funnel-plot to illustrate publication biases for studies included in comparing heart rate between POTS versus healthy groups. Medium sample sized studies at the middle of the funnel-plot were more symmetrically distributed as compared to large sample sized studies at the top. This suggest heterogeneity of the study estimates as well as likely publication bias favoring studies with medium sample sizes than large sample sizes.