| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
(Hypertension. 2007;50:143.)
© 2007 American Heart Association, Inc.
Original Articles |
From the Inquiry Methodology Program (K.K.), Indiana University, Bloomington; Division of Nephrology (R.P.L., R.A.), Department of Medicine, Indiana University School of Medicine, Indianapolis; and Richard L. Roudebush VA Medical Center (R.P.L., R.A.), Indianapolis, Ind.
Correspondence to Rajiv Agarwal, VA Medical Center, 111N, 1481 West 10th St, Indianapolis IN 46202. E-mail ragarwal{at}iupui.edu
| Abstract |
|---|
|
|
|---|
Key Words: hemodialysis hypertension ambulatory blood pressure monitoring circadian rhythms statistical modeling
| Introduction |
|---|
|
|
|---|
In the renoprival state, gain of volume over an interdialytic period leads to increase in BP. The traditional cosinor model requires midline estimating statistic of rhythm (MESOR), defined as the average value of the rhythmic function fitted to the data to be flat.14 Gain in volume over an interdialytic interval may cause a steady increase in systemic arterial pressure between 2 dialysis treatments, and the above assumption of a flat MESOR may no longer be tenable. Thus, a more complex model that includes terms for cosinor, as well as a straight-line change, would be more appropriate to describe the data. A trended cosinor model that includes both parameters has not been developed but can be useful to evaluate the impact of interventions on not only the mean BP but also on the rate of change in BP between dialyses and the amplitude of these variations.
The purpose of this report is to describe the development of a trended cosinor model. We then examine the association of the impact of antihypertensive drug therapy with BP patterns in hemodialysis patients using this model.
| Methods |
|---|
|
|
|---|
Subjects
Patients
18 years of age who had been on chronic hemodialysis for >3 months and were free of vascular, infectious, or bleeding complication within 1 month were enrolled in the study. Those who missed
2 hemodialysis treatments over 1 month, used illicit drugs, and had chronic atrial fibrillation or body mass index of
40 kg/m2 were excluded. Patients who had a change in dry weight or change in antihypertensive drugs within 2 weeks were also excluded. Presence or absence of hypertension was not a selection criterion. All of the patients underwent standard dialysis 3 times a week.
Anthropometric and demographic characteristics and antihypertensive medications actually taken by the patient were recorded. The study was approved by the institutional review board of Indiana University and the research and development committee of the Roudebush VA Medical Center, and all of the subjects gave written informed consent.
Ambulatory BP Monitoring
Ambulatory BP monitoring was performed after the midweek hemodialysis session for 44 hours. Ambulatory BPs were recorded every 20 minutes during the day (6:00 AM to 10:00 PM) and every 30 minutes during the night (10:00 PM to 6:00 AM) using a Spacelab 90207 ABP monitor (SpaceLabs Medical Inc) in the nonaccess arm, as done previously.13 Patients with <16 ambulatory BP recordings were excluded, because pattern recognition was not possible with a limited number of recordings. The remaining 136 patients had a combined 12 750 BP measurements. These data were exported to a relational database to allow for data management, as well as centering the time to that elapsed after dialysis using standard programming tools.
Analysis
Original oscillometric data from each BP series were first synchronized for each individual by recomputing all of the times of sampling in hours from the end of dialysis to avoid differences among subjects dialyzing on different dialysis shifts. We used both clock time and resynchronized time postdialysis for model fit. We compared models that used time postdialysis with models that used clock time when modeling hemodynamic parameters over an interdialytic interval. In addition, we used a composite model where we used clock time for cosinor and time elapsed after dialysis for linear trends.
Statistical Methods
To detect the presence of any diurnal pattern in hemodynamic measures in an interdialytic interval, we used the cosinor model.14,15 This method entails fitting an oscillating curve to temporal hemodynamic variables, such as BP and heart rate, using some specified (eg, 24-hour) periodicity. The hemodynamic variables were obtained for
44 hours after hemodialysis. The cosinor model first considered to be describing the rhythmic cycle can be expressed as y = b0+b1xCos[(2
/24)t]+b2xSin[(2
/24)t] where y represents the observed systolic BP, diastolic BP, pulse pressure or heart rate; b0, b1, and b2 are regression coefficients; and t represents time (eg, elapsed time after dialysis or clock time). The constant 2
/24 represents the 24-hour periodicity of BP. The coefficient b0 represents the 24-hour rhythm-adjusted mean BP also called the MESOR, defined as the average value of hemodynamic measure (BP, pulse pressure, or heart rate). The regression coefficients b1 and b2 are the coefficients for the cosine and sine component, respectively, and collectively describe the amplitude of the cosine curve, which is defined as amplitude=
|
|
Although the cosinor model described is adequate for many biological processes that cycle over a 24-hour period, a fundamental issue not addressed with the cosinor model arises. This issue is that, over an interdialytic interval, not only are there oscillations in BP, there is also a systematic positive increasing trend because of the interdialytic volume expansion (see observed and modeled BP and heart rate, eg, in Figure 1). Thus, it was necessary to add an additional component to the cosinor model that allowed for the systematic change in BP in addition to the oscillation because of the BP cycle, which is modeled with the cosinor model. This generalized cosinor model, termed the trended cosinor model, is given as follows:
|
|
|
which is the same as that given previously with the additional parameter accounting for any systematic linear change over time. Thus, the change model explicitly considers 2 types of change in a unified manner: change that has a systematic linear component and change that oscillates.
In the trended cosinor model, the meaning of the amplitude is slightly modified to the extent of rhythmic change over and above the linear trend. The meaning of b0 changes considerably from the MESOR in the standard cosinor model to the predicted value of BP1 at t=0 (the 1 is necessary because cosine[0]=1). The mean BP is now conditional on the value of time in the generalized cosinor model because of the (assumed) non-0 slope.
When models were nested, comparisons of goodness of fit between models were made by testing the improvement in the log-likelihood ratio using a
2 test. When models were not nested, comparisons of goodness of fit were make by comparing Akaike Information Criteria and the Bayesian Information Criteria. For Akaike Information Criteria and Bayesian Information Criteria comparisons, smaller numbers indicate a better model fit.
The mixed-effects models were fit in R using the nonlinear mixed-effects package (nlme).16,17 Full information maximum likelihood was used for parameter estimation, and there was a random component associated with each fixed effect. The delta method with first-order Taylor expansion was used to estimate the SEs (and, thus, confidence limits and P values) of functions of fixed effect parameters (eg, the amplitude).18
After development of a model where the parameters were unconditional, conditional models were developed where the number of BP medications (0, 1, 2, 3, 4, or more) was used to explain interindividual differences in change.
| Results |
|---|
|
|
|---|
|
Table 2 shows the sleepawake BP and heart rate as a function of BP medications. Systolic BP fell 2.3 mm Hg and diastolic BP by twice as much; thus, the pulse pressure increased by 2.3 mm Hg. There was a strong direct relationship between antihypertensive drug use and systolic BP. The diurnal fall in BP was influenced by the number of antihypertensive drugs. Those who took a greater number of antihypertensive drugs had, in fact, a paradoxical increase in systolic and diastolic BP.
|
Table 3 shows the modeled linear and nonlinear effects for systolic BP. The average systolic ambulatory BP was 129.7 mm Hg (model 1). Just after dialysis, systolic BP was 124.4 mm Hg, and BP increased linearly at a rate of 0.26 mm Hg per hour elapsed postdialysis (model 2). The model fit was considerably improved with the linear term, as judged by the improvement in the log-likelihood function. Fitting the traditional cosinor model (model 3) improved the fit compared with model 1 (likelihood ratio improved) but deteriorated the fit compared with model 2 (higher Akaike Information Criteria and Bayesian Information Criteria). That is, the cosinor model was better than the no-slope model, but the model with a slope was better than the cosinor model. However, the straight-line model treats the oscillation about the line as error, but such oscillations may be important from a clinical perspective. The trended cosinor model discussed previously includes simultaneously oscillatory terms from the cosinor model and a slope from the straight-line change model. The trended cosinor model (model 4) was superior to any of the other models. Although sine and cosine coefficients were not significant, removal of either of these components deteriorated the model fit (data not shown). That is to say, although the mean of these coefficients did not differ significantly from 0, there was interindividual variation in these change coefficients that lead to a better fit when the model allowed the terms to be fit and to vary among the individuals. Thus, allowing for variation in systolic BP in a linear and oscillating fashion best described the data. Model 5 represents fitting the linear trend to time elapsed after dialysis and cosinor to clock time. Model fit was similar to model 4; therefore, we decided to use model 4 in all of the subsequent analyses.
|
To examine the impact of dialysis time on BP changes, we reanalyzed the data first by clock time starting at midnight and then centering the time such that the clock started at the end of dialysis for each patient. Table 4 shows that the model fit improved considerably for the straight-line model when time was counted from that elapsed after dialysis. In the case of the cosinor model, the clock time allowed for a slightly better model fit compared with the time elapsed after dialysis. (This inconsistency is why we fit model 5). The trended cosinor model had the best fit when time was counted after dialysis.
|
Circadian parameters for other BP components and heart rate are shown in Table 5. The amplitude of change in systolic BP was <1 mm Hg and was marginally statistically significant (P=0.07). On the other hand, increase in slope of 0.26 mm Hg/h was highly significant. For diastolic BP, the rate of change was approximately half that seen for systolic BP, and the amplitude of variation was twice as large as systolic BP. Pulse pressure increased at a rate that would lead to 6-mm Hg amplification over the 2-day interdialytic interval. No significant linear change was seen in heart rate. The peak-to-trough variation in amplitude in heart rate would average 5 mm Hg based on our data.
|
Figure 2 shows the representation of the modeled BP and heart rate. Undulations in systolic BP are gentle and nonsignificant, but a steady increase is notable. Diastolic BP has more notable undulations and a less steep trajectory compared with systolic BP. Heart rate has a flat trajectory and circadian variation that is most notable.
|
The impact of number of antihypertensive drugs on BP characteristics is shown in Table 6. Each medication increment was associated with an increase in the systolic BP of 6.3 mm Hg and diastolic BP of 2.7 mm Hg. Significant blunting in the rate of change of systolic and diastolic BP was seen with antihypertensive medication intake. No such blunting in pulse pressure was seen with medications. The amplitude of variation in systolic BP followed a U-shaped curve for systolic BP, reduction for diastolic BP, and increase for pulse pressure (Figure 3). The lower limit of 95% CI for patients on 2 or 3 antihypertensive medications cross 0 in the case of systolic BP. Therefore, these variations were not significant. There were significant differences between mean amplitudes for each parameter (systolic pulse pressure, diastolic pulse pressure, or heart rate) for the number of antihypertensive drugs.
|
|
Figure 4 shows the modeled curves of circadian systolic and diastolic BP by intake of antihypertensive agents. The linear slope in systolic BP is more than twice that of diastolic BP. Furthermore, the linear slope of change becomes flatter with increasing antihypertensive medication intake. The U-shaped nature in variation in amplitude is evident for systolic BPs, and blunting in variation in diastolic BP with increasing number of medications is seen.
|
| Discussion |
|---|
|
|
|---|
The clinical significance of the trended cosinor model is illustrated by comparing BP and heart rate data in 136 hemodialysis patients. For example, heart rate showed no straight-line trend and demonstrated substantial amplitude, but systolic BP had a significant straight-line trend and little amplitude (Table 4 and Figure 1). In contrast, diastolic BP had straight-line trend that was less steep than systolic BP and amplitude that was twice that of systolic BP. Thus, hemodialysis impacts systolic BP, diastolic BP, and heart rate differently. It is possible that volume factors predominantly affect systolic BP, whereas vascular tone may impact diastolic BP, potentially explaining the findings of our study.
In addition, we could detect important differences related to the number of drugs used to treat hypertension. Although the recognition of greater average systolic and diastolic BP with more medications is well documented, the analysis of the impact of the number of drugs and BP trends has not previously been possible.19 The increasing number of antihypertensive drugs blunted the rate of rise in systolic BP and diastolic BP in the interdialytic period. In other words, those who received greater numbers of antihypertensive medications experienced less interdialytic BP increase; this increment was seen more for systolic than diastolic BP. The amplitude of circadian changes in systolic BP was no different from 0 for those who received 2 or 3 drugs. In contrast, those who were on no drugs, 1 drug, or >3 drugs had non-0 amplitudes (Figure 2). For diastolic BP, reduction in amplitude was seen for increasing number of medications, such that pulse pressure amplitude increased with increasing number of medications. We were able to detect these patterns of changes as a direct result of the development of this statistical model.
Rodby et al15 have previously described circadian rhythm in systolic and diastolic BP and heart rate in hemodialysis patients. In fact, the authors found that, despite interdialytic weight gain, there was no apparent increase in BP from day 1 to day 2 in dialysis patients. Fitting the cosinor model that assumes a flat MESOR may have obscured the linear trend in their data. Evidently, the cosinor model fit is reasonable when the linear trend is ignored in part because the misspecification recovers some of the linear change. Comparison of the amplitude of variation and linear changes in systolic BP using the trended cosinor model (model 4; Table 2 and Figure 1) demonstrates that the linear change is strong and overwhelms the relatively small variation in amplitude. Although population averaged circadian variation in diastolic BP and heart rate was detected, we found marginal circadian variation in systolic BP (P=0.07) despite a much larger study compared with Rodby et al.15 This may be related to fitting fundamentally different models, as noted above, but may in part be related to the number of antihypertensive drugs. Hemodialysis patients using 0 to 1 antihypertensive drugs constituted 75% of the population reported by Rodby et al,15 whereas these patients composed only 45% of our sample. Analysis of circadian variation in systolic BP was seen in patients taking 0 to 1 drugs but not in those taking 2 to 3 drugs. Thus, dampening of circadian variation in systolic BP seen in our sample could be related to patient and model differences.
Pulse pressure has emerged as 1 of the strongest predictors of total mortality in hemodialysis patients.20,21 An increase in pulse pressure over the interdialytic interval was noted amounting to
6 mm Hg over 44 hours. This increase in pulse pressure may reach 10 mm Hg over the weekend. These data may explain, in part, why mortality is seen to be greatest just before the first dialysis of the week.22 Increasing the number of antihypertensive drugs did not blunt the pulse pressure. Although these are observational data, the model developed can be used to study the impact of novel drugs or reduction in dry weight in randomized trials on pulse pressure and/or other hemodynamic measures.
Perspectives
One practical application of this model could be to analyze patterns of BP changes with various antihypertensive drugs. Recent data in patients without chronic kidney disease demonstrate differences between ß-blocker- and angiotensin-converting enzyme inhibitorbased therapies on central pressures, which were thought to be of prognostic importance.23 Whether prognostic information can be derived from analysis of patterns of ambulatory BP changes obtained with the trended cosinor model remains to be seen. Drugs such as cyclosporine and erythropoietin often lead to an increase in BP in patients with chronic kidney disease.24,25 Whether changes in amplitude or slopes predate the development of absolute increase in BP may shed light on the pathophysiology of these disorders. Future studies will analyze the impact of changes in these parameters in randomized, controlled trials of drugs, diets, and dry weight changes in hemodialysis patients. The development of the trended cosinor statistical model allows us to study the diagnostic and prognostic significance of the hemodynamic parameters in hemodialysis patients.
| Acknowledgments |
|---|
Sources of Funding
This work was supported by the National Institutes of Health (National Institutes of Health/National Institute of Diabetes and Digestive and Kidney Diseases 5RO1 DK062030-04).
Disclosures
None.
Received March 28, 2007; first decision April 16, 2007; accepted April 21, 2007.
| References |
|---|
|
|
|---|
2. Agarwal R, Andersen MJ. Prognostic importance of clinic and home blood pressure recordings in patients with chronic kidney disease. Kidney Int. 2006; 69: 406411.[CrossRef][Medline] [Order article via Infotrieve]
3. Agarwal R. Role of home blood pressure monitoring in hemodialysis patients. Am J Kidney Dis. 1999; 33: 682687.[Medline] [Order article via Infotrieve]
4. Perloff D, Sokolow M, Cowan R. The prognostic value of ambulatory blood pressures. JAMA. 1983; 249: 27922798.
5. Townsend RR, Ford V. Ambulatory blood pressure monitoring: coming of age in nephrology. J Am Soc Nephrol. 1996; 7: 22792287.[Abstract]
6. Thompson AM, Pickering TG. The role of ambulatory blood pressure monitoring in chronic and end-stage renal disease. Kidney Int. 2006; 70: 10001007.[CrossRef][Medline] [Order article via Infotrieve]
7. Lurbe E, Redon J. Assessing ambulatory blood pressure in renal diseases: facts and concerns. Nephrol Dial Transplant. 1999; 14: 25642568.
8. Nielsen FS, Rossing P, Bang LE, Svendsen TL, Gall MA, Smidt UM, Parving HH. On the mechanisms of blunted nocturnal decline in arterial blood pressure in NIDDM patients with diabetic nephropathy. Diabetes. 1995; 44: 783789.[Abstract]
9. Portaluppi F, Montanari L, Massari M, Di Chiara V, Capanna M. Loss of nocturnal decline of blood pressure in hypertension due to chronic renal failure. Am J Hypertens. 1991; 4: 2026.[Medline] [Order article via Infotrieve]
10. Agarwal R, Andersen MJ. Correlates of systolic hypertension in patients with chronic kidney disease. Hypertension. 2005; 46: 514520.
11. Coomer RW, Schulman G, Breyer JA, Shyr Y. Ambulatory blood pressure monitoring in dialysis patients and estimation of mean interdialytic blood pressure [see comments]. Am J Kidney Dis. 1997; 29: 678684.[Medline] [Order article via Infotrieve]
12. Rahman M, Griffin V, Heyka R, Hoit B. Diurnal variation of blood pressure; reproducibility and association with left ventricular hypertrophy in hemodialysis patients. Blood Press Monit. 2005; 10: 2532.[CrossRef][Medline] [Order article via Infotrieve]
13. Agarwal R. Supervised atenolol therapy in the management of hemodialysis hypertension. Kidney Int. 1999; 55: 15281535.[CrossRef][Medline] [Order article via Infotrieve]
14. Nelson W, Tong YL, Lee JK, Halberg F. Methods for cosinor-rhythmometry. Chronobiologia. 1979; 6: 305323.[Medline] [Order article via Infotrieve]
15. Rodby RA, Vonesh EF, Korbet SM. Blood pressures in hemodialysis and peritoneal dialysis using ambulatory blood pressure monitoring. Am J Kidney Dis. 1994; 23: 401411.[Medline] [Order article via Infotrieve]
16. Pinheiro JC, Bates DM. Mixed-effects models in S and S-PLUS. New York, NY: Springer; 2000.
17. Pinheiro J, Bates D, DeloRoy S, Sarkar D. nlme: Linear and Nonlinear Mixed Effects Models [computer program]. Available at: www.r-project.org. Accessed May 5, 2007.
18. Oehlert GW. A note on the delta method. Am Stat. 1992; 46: 2729.[CrossRef]
19. Agarwal R, Nissenson AR, Batlle D, Coyne DW, Trout JR, Warnock DG. Prevalence, treatment, and control of hypertension in chronic hemodialysis patients in the United States. Am J Med. 2003; 115: 291297.[CrossRef][Medline] [Order article via Infotrieve]
20. Klassen PS, Lowrie EG, Reddan DN, DeLong ER, Coladonato JA, Szczech LA, Lazarus JM, Owen WF Jr. Association between pulse pressure and mortality in patients undergoing maintenance hemodialysis. JAMA. 2002; 287: 15481555.
21. Tozawa M, Iseki K, Iseki C, Takishita S. Pulse pressure and risk of total mortality and cardiovascular events in patients on chronic hemodialysis. Kidney Int. 2002; 61: 717726.[CrossRef][Medline] [Order article via Infotrieve]
22. Bleyer AJ, Russell GB, Satko SG. Sudden and cardiac death rates in hemodialysis patients. Kidney Int. 1999; 55: 15531559.[CrossRef][Medline] [Order article via Infotrieve]
23. Williams B, Lacy PS, Thom SM, Cruickshank K, Stanton A, Collier D, Hughes AD, Thurston H, ORourke M. Differential impact of blood pressure-lowering drugs on central aortic pressure and clinical outcomes: principal results of the Conduit Artery Function Evaluation (CAFE) study. Circulation. 2006; 113: 12131225.
24. Berns JS, Rudnick MR, Cohen RM, Bower JD, Wood BC. Effects of normal hematocrit on ambulatory blood pressure in epoetin- treated hemodialysis patients with cardiac disease. Kidney Int. 1999; 56: 253260.[CrossRef][Medline] [Order article via Infotrieve]
25. Agarwal R. Hypertension in chronic kidney disease and dialysis: pathophysiology and management. Cardiol Clin. 2005; 23: 237248.[CrossRef][Medline] [Order article via Infotrieve]
This article has been cited by other articles:
![]() |
R. Agarwal and R. P. Light GFR, proteinuria and circadian blood pressure Nephrol. Dial. Transplant., February 27, 2009; (2009) gfp074v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Agarwal, T. Metiku, G. G. Tegegne, R. P. Light, Z. Bunaye, D. M. Bekele, and K. Kelley Diagnosing Hypertension by Intradialytic Blood Pressure Recordings Clin. J. Am. Soc. Nephrol., September 1, 2008; 3(5): 1364 - 1372. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Agarwal and R. P. Light Arterial stiffness and interdialytic weight gain influence ambulatory blood pressure patterns in hemodialysis patients Am J Physiol Renal Physiol, February 1, 2008; 294(2): F303 - F308. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Hypertension Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2007 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |