Abstract

BACKGROUND

Although blood pressure (BP) is regulated by the autonomic nervous system, it is not fully understood how autonomic activity affects BP at home in the general population.

METHODS

Subjects were enrolled from 2009 to 2012 and included 1,888 men and women aged 30–79 years. We measured casual BP in the morning during health checkups and asked participants to monitor BP at home twice in the morning and evening for 1 week. The mean of the two measurements of mean arterial pressure (MAP) was calculated. Five-minute recordings of the pulse wave from a fingertip sensor were used to determine the following indices of heart rate variability (HRV): standard deviation of normal-to-normal RR intervals (SDNN), root mean square of successive differences in RR intervals (RMSSD), high frequency (HF) power, low frequency (LF) power, and LF/HF.

RESULTS

Sex- and age-adjusted means of casual MAP, and morning and evening MAP at home were significantly different among quartiles of SDNN, RMSSD, and HF. When further adjusted for smoking, alcohol drinking, medication for hypertension, diabetes, sleeping hours, snoring, and mental health status, the associations were somewhat attenuated. Inverse relationships were found between the means of morning home MAP, and RMSSD (P = 0.02) and HF (P = 0.051) after adjustment for confounders. The association between MAP and RMSSD, or MAP and HF was evident in individuals <65 years old.

CONCLUSION

Low HF and RMSSD, which reflect impaired parasympathetic nervous system activity, were associated with increased home MAP in the morning rather than in the evening.

The autonomic nervous system plays an important role in the regulation of blood pressure (BP),1 and autonomic dysfunction can cause hypertension that is associated with reduced baroreflex sensitivity.2–4 Recently, the joint scientific statement on hypertension by several professional societies strongly recommended the use of home BP monitoring both to prevent cardiovascular outcomes and to control hypertension.5,6 Given the pathological mechanisms that can raise BP, sympathovagal imbalance may have an important effect on home BP; however, there are few epidemiological studies on the association between sympathovagal imbalance and home BP in the general population.

Heart rate variability (HRV) is regulated predominantly by cardiac vagal tone, and lower HRV is associated with an increased risk of diabetes, hypertension, and cardiovascular disease (CVD) in epidemiologic studies.7–10 We have reported that reduced HRV was associated with insulin resistance,11 metabolic syndrome,12 and elevated C-reactive protein concentration.13 In addition, we found that impaired parasympathetic function was associated with diastolic BP but not systolic BP at health checkups.14 This was explained to some extent by baroreflex sensitivity, which was inversely related to mean arterial pressure (MAP) and diastolic BP rather than systolic BP.2

Our study hypothesis was that autonomic dysfunction based on standardized indices of HRV15 obtained during health checkups may be associated with an increased MAP at home. HRV assessment may also help to manage BP at home.

METHODS

Study subjects

The Toon Health Study enrolled 2,032 men and women, 30–79 years of age between 2009 and 2012. The study included subjects who had HRV assessment (n = 2,030) and did not have atrial fibrillation (n = 15) on an electrocardiogram (ECG). Among 2,019 participants who had BP recorded at home, we included data from 1,905 who had morning and evening BP recorded for more than 3 days. Based on these criteria, 1,888 individuals were included in the analysis.

Written informed consent was obtained from all participants. The study protocol was approved by the Human Ethics Review Committees of Ehime University Graduate School of Medicine (no: 1705011).

Measurements

Overnight fasting blood samples were drawn from the antecubital vein into vacuum tubes containing a serum separator gel (for glucose and blood chemistry). The serum tube was centrifuged immediately at 3,000 rpm for 15 minutes, and the separated serum was sent to the laboratory for analysis.

BP was measured twice with individuals in the sitting position after a rest of at least 5 minutes using an automatic sphygmomanometer (BP-103iII; OMRON Colin Co., Ltd., Tokyo, Japan). The mean of the two measurements was used for analysis. All measurements were done from 8 to 11 am. Participants were asked to monitor their BP at home twice in the morning (within 1 hour after waking) and twice in the evening (before going to bed) using an automated device (UA-767; A&D Co., Ltd., Tokyo, Japan) that stored all BP measurements in the device memory. Subjects were required to measure their home BP for 1 week. Home BP was defined as the mean of the first two measurements of systolic and diastolic BP during each time period. Hypertension was defined as a systolic BP ≥140 mm Hg, and/or a diastolic BP ≥90 mm Hg, and/or the current use of any antihypertensive medication. Type 2 diabetes was defined as a fasting plasma glucose ≥7.0 mmol/l, a 2-hour-postload glucose ≥11.1 mmol/l, or the current use of antihyperglycemic agents. Body mass index (BMI) was calculated as weight (kg) divided by height (m) squared.

A self-administered questionnaire was used to assess medical history (presence of hypertension, dyslipidemia, or diabetes), smoking habits, and alcohol consumption. The amount of alcohol consumed each week was evaluated by measuring the weekly frequency of drinking and the type of alcoholic beverage consumed (beer, sake, whiskey, shochu, or wine). A regular alcohol drinker was defined as an individual with alcohol consumption ≥1 g/wk. Snoring was defined as a “Yes” response among three possible answers (“Yes,” “No,” or “Unknown”) to a question about snoring (“Do you snore?”). Self-reported snoring has been shown to have a 94% sensitivity and 58% specificity compared with recordings of snoring sounds as a gold standard.16 Furthermore, self-reported snoring has also been validated using pulse oximetry and cardiovascular disease outcomes.17 Physical activity levels were assessed using a validated questionnaire that consisted of 14 questions on occupation, locomotion, housework, sleep time, and leisure time physical activities.18 Responses for each physical activity category were converted to metabolic equivalents (METs), according to the Compendium by Ainsworth et al.,19 and expressed as METs·h/day. We identified mental health status using the Short Form-8 (SF-8) mental component score (MCS).20

Assessment of autonomic function

HRV was measured with a fingertip pulse wave sensor (TAS9; YKC, Co., Ltd., Tokyo, Japan). The 5-minute pulse rate was recorded, and then the standard deviation of normal-to-normal RR intervals (SDNN) and root mean square of successive differences in RR intervals (RMSSD) were determined. Power spectral analysis of the pulse recordings was also used to obtain frequency-domain measures of HRV. The power spectrum was decomposed into its frequency components and quantified in terms of the relative intensity (power) of each component. The power spectrum was divided into frequency bands, and we determined the high frequency band (HF) (0.15–0.40 Hz) and low frequency band (LF) (0.04–0.15 Hz). The HF and LF power and the LF/HF ratio were used for further analysis.

We performed 24-hour Holter ECG monitoring (PMP400; Pacific Medico, Co., Ltd., Tokyo, Japan) in 140 volunteers among 557 study participants in 2015. To validate the 5-minute HRV measurement using the fingertip device, power spectral analysis of RR intervals from the ECG was performed every 5 minutes over 24 hours in 139 individuals who had complete records. For each 5-minute interval, we calculated LF, HF, and their ratio (LF/HF) in the same frequency bands.

To assess the reliability of the HRV parameters, each parameter was measured twice in the same individual at an interval of 2 months (n = 37). Pearson’s correlation coefficients between the two measurements within subjects were 0.67 for SDNN, 0.51 for RMSSD, 0.63 for LF, 0.54 for HF, and 0.21 for LF/HF. Coefficients of variation were 0.15, 0.20, 0.29, 0.27, and 0.20, respectively, which were similar to the values reported in a previously published validation study.21

Statistical analysis

Because of skewed distributions, SDNN, RMSSD, LF, and HF were log-transformed before analysis. The LF/HF ratio was calculated using log-transformed values of LF and HF. Sex- and age-adjusted means were calculated using analysis of covariance. The trend test was performed using log-transformed continuous values of SDNN, RMSSD, LF, HF, and LF/HF. Multivariable models were further adjusted for smoking, alcohol drinking, the use of hypertensive agents, diabetes, sleep duration, snoring, physical activity, and SF-8 MCS. A test for interaction of the stratified variables was done using the interaction terms that consisted of log-transformed values of SDNN, RMSSD, LF, HF, and LF/HF ratio as continuous variables adjusted for sex and age. Statistical significance was assumed at P <0.05. All statistical analyses were performed using SAS software, version 9.4 (SAS Institute, Inc., Cary, NC).

RESULTS

Table 1 shows the characteristics of all subjects included in the study. The mean values of systolic BP, diastolic BP, and MAP were lower when they were measured at home than in the health checkups.

Table 1.

Subject characteristics

VariablesValues
(n = 1,888)
Age, years old58.1 ± 12.5
Men, %34.3
BMI, kg/m223.1 ± 3.2
Waist circumference, cm83.2 ± 9.2
Casual systolic BP, mm Hg126.3 ± 20.0
Casual diastolic BP, mm Hg76.1 ± 11.8
Casual MAP, mm Hg92.8 ± 13.9
Systolic BP at home
 Morning, mm Hg120.0 ± 17.2
 Evening, mm Hg116.3 ± 15.7
Diastolic BP at home
 Morning, mm Hg75.3 ± 10.7
 Evening, mm Hg70.6 ± 10.0
MAP at home
 Morning, mm Hg90.2 ± 12.1
 Evening, mm Hg85.8 ± 11.2
Hypertension, %37.3
Diabetes, %10.2
Sleep, hours6.4 ± 1.1
Snoring, %39.5
Current smoker, %8.3
Regular alcohol drinker, %51.1
Physical activity, METs·h/day35.6 ± 4.5
VariablesValues
(n = 1,888)
Age, years old58.1 ± 12.5
Men, %34.3
BMI, kg/m223.1 ± 3.2
Waist circumference, cm83.2 ± 9.2
Casual systolic BP, mm Hg126.3 ± 20.0
Casual diastolic BP, mm Hg76.1 ± 11.8
Casual MAP, mm Hg92.8 ± 13.9
Systolic BP at home
 Morning, mm Hg120.0 ± 17.2
 Evening, mm Hg116.3 ± 15.7
Diastolic BP at home
 Morning, mm Hg75.3 ± 10.7
 Evening, mm Hg70.6 ± 10.0
MAP at home
 Morning, mm Hg90.2 ± 12.1
 Evening, mm Hg85.8 ± 11.2
Hypertension, %37.3
Diabetes, %10.2
Sleep, hours6.4 ± 1.1
Snoring, %39.5
Current smoker, %8.3
Regular alcohol drinker, %51.1
Physical activity, METs·h/day35.6 ± 4.5

Values are expressed as mean ± standard deviation or percentage. Abbreviations: BMI, body mass index; BP, blood pressure; MAP, mean arterial pressure; METs, metabolic equivalents.

Table 1.

Subject characteristics

VariablesValues
(n = 1,888)
Age, years old58.1 ± 12.5
Men, %34.3
BMI, kg/m223.1 ± 3.2
Waist circumference, cm83.2 ± 9.2
Casual systolic BP, mm Hg126.3 ± 20.0
Casual diastolic BP, mm Hg76.1 ± 11.8
Casual MAP, mm Hg92.8 ± 13.9
Systolic BP at home
 Morning, mm Hg120.0 ± 17.2
 Evening, mm Hg116.3 ± 15.7
Diastolic BP at home
 Morning, mm Hg75.3 ± 10.7
 Evening, mm Hg70.6 ± 10.0
MAP at home
 Morning, mm Hg90.2 ± 12.1
 Evening, mm Hg85.8 ± 11.2
Hypertension, %37.3
Diabetes, %10.2
Sleep, hours6.4 ± 1.1
Snoring, %39.5
Current smoker, %8.3
Regular alcohol drinker, %51.1
Physical activity, METs·h/day35.6 ± 4.5
VariablesValues
(n = 1,888)
Age, years old58.1 ± 12.5
Men, %34.3
BMI, kg/m223.1 ± 3.2
Waist circumference, cm83.2 ± 9.2
Casual systolic BP, mm Hg126.3 ± 20.0
Casual diastolic BP, mm Hg76.1 ± 11.8
Casual MAP, mm Hg92.8 ± 13.9
Systolic BP at home
 Morning, mm Hg120.0 ± 17.2
 Evening, mm Hg116.3 ± 15.7
Diastolic BP at home
 Morning, mm Hg75.3 ± 10.7
 Evening, mm Hg70.6 ± 10.0
MAP at home
 Morning, mm Hg90.2 ± 12.1
 Evening, mm Hg85.8 ± 11.2
Hypertension, %37.3
Diabetes, %10.2
Sleep, hours6.4 ± 1.1
Snoring, %39.5
Current smoker, %8.3
Regular alcohol drinker, %51.1
Physical activity, METs·h/day35.6 ± 4.5

Values are expressed as mean ± standard deviation or percentage. Abbreviations: BMI, body mass index; BP, blood pressure; MAP, mean arterial pressure; METs, metabolic equivalents.

Table 2 shows the correlation coefficients between HRV parameters estimated from 5-minute fingertip pulse recordings and the same HRV parameters obtained from 24-hour Holter ECG recordings. The mean HRV parameters from successive 5-minute ECG recordings were determined during five different time periods: morning, awake, evening, asleep, and all 24 hours. Heart rate, LF, HF, and LF/HF measured from pulse recordings were moderately correlated with the means of these same parameters measured from the ECG during each time period. The means of heart rate, ln LF, ln HF, and ln LF/HF in the 24-hour Holter ECG are shown by time period in Supplementary Table S1. Ln HF was higher and ln LF/HF was lower in the evening than in the morning.

Table 2.

Pearson’s correlation coefficients between HRV parameters from 5-minute measurements during a health checkup and the means of each time period obtained from 24-hour Holter electrocardiographic recordings (n = 139)

HRV parametersTime period
Morning
(7–8 am)
Awake
(9 am–9 pm)
Evening
(10–11 pm)
Asleep
(12–6 am)
All 24 hours
rPrPrPrPrP
Heart rate0.56<0.0010.74<0.0010.54<0.0010.53<0.0010.73<0.001
LF0.34<0.0010.50<0.0010.47<0.0010.48<0.0010.53<0.001
HF0.43<0.0010.57<0.0010.51<0.0010.53<0.0010.59<0.001
LF/HF0.26<0.010.35<0.0010.35<0.0010.44<0.0010.42<0.001
HRV parametersTime period
Morning
(7–8 am)
Awake
(9 am–9 pm)
Evening
(10–11 pm)
Asleep
(12–6 am)
All 24 hours
rPrPrPrPrP
Heart rate0.56<0.0010.74<0.0010.54<0.0010.53<0.0010.73<0.001
LF0.34<0.0010.50<0.0010.47<0.0010.48<0.0010.53<0.001
HF0.43<0.0010.57<0.0010.51<0.0010.53<0.0010.59<0.001
LF/HF0.26<0.010.35<0.0010.35<0.0010.44<0.0010.42<0.001

Abbreviations: HRV, heart rate variability; HF, high frequency; LF, low frequency.

Table 2.

Pearson’s correlation coefficients between HRV parameters from 5-minute measurements during a health checkup and the means of each time period obtained from 24-hour Holter electrocardiographic recordings (n = 139)

HRV parametersTime period
Morning
(7–8 am)
Awake
(9 am–9 pm)
Evening
(10–11 pm)
Asleep
(12–6 am)
All 24 hours
rPrPrPrPrP
Heart rate0.56<0.0010.74<0.0010.54<0.0010.53<0.0010.73<0.001
LF0.34<0.0010.50<0.0010.47<0.0010.48<0.0010.53<0.001
HF0.43<0.0010.57<0.0010.51<0.0010.53<0.0010.59<0.001
LF/HF0.26<0.010.35<0.0010.35<0.0010.44<0.0010.42<0.001
HRV parametersTime period
Morning
(7–8 am)
Awake
(9 am–9 pm)
Evening
(10–11 pm)
Asleep
(12–6 am)
All 24 hours
rPrPrPrPrP
Heart rate0.56<0.0010.74<0.0010.54<0.0010.53<0.0010.73<0.001
LF0.34<0.0010.50<0.0010.47<0.0010.48<0.0010.53<0.001
HF0.43<0.0010.57<0.0010.51<0.0010.53<0.0010.59<0.001
LF/HF0.26<0.010.35<0.0010.35<0.0010.44<0.0010.42<0.001

Abbreviations: HRV, heart rate variability; HF, high frequency; LF, low frequency.

Sex- and age-adjusted, and multivariable-adjusted means of casual and home MAP by quartile of HRV parameters are shown in Table 3. Sex- and age-adjusted means of casual, and home-measured morning and evening MAP were significantly different among quartiles of SDNN, RMSSD, and HF. When further adjusted for smoking, alcohol drinking, medication for hypertension, diabetes, sleeping hours, snoring, and SF-8 MCS, the associations were attenuated. Inverse relationships were found between the mean of morning MAP at home and RMSSD (P = 0.02) and HF (P = 0.051) after adjustment for confounders. The means of casual, and home systolic and diastolic BPs adjusted for covariates are shown by quartile of the HRV parameters in Supplementary Tables S2 and S3. Morning diastolic BP at home was strongly associated with both RMSSD and HF as well as MAP, and these associations were independent of covariates.

Table 3.

Multivariable-adjusted means of casual and home blood pressures grouped by HRV parameter quartiles

HRV parametersQuartileMAP, mm Hg
Sex- and age-adjusted meansMultivariable-adjusted means
CasualHomeCasualHome
MorningEveningMorningEvening
SDNNQ194.992.386.995.692.386.9
Q293.991.887.295.192.387.8
Q393.190.786.294.991.887.3
Q492.690.686.194.391.487.2
P for difference0.0480.060.300.450.500.52
P for trend<0.01<0.010.090.100.070.77
RMSSDQ195.893.088.196.493.088.0
Q293.991.386.495.291.987.1
Q392.490.786.293.891.386.8
Q492.490.585.894.391.587.0
P for difference<0.001<0.01<0.01<0.010.070.26
P for trend<0.001<0.01<0.01<0.010.020.09
LFQ194.292.387.095.192.587.1
Q293.691.486.894.891.887.2
Q393.691.286.795.492.287.9
Q493.190.585.994.691.386.8
P for difference0.620.120.380.770.340.39
P for trend0.14<0.010.080.500.070.51
HFQ195.593.088.195.992.787.7
Q294.091.286.395.591.987.1
Q392.490.986.294.291.987.3
Q492.590.485.894.291.386.9
P for difference<0.01<0.01<0.010.070.250.57
P for trend<0.01<0.01<0.010.040.0510.21
LF/HFQ193.191.186.294.691.887.1
Q293.591.686.494.992.287.1
Q392.990.986.694.491.587.3
Q494.991.887.196.092.187.5
P for difference0.060.550.620.140.750.92
P for trend0.010.470.120.040.960.43
HRV parametersQuartileMAP, mm Hg
Sex- and age-adjusted meansMultivariable-adjusted means
CasualHomeCasualHome
MorningEveningMorningEvening
SDNNQ194.992.386.995.692.386.9
Q293.991.887.295.192.387.8
Q393.190.786.294.991.887.3
Q492.690.686.194.391.487.2
P for difference0.0480.060.300.450.500.52
P for trend<0.01<0.010.090.100.070.77
RMSSDQ195.893.088.196.493.088.0
Q293.991.386.495.291.987.1
Q392.490.786.293.891.386.8
Q492.490.585.894.391.587.0
P for difference<0.001<0.01<0.01<0.010.070.26
P for trend<0.001<0.01<0.01<0.010.020.09
LFQ194.292.387.095.192.587.1
Q293.691.486.894.891.887.2
Q393.691.286.795.492.287.9
Q493.190.585.994.691.386.8
P for difference0.620.120.380.770.340.39
P for trend0.14<0.010.080.500.070.51
HFQ195.593.088.195.992.787.7
Q294.091.286.395.591.987.1
Q392.490.986.294.291.987.3
Q492.590.485.894.291.386.9
P for difference<0.01<0.01<0.010.070.250.57
P for trend<0.01<0.01<0.010.040.0510.21
LF/HFQ193.191.186.294.691.887.1
Q293.591.686.494.992.287.1
Q392.990.986.694.491.587.3
Q494.991.887.196.092.187.5
P for difference0.060.550.620.140.750.92
P for trend0.010.470.120.040.960.43

MAP was adjusted for sex, age, smoking, alcohol drinking, the use of hypertensive agents, diabetes, sleep duration, snoring, physical activity, and SF-8 mental component score by ANCOVA. Abbreviations: MAP, mean arterial pressure; HRV, heart rate variability; SDNN, standard deviation of normal-to-normal RR intervals; RMSSD, root mean square of successive RR interval differences; LF, low frequency; HF, high frequency.

Table 3.

Multivariable-adjusted means of casual and home blood pressures grouped by HRV parameter quartiles

HRV parametersQuartileMAP, mm Hg
Sex- and age-adjusted meansMultivariable-adjusted means
CasualHomeCasualHome
MorningEveningMorningEvening
SDNNQ194.992.386.995.692.386.9
Q293.991.887.295.192.387.8
Q393.190.786.294.991.887.3
Q492.690.686.194.391.487.2
P for difference0.0480.060.300.450.500.52
P for trend<0.01<0.010.090.100.070.77
RMSSDQ195.893.088.196.493.088.0
Q293.991.386.495.291.987.1
Q392.490.786.293.891.386.8
Q492.490.585.894.391.587.0
P for difference<0.001<0.01<0.01<0.010.070.26
P for trend<0.001<0.01<0.01<0.010.020.09
LFQ194.292.387.095.192.587.1
Q293.691.486.894.891.887.2
Q393.691.286.795.492.287.9
Q493.190.585.994.691.386.8
P for difference0.620.120.380.770.340.39
P for trend0.14<0.010.080.500.070.51
HFQ195.593.088.195.992.787.7
Q294.091.286.395.591.987.1
Q392.490.986.294.291.987.3
Q492.590.485.894.291.386.9
P for difference<0.01<0.01<0.010.070.250.57
P for trend<0.01<0.01<0.010.040.0510.21
LF/HFQ193.191.186.294.691.887.1
Q293.591.686.494.992.287.1
Q392.990.986.694.491.587.3
Q494.991.887.196.092.187.5
P for difference0.060.550.620.140.750.92
P for trend0.010.470.120.040.960.43
HRV parametersQuartileMAP, mm Hg
Sex- and age-adjusted meansMultivariable-adjusted means
CasualHomeCasualHome
MorningEveningMorningEvening
SDNNQ194.992.386.995.692.386.9
Q293.991.887.295.192.387.8
Q393.190.786.294.991.887.3
Q492.690.686.194.391.487.2
P for difference0.0480.060.300.450.500.52
P for trend<0.01<0.010.090.100.070.77
RMSSDQ195.893.088.196.493.088.0
Q293.991.386.495.291.987.1
Q392.490.786.293.891.386.8
Q492.490.585.894.391.587.0
P for difference<0.001<0.01<0.01<0.010.070.26
P for trend<0.001<0.01<0.01<0.010.020.09
LFQ194.292.387.095.192.587.1
Q293.691.486.894.891.887.2
Q393.691.286.795.492.287.9
Q493.190.585.994.691.386.8
P for difference0.620.120.380.770.340.39
P for trend0.14<0.010.080.500.070.51
HFQ195.593.088.195.992.787.7
Q294.091.286.395.591.987.1
Q392.490.986.294.291.987.3
Q492.590.485.894.291.386.9
P for difference<0.01<0.01<0.010.070.250.57
P for trend<0.01<0.01<0.010.040.0510.21
LF/HFQ193.191.186.294.691.887.1
Q293.591.686.494.992.287.1
Q392.990.986.694.491.587.3
Q494.991.887.196.092.187.5
P for difference0.060.550.620.140.750.92
P for trend0.010.470.120.040.960.43

MAP was adjusted for sex, age, smoking, alcohol drinking, the use of hypertensive agents, diabetes, sleep duration, snoring, physical activity, and SF-8 mental component score by ANCOVA. Abbreviations: MAP, mean arterial pressure; HRV, heart rate variability; SDNN, standard deviation of normal-to-normal RR intervals; RMSSD, root mean square of successive RR interval differences; LF, low frequency; HF, high frequency.

Regression analysis was performed using casual, home morning, or home evening MAP as dependent variables, and log-transformed values of RMSSD or HF as independent variables (Table 4). Partial regression coefficients were determined after subjects were stratified into two groups for each of the following six variables: sex, age, overweight, alcohol drinking, sleeping duration, and snoring. For the association of morning MAP at home with RMSSD or HF, there were no significant differences in the partial regression coefficients between the two groups (P > 0.05 for interaction) for any of the six variables included in the analyses. However, for the association of evening MAP at home with RMSSD or HF, there was a significant difference in the partial regression coefficients between the two age groups (P < 0.05 for interaction). A significant association between MAP and RMSSD or MAP and HF was evident in individuals <65 years old.

Table 4.

Partial regression coefficients for the association of MAP (casual and home) with RMSSD or HF in subjects stratified by selected variables

VariablesStratificationIndependent: ln RMSSDIndependent: ln HF
Dependent: MAP, mm HgDependent: MAP, mm Hg
CasualHomeCasualHome
MorningEveningMorningEvening
β (SE)β (SE)β (SE)β (SE)β (SE)β (SE)
SexMen−1.48 (0.77)−1.63 (0.70)*−0.86 (0.66)−0.54 (0.39)−0.80 (0.35)*−0.42 (0.33)
Women−1.95 (0.59)**−1.19 (0.49)*−1.28 (0.48)**−0.80 (0.30)**−0.50 (0.25)*−0.55 (0.25)*
P for interaction0.360.940.300.220.970.22
Age group<65 years−2.63 (0.65)**−1.74 (0.57)**−1.72 (0.55)**−1.04 (0.33)**−0.85 (0.29)**−0.79 (0.28)**
≥65 years−0.31 (0.68)−0.40 (0.57)−0.09 (0.55)−0.09 (0.34)−0.11 (0.28)−0.01 (0.27)
P for interaction0.020.120.0470.040.060.04
OverweightYes0.01 (0.90)−0.61 (0.77)−0.41 (0.77)0.06 (0.46)−0.41 (0.39)−0.35 (0.39)
No−1.91 (0.53)**−1.10 (0.46)*−0.84 (0.43)−0.72 (0.27)**−0.42 (0.23)−0.27 (0.22)
P for interaction0.040.450.480.040.590.76
Alcohol drinkingYes−2.52 (0.65)**−1.93 (0.58)**−1.57 (0.56)**−1.19 (0.33)**−1.00 (0.30)**−0.82 (0.29)**
No−1.21 (0.68)−0.89 (0.57)−0.85 (0.55)−0.35 (0.34)−0.31 (0.29)−0.30 (0.28)
P for interaction0.190.220.450.100.100.29
SnoringYes−1.95 (0.72)**−1.30 (0.65)*−1.14 (0.63)−0.70 (0.37)−0.61 (0.34)−0.46 (0.33)
No−1.69 (0.61)**−1.41 (0.51)**−1.19 (0.50)*−0.72 (0.31)*−0.63 (0.25)*−0.57 (0.25)*
P for interaction0.870.690.790.850.690.59
Sleep duration<7h−2.16 (0.67)**−1.35 (0.59)*−0.94 (0.56)−1.12 (0.34)**−0.83 (0.30)**−0.61 (0.28)*
≥7h−1.57 (0.66)*−1.45 (0.56)**−1.41 (0.56)*−0.42 (0.33)−0.48 (0.28)−0.47 (0.28)
P for interaction0.550.900.540.170.400.75
VariablesStratificationIndependent: ln RMSSDIndependent: ln HF
Dependent: MAP, mm HgDependent: MAP, mm Hg
CasualHomeCasualHome
MorningEveningMorningEvening
β (SE)β (SE)β (SE)β (SE)β (SE)β (SE)
SexMen−1.48 (0.77)−1.63 (0.70)*−0.86 (0.66)−0.54 (0.39)−0.80 (0.35)*−0.42 (0.33)
Women−1.95 (0.59)**−1.19 (0.49)*−1.28 (0.48)**−0.80 (0.30)**−0.50 (0.25)*−0.55 (0.25)*
P for interaction0.360.940.300.220.970.22
Age group<65 years−2.63 (0.65)**−1.74 (0.57)**−1.72 (0.55)**−1.04 (0.33)**−0.85 (0.29)**−0.79 (0.28)**
≥65 years−0.31 (0.68)−0.40 (0.57)−0.09 (0.55)−0.09 (0.34)−0.11 (0.28)−0.01 (0.27)
P for interaction0.020.120.0470.040.060.04
OverweightYes0.01 (0.90)−0.61 (0.77)−0.41 (0.77)0.06 (0.46)−0.41 (0.39)−0.35 (0.39)
No−1.91 (0.53)**−1.10 (0.46)*−0.84 (0.43)−0.72 (0.27)**−0.42 (0.23)−0.27 (0.22)
P for interaction0.040.450.480.040.590.76
Alcohol drinkingYes−2.52 (0.65)**−1.93 (0.58)**−1.57 (0.56)**−1.19 (0.33)**−1.00 (0.30)**−0.82 (0.29)**
No−1.21 (0.68)−0.89 (0.57)−0.85 (0.55)−0.35 (0.34)−0.31 (0.29)−0.30 (0.28)
P for interaction0.190.220.450.100.100.29
SnoringYes−1.95 (0.72)**−1.30 (0.65)*−1.14 (0.63)−0.70 (0.37)−0.61 (0.34)−0.46 (0.33)
No−1.69 (0.61)**−1.41 (0.51)**−1.19 (0.50)*−0.72 (0.31)*−0.63 (0.25)*−0.57 (0.25)*
P for interaction0.870.690.790.850.690.59
Sleep duration<7h−2.16 (0.67)**−1.35 (0.59)*−0.94 (0.56)−1.12 (0.34)**−0.83 (0.30)**−0.61 (0.28)*
≥7h−1.57 (0.66)*−1.45 (0.56)**−1.41 (0.56)*−0.42 (0.33)−0.48 (0.28)−0.47 (0.28)
P for interaction0.550.900.540.170.400.75

Partial regression coefficients were calculated after adjusting MAP for sex and age. Abbreviations: MAP, mean arterial pressure; SE, standard error of mean; RMSSD, root mean square of successive RR interval differences; HF, high frequency.

*P < 0.05; **P < 0.01.

Table 4.

Partial regression coefficients for the association of MAP (casual and home) with RMSSD or HF in subjects stratified by selected variables

VariablesStratificationIndependent: ln RMSSDIndependent: ln HF
Dependent: MAP, mm HgDependent: MAP, mm Hg
CasualHomeCasualHome
MorningEveningMorningEvening
β (SE)β (SE)β (SE)β (SE)β (SE)β (SE)
SexMen−1.48 (0.77)−1.63 (0.70)*−0.86 (0.66)−0.54 (0.39)−0.80 (0.35)*−0.42 (0.33)
Women−1.95 (0.59)**−1.19 (0.49)*−1.28 (0.48)**−0.80 (0.30)**−0.50 (0.25)*−0.55 (0.25)*
P for interaction0.360.940.300.220.970.22
Age group<65 years−2.63 (0.65)**−1.74 (0.57)**−1.72 (0.55)**−1.04 (0.33)**−0.85 (0.29)**−0.79 (0.28)**
≥65 years−0.31 (0.68)−0.40 (0.57)−0.09 (0.55)−0.09 (0.34)−0.11 (0.28)−0.01 (0.27)
P for interaction0.020.120.0470.040.060.04
OverweightYes0.01 (0.90)−0.61 (0.77)−0.41 (0.77)0.06 (0.46)−0.41 (0.39)−0.35 (0.39)
No−1.91 (0.53)**−1.10 (0.46)*−0.84 (0.43)−0.72 (0.27)**−0.42 (0.23)−0.27 (0.22)
P for interaction0.040.450.480.040.590.76
Alcohol drinkingYes−2.52 (0.65)**−1.93 (0.58)**−1.57 (0.56)**−1.19 (0.33)**−1.00 (0.30)**−0.82 (0.29)**
No−1.21 (0.68)−0.89 (0.57)−0.85 (0.55)−0.35 (0.34)−0.31 (0.29)−0.30 (0.28)
P for interaction0.190.220.450.100.100.29
SnoringYes−1.95 (0.72)**−1.30 (0.65)*−1.14 (0.63)−0.70 (0.37)−0.61 (0.34)−0.46 (0.33)
No−1.69 (0.61)**−1.41 (0.51)**−1.19 (0.50)*−0.72 (0.31)*−0.63 (0.25)*−0.57 (0.25)*
P for interaction0.870.690.790.850.690.59
Sleep duration<7h−2.16 (0.67)**−1.35 (0.59)*−0.94 (0.56)−1.12 (0.34)**−0.83 (0.30)**−0.61 (0.28)*
≥7h−1.57 (0.66)*−1.45 (0.56)**−1.41 (0.56)*−0.42 (0.33)−0.48 (0.28)−0.47 (0.28)
P for interaction0.550.900.540.170.400.75
VariablesStratificationIndependent: ln RMSSDIndependent: ln HF
Dependent: MAP, mm HgDependent: MAP, mm Hg
CasualHomeCasualHome
MorningEveningMorningEvening
β (SE)β (SE)β (SE)β (SE)β (SE)β (SE)
SexMen−1.48 (0.77)−1.63 (0.70)*−0.86 (0.66)−0.54 (0.39)−0.80 (0.35)*−0.42 (0.33)
Women−1.95 (0.59)**−1.19 (0.49)*−1.28 (0.48)**−0.80 (0.30)**−0.50 (0.25)*−0.55 (0.25)*
P for interaction0.360.940.300.220.970.22
Age group<65 years−2.63 (0.65)**−1.74 (0.57)**−1.72 (0.55)**−1.04 (0.33)**−0.85 (0.29)**−0.79 (0.28)**
≥65 years−0.31 (0.68)−0.40 (0.57)−0.09 (0.55)−0.09 (0.34)−0.11 (0.28)−0.01 (0.27)
P for interaction0.020.120.0470.040.060.04
OverweightYes0.01 (0.90)−0.61 (0.77)−0.41 (0.77)0.06 (0.46)−0.41 (0.39)−0.35 (0.39)
No−1.91 (0.53)**−1.10 (0.46)*−0.84 (0.43)−0.72 (0.27)**−0.42 (0.23)−0.27 (0.22)
P for interaction0.040.450.480.040.590.76
Alcohol drinkingYes−2.52 (0.65)**−1.93 (0.58)**−1.57 (0.56)**−1.19 (0.33)**−1.00 (0.30)**−0.82 (0.29)**
No−1.21 (0.68)−0.89 (0.57)−0.85 (0.55)−0.35 (0.34)−0.31 (0.29)−0.30 (0.28)
P for interaction0.190.220.450.100.100.29
SnoringYes−1.95 (0.72)**−1.30 (0.65)*−1.14 (0.63)−0.70 (0.37)−0.61 (0.34)−0.46 (0.33)
No−1.69 (0.61)**−1.41 (0.51)**−1.19 (0.50)*−0.72 (0.31)*−0.63 (0.25)*−0.57 (0.25)*
P for interaction0.870.690.790.850.690.59
Sleep duration<7h−2.16 (0.67)**−1.35 (0.59)*−0.94 (0.56)−1.12 (0.34)**−0.83 (0.30)**−0.61 (0.28)*
≥7h−1.57 (0.66)*−1.45 (0.56)**−1.41 (0.56)*−0.42 (0.33)−0.48 (0.28)−0.47 (0.28)
P for interaction0.550.900.540.170.400.75

Partial regression coefficients were calculated after adjusting MAP for sex and age. Abbreviations: MAP, mean arterial pressure; SE, standard error of mean; RMSSD, root mean square of successive RR interval differences; HF, high frequency.

*P < 0.05; **P < 0.01.

DISCUSSION

Low RMSSD and HF were significantly associated with an increased home MAP in the morning rather than in the evening, and these associations were independent of sex, age, BMI, smoking, alcohol drinking, use of hypertensive agents, diabetes, sleeping duration, snoring, physical activity, and SF-8 MCS. Our validation study suggested that 5-minute HRV parameters measured from fingertip pulse recordings at health checkups were moderately associated with HRV parameters from 24-hour Holter ECG recordings. Thus, HRV assessed from pulse recordings was fairly representative of HRV over a 24-hour period and influenced home MAP in the morning regardless of daily activities and sleep conditions.

Both home BP monitoring and ambulatory BP monitoring (ABPM) are thought to be superior to BP measurement in the clinic for the prediction of future CVD.22 In addition, day-to-day variability of home BP was shown to independently predict cardiovascular events and total mortality.23 Furthermore, self-monitoring of BP in conjunction with other interventions (medication titration, patient education, and lifestyle counseling) in patients with hypertension led to a clinically significant BP reduction that persisted for at least 12 months.24 Therefore, it is important to identify the factors influencing home BP to obtain better BP control throughout the day. Using ABPM, Ohira et al.25 found that habitual alcohol intake raised morning BP, and this was linked with increased sympathetic activity during sleep.

HRV is a noninvasive tool to evaluate autonomic nervous system function.15 SDNN and LF are considered to represent overall sympathetic and parasympathetic activity. RMSSD and HF are indices of parasympathetic activity, and LF/HF is an indicator of sympathovagal balance. We confirmed that there are daily variations in autonomic nervous activity in the present study; there was a relative increase in parasympathetic activity in the evening and during sleep, which contributed to a lower evening BP at home. Consequently, a relative increase in parasympathetic tone and a decrease in sympathetic tone in the evening might weaken the associations between MAP and the HRV indices.

An association between autonomic nervous function and hypertension can be explained in part by the arterial baroreflex.1 The baroreflex is believed to be a short-term controller of BP reflecting parasympathetic activation and sympathetic inhibition; however, Lohmeier et al.3 reported that the baroreflex has a long-term effect on BP. Hesse et al.2 demonstrated that baroreflex sensitivity was inversely associated with the means of 24-hour or daytime MAP. Consistent with this, the Rotterdam Study showed that impaired baroreflex sensitivity reflected arterial stiffness and orthostatic BP changes in older adults.4

In the present study, we focused on MAP because of the role of baroreflex sensitivity in regulating MAP, as mentioned above. When systolic and diastolic BPs were analyzed separately, RMSSD and HF were inversely associated with morning home diastolic BP, rather than morning home systolic BP. This finding is consistent with a previous study that showed baroreflex sensitivity was highly correlated with diastolic but not systolic BP.2

The associations between RMSSD, HF, and MAP were obvious in individuals <65 years old. Since baroreflex sensitivity depends on age,26 the impact of reduced parasympathetic activity on baroreflex sensitivity might be eliminated in the elderly.

The strength of the present study is that it included standardized HRV assessment and a large number of lifestyle variables. Nonetheless, several potential limitations should be noted. First, because this was a cross-sectional study, we could not confirm a causal relationship between decreased HRV and increased MAP. Second, we could not exclude the effects of agents that act on autonomic cardiac function, such as antihypertensive drugs. Our data on medication use were based on a self-administered questionnaire that did not involve information on specific classes of drugs. Furthermore, sensitivity analysis verified that the association between HRV and MAP did not change when patients on antihypertensive drugs were excluded. Third, because our subjects were volunteers from a single community, our findings may not be representative of the general Japanese population.

CONCLUSIONS

Low HF and RMSSD due to impaired parasympathetic nervous activity were associated with increased MAP levels in the morning rather than in the evening. Our findings suggest that factors which induce sympathovagal imbalance (physical inactivity, insomnia, and socioeconomic stress) should be eliminated or reduced to lower home BP in the morning.

SUPPLEMENTARY MATERIAL

Supplementary data are available at American Journal of Hypertension online.

ACKNOWLEDGMENTS

This study was supported, in part, by Grants-in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, Science, and Technology of Japan (Grants-in-Aid for Research B, no. 22390134 in 2010–2012 and 25293142 from 2013, Grants-in-Aid for Young Scientists (B), no. 25860443 and 25860441 from 2013, and Grant-in-Aid for Research C, no. 26460767) and Health and Labor Sciences Research Grants from the Ministry of Health, Welfare, and Labor, Japan (Comprehensive Research on Life-Style Related Diseases including Cardiovascular Diseases and Diabetes Mellitus, no. 201021038A in 2010–2012). We thank Takuma Akita, other staff and participants of the Toon Health Study and the municipal authorities, officers, and health professionals of Toon City for their valuable contributions.

DISCLOSURE

The authors declared no conflict of interest.

REFERENCES

1.

Carthy
ER
.
Autonomic dysfunction in essential hypertension: a systematic review
.
Ann Med Surg (Lond)
 
2014
;
3
:
2
7
.

2.

Hesse
C
,
Charkoudian
N
,
Liu
Z
,
Joyner
MJ
,
Eisenach
JH
.
Baroreflex sensitivity inversely correlates with ambulatory blood pressure in healthy normotensive humans
.
Hypertension
 
2007
;
50
:
41
46
.

3.

Lohmeier
TE
,
Irwin
ED
,
Rossing
MA
,
Serdar
DJ
,
Kieval
RS
.
Prolonged activation of the baroreflex produces sustained hypotension
.
Hypertension
 
2004
;
43
:
306
311
.

4.

Mattace-Raso
FU
,
van den Meiracker
AH
,
Bos
WJ
,
van der Cammen
TJ
,
Westerhof
BE
,
Elias-Smale
S
,
Reneman
RS
,
Hoeks
AP
,
Hofman
A
,
Witteman
JC
.
Arterial stiffness, cardiovagal baroreflex sensitivity and postural blood pressure changes in older adults: the Rotterdam Study
.
J Hypertens
 
2007
;
25
:
1421
1426
.

5.

Pickering
TG
,
Miller
NH
,
Ogedegbe
G
,
Krakoff
LR
,
Artinian
NT
,
Goff
D
;
American Heart Association; American Society of Hypertension; Preventive Cardiovascular Nurses Association
.
Call to action on use and reimbursement for home blood pressure monitoring: executive summary: a joint scientific statement from the American Heart Association, American Society of Hypertension, and Preventive Cardiovascular Nurses Association
.
Hypertension
 
2008
;
52
:
1
9
.

6.

Parati
G
,
Omboni
S
,
Bilo
G
.
Why is out-of-office blood pressure measurement needed
?
Hypertension
 
2009
;
54
:
181
187
.

7.

Liao
D
,
Carnethon
M
,
Evans
GW
,
Cascio
WE
,
Heiss
G
.
Lower heart rate variability is associated with the development of coronary heart disease in individuals with diabetes: the Atherosclerosis Risk in Communities (ARIC) study
.
Diabetes
 
2002
;
51
:
3524
3531
.

8.

Singh
JP
,
Larson
MG
,
Tsuji
H
,
Evans
JC
,
O’Donnell
CJ
,
Levy
D
.
Reduced heart rate variability and new-onset hypertension: insights into pathogenesis of hypertension: the Framingham Heart Study
.
Hypertension
 
1998
;
32
:
293
297
.

9.

Schroeder
EB
,
Chambless
LE
,
Liao
D
,
Prineas
RJ
,
Evans
GW
,
Rosamond
WD
,
Heiss
G
;
Atherosclerosis Risk in Communities (ARIC) study
.
Diabetes, glucose, insulin, and heart rate variability: the Atherosclerosis Risk in Communities (ARIC) study
.
Diabetes Care
 
2005
;
28
:
668
674
.

10.

Tsuji
H
,
Larson
MG
,
Venditti
FJ
Jr
,
Manders
ES
,
Evans
JC
,
Feldman
CL
,
Levy
D
.
Impact of reduced heart rate variability on risk for cardiac events. The Framingham Heart Study
.
Circulation
 
1996
;
94
:
2850
2855
.

11.

Saito
I
,
Hitsumoto
S
,
Maruyama
K
,
Nishida
W
,
Eguchi
E
,
Kato
T
,
Kawamura
R
,
Takata
Y
,
Onuma
H
,
Osawa
H
,
Tanigawa
T
.
Heart rate variability, insulin resistance, and insulin sensitivity in Japanese adults: the Toon Health Study
.
J Epidemiol
 
2015
;
25
:
583
591
.

12.

Saito
I
,
Maruyama
K
,
Eguchi
E
,
Kato
T
,
Kawamura
R
,
Takata
Y
,
Onuma
H
,
Osawa
H
,
Tanigawa
T
.
Low heart rate variability and sympathetic dominance modifies the association between insulin resistance and metabolic syndrome—the Toon Health Study
.
Circ J
 
2017
;
81
:
1447
1453
.

13.

Saito
I
,
Hitsumoto
S
,
Maruyama
K
,
Eguchi
E
,
Kato
T
,
Okamoto
A
,
Kawamura
R
,
Takata
Y
,
Nishida
W
,
Nishimiya
T
,
Onuma
H
,
Osawa
H
,
Tanigawa
T
.
Impact of heart rate variability on C-reactive protein concentrations in Japanese adult nonsmokers: the Toon Health Study
.
Atherosclerosis
 
2016
;
244
:
79
85
.

14.

Mori
H
,
Saito
I
,
Eguchi
E
,
Maruyama
K
,
Kato
T
,
Tanigawa
T
.
Heart rate variability and blood pressure among Japanese men and women: a community-based cross-sectional study
.
Hypertens Res
 
2014
;
37
:
779
784
.

15.

Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology
.
Heart rate variability. Standards of measurement, physiological interpretation, and clinical use. Task force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology
.
Eur Heart J
 
1996
;
17
:
354
381
.

16.

Jennum
P
,
Hein
HO
,
Suadicani
P
,
Gyntelberg
F
.
Cardiovascular risk factors in snorers. A cross-sectional study of 3,323 men aged 54 to 74 years: the Copenhagen Male Study
.
Chest
 
1992
;
102
:
1371
1376
.

17.

Nagayoshi
M
,
Tanigawa
T
,
Yamagishi
K
,
Sakurai
S
,
Kitamura
A
,
Kiyama
M
,
Okada
T
,
Maeda
K
,
Ohira
T
,
Imano
H
,
Sato
S
,
Iso
H
;
CIRCS Investigators
.
Self-reported snoring frequency and incidence of cardiovascular disease: the Circulatory Risk in Communities Study (CIRCS)
.
J Epidemiol
 
2012
;
22
:
295
301
.

18.

Ishikawa-Takata
K
,
Naito
Y
,
Tanaka
S
,
Ebine
N
,
Tabata
I
.
Use of doubly labeled water to validate a physical activity questionnaire developed for the Japanese population
.
J Epidemiol
 
2011
;
21
:
114
121
.

19.

Ainsworth
BE
,
Haskell
WL
,
Whitt
MC
,
Irwin
ML
,
Swartz
AM
,
Strath
SJ
,
O’Brien
WL
,
Bassett
DR
Jr
,
Schmitz
KH
,
Emplaincourt
PO
,
Jacobs
DR
Jr
,
Leon
AS
.
Compendium of physical activities: an update of activity codes and MET intensities
.
Med Sci Sports Exerc
 
2000
;
32
:
S498
S504
.

20.

Tokuda
Y
,
Okubo
T
,
Ohde
S
,
Jacobs
J
,
Takahashi
O
,
Omata
F
,
Yanai
H
,
Hinohara
S
,
Fukui
T
.
Assessing items on the SF-8 Japanese version for health-related quality of life: a psychometric analysis based on the nominal categories model of item response theory
.
Value Health
 
2009
;
12
:
568
573
.

21.

Hallman
DM
,
Srinivasan
D
,
Mathiassen
SE
.
Short- and long-term reliability of heart rate variability indices during repetitive low-force work
.
Eur J Appl Physiol
 
2015
;
115
:
803
812
.

22.

Shimbo
D
,
Abdalla
M
,
Falzon
L
,
Townsend
RR
,
Muntner
P
.
Role of ambulatory and home blood pressure monitoring in clinical practice: a narrative review
.
Ann Intern Med
 
2015
;
163
:
691
700
.

23.

Stergiou
GS
,
Ntineri
A
,
Kollias
A
,
Ohkubo
T
,
Imai
Y
,
Parati
G
.
Blood pressure variability assessed by home measurements: a systematic review
.
Hypertens Res
 
2014
;
37
:
565
572
.

24.

Tucker
KL
,
Sheppard
JP
,
Stevens
R
,
Bosworth
HB
,
Bove
A
,
Bray
EP
,
Earle
K
,
George
J
,
Godwin
M
,
Green
BB
,
Hebert
P
,
Hobbs
FDR
,
Kantola
I
,
Kerry
SM
,
Leiva
A
,
Magid
DJ
,
Mant
J
,
Margolis
KL
,
McKinstry
B
,
McLaughlin
MA
,
Omboni
S
,
Ogedegbe
O
,
Parati
G
,
Qamar
N
,
Tabaei
BP
,
Varis
J
,
Verberk
WJ
,
Wakefield
BJ
,
McManus
RJ
.
Self-monitoring of blood pressure in hypertension: a systematic review and individual patient data meta-analysis
.
PLoS Med
 
2017
;
14
:
e1002389
.

25.

Ohira
T
,
Tanigawa
T
,
Tabata
M
,
Imano
H
,
Kitamura
A
,
Kiyama
M
,
Sato
S
,
Okamura
T
,
Cui
R
,
Koike
KA
,
Shimamoto
T
,
Iso
H
.
Effects of habitual alcohol intake on ambulatory blood pressure, heart rate, and its variability among Japanese men
.
Hypertension
 
2009
;
53
:
13
19
.

26.

Laitinen
T
,
Hartikainen
J
,
Vanninen
E
,
Niskanen
L
,
Geelen
G
,
Länsimies
E
.
Age and gender dependency of baroreflex sensitivity in healthy subjects
.
J Appl Physiol (1985)
 
1998
;
84
:
576
583
.

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic.oup.com/pages/standard-publication-reuse-rights)

Supplementary data