Original Research ARTICLE
Home Monitoring of Heart Rate as a Predictor of Imminent Cardiovascular Events
- Research Unit of Internal Medicine, Medical Research Center Oulu, Oulu University Hospital, University of Oulu, Oulu, Finland
Introduction: Previous studies have documented that day-to-day variability of heart rate (HR) has prognostic significance for cardiovascular (CV) events in general population. It is unknown how HR dynamics variate before imminent CV event in patients with coronary artery disease (CAD). Our aim was to study day-to-day variation in HR dynamics before the occurrence of CV event in patients with initially stable CAD.
Methods: Forty-four patients with angiographically documented CAD from ARTEMIS study measured R-R intervals on a weekly basis at home for 2 years. Home measurements were performed in controlled conditions (3 min at supine and sitting) 1–2 times per week. Eleven patients had a CV event (7 acute coronary syndromes, 1 cardiac death, 2 new onset of arrhythmia needing hospitalization and 1 stroke), which occurred 11 ± 7 months after enrolment. Mean R-R interval was analyzed prospectively from the home measurements. For the patients with new CV event, average, and standard deviation (SD) of the mean R-R interval over 8 weeks preceding the CV event were calculated. For the patients without new CV event, corresponding period was determined by the median follow-up at the occurrence of new CV event.
Results: There were no differences in the mean R-R interval analyzed over 8 weeks between the patients with and without new CV event. The variability of mean R-R interval over 8 weeks was greater in the patients with new CV event compared to the patients without new CV event at the supine (95 ± 34 vs. 59 ± 26 ms, p < 0.001) and sitting positions (92 ± 28 vs. 62 ± 24 ms, p < 0.001).
Conclusion: Day-to-day variability of mean R-R interval is greater before the new CV event in CAD patients suggesting to a more unstable cardiac autonomic regulation preceding these events.
Previous studies have demonstrated prognostic significance of day-to-day blood pressure and heart rate (HR) variability for cardiovascular (CV) events in general population (Kikuya et al., 2008; Johansson et al., 2012). Kikuya et al. (2008) showed that increased day-to-day variability of systolic blood pressure and HR was associated with greater CV mortality during long-term follow-up even after adjustment for traditional CV risk factors. Also Johansson et al. (2012) reported that higher variability of morning day-to-day blood pressure and HR predicted fatal and non-fatal CV events in a middle-aged general population, whereas the variability of evening day-by-day blood pressure and HR did not.
Among patient populations, increased day-to-day or visit-to-visit systolic blood pressure variability has been associated with future risk for CV events in patients with hypertension and recent transient ischemic attack or stroke (Chowdhury et al., 2014; Yu et al., 2014; Webb et al., 2018). However, the association between day-to-day variability of HR and CV events is unclear in stable coronary artery disease (CAD) patients. It is well-known that HR variability is altered among patients with ischemic heart disease compared with healthy subjects (Airaksinen et al., 1987; Huikuri et al., 1994; Bigger et al., 1995). Secondly, Adamson et al. (2004) demonstrated that HR variability, measured continuously from an implanted cardiac resynchronization device, declined a median of 16 days before the hospitalization in patients with New York Heart Association class III or IV heart failure. However, the previous studies have not reported to the role of day-to-day variation in mean HR or HR variability before imminent CV event in patients with initially stable CAD during the current treatment era. Therefore, the purpose of present study was to examine the day-to-day HR dynamics in patients with stable CAD before imminent CV event.
Materials and Methods
Subjects and Study Protocol
The work described here is part of a larger ARTEMIS (Innovation to Reduce Cardiovascular Complications of Diabetes at the Intersection) study taking place in the Division of Cardiology at Oulu University Hospital at 2007–2014 (Oulu, Finland). The ARTEMIS study is registered at ClinicalTrials.gov, Record 1539/31/06 identifier number NCT01426685. The study was performed according to the Declaration of Helsinki, the local research ethics committee of the Northern Ostrobothnia Hospital District approved the protocol, and all the subjects gave their written informed consent. The primary outcomes of the study have been published recently (Junttila et al., 2018).
One-hundred sixty patients with angiographically documented CAD participated a 2-years home-based exercise training intervention. This ARTEMIS sub-study is previously described elsewhere in detail (Karjalainen et al., 2015). Forty-four patients were willing to measure R-R intervals (Polar RS800; Polar Electro Oy, Kempele, Finland) weekly basis at home during follow-up. More men than women were willing to measure R-R interval at home condition. During the 2-year follow-up, patients were asked to measure R-R intervals once or twice a week immediately after waking (before breakfast and washing up) using HR monitor. The home measurement protocol started with supine position (3 min) followed by sitting position (3 min), while the patients were advised to avoid additional movements and talking. At the baseline left systolic (ejection fraction, LVEF) and diastolic function (ratio of early transmitral flow velocity to early diastolic mitral annulus velocity, E/E’) were measured with two-dimensional tissue Doppler echocardiography (Vivid 7, GE Healthcare, Wauwatosa, WI, United States). Additionally, fasting blood samples were obtained at the baseline for analysis of inflammation, lipid, and glucose metabolism.
Analysis of R-R Interval Data
The R-R interval data were edited and analyzed in Hearts program (Heart Signal Co., Oulu, Finland). Ectopic beats and artifacts were removed from the tachogram based on visual inspection. Last 2 min for each phase was analyzed to secure stationarity of the data. Mean R-R intervals and the SD of normal-to-normal R-R intervals (SDNN) were used as time domain measures. Additionally, an autoregressive model (model 20) was performed to quantify the low- (LF, 0.04–0.15 Hz) and high-frequency (HF, 0.15–0.40 Hz) power of R-R interval oscillation from each phase.
In CAD patients with new CV event, measurements for 8 weeks preceding new CV event were included in the analyses. In patients without event, 8-weeks period corresponding to the median follow-up at the occurrence of new CV event in patients who encountered this event was included. The mean R-R interval and HR variability indexes were averaged over all recordings during the 8-weeks period within each patient. Also, the within-patient standard deviation (SD) of mean R-R interval and HR variability indexes were calculated. Finally, the mean R-R interval and HR variability indexes from the first and last two measurements of the 8-weeks period were averaged to evaluate possible trends occurring before the CV event.
New Cardiovascular Event
The participants were contacted through a mailed questionnaire and by telephone to inquire about the possible interim hospitalization after 2 years of follow-up. The final adjudication of the reason for hospitalization was ascertained from the primary discharge diagnosis from medical records. The deaths during the 2-year follow-up were registered using emergency rescue reports, hospital and physician records, autopsy data, death certificates, and interviews with the next of kin. The cause and mode of death were reviewed and adjudicated by two independent investigators; and if needed, disagreement or uncertainty was resolved in consultation with the investigators (MJJ and HH). The CV events prospectively included cardiovascular death, resuscitation from sudden cardiac arrest, hospitalization for either heart failure, acute coronary syndrome, stroke, or arrhythmia (atrial fibrillation/flutter or ventricular tachycardia) diagnosed according to the current guidelines.
Shapiro–Wilk’s test was used to examine the Gaussian distribution of the data. The variables that have non-Gaussian distribution were transformed into natural logarithm before the parametric statistical tests. Independent t-test was used for the between-group comparisons for continuous and chi-square test for categorical variables. Significant differences in variation of the mean R-R interval and HR variability between the groups were adjusted for E/E’, use of nitrates, use of alcohol, type 2 diabetes, sex, and sex × CV event interaction. based on differences in these baseline parameters between the groups with and without new CV event, as well as significant determinants of the variation in the mean R-R interval and HR variability. Pearson correlation analyses and independent t-test, followed by stepwise multiple linear regression, was used to assess the main determinants of variation in the mean R-R interval and HR variability. Two-way ANOVA for repeated measurement was used to assess possible differences in the mean R-R interval and HR variability (average of first vs. last 2 measurements) between the groups with and without CV event (time, group and time × group interaction) and was adjusted for E/E’, use of nitrates, use of alcohol, type 2 diabetes, and sex. The statistical analyses were performed using SPSS software, version 22.0 (SPSS Inc., Chicago, IL, United States). A p-value < 0.05 was considered statistically significant.
Eleven out of 44 patients had a new CV event (7 acute coronary syndromes, 1 sudden cardiac death, 2 new onset of arrhythmia needing hospitalization, and 1 stroke) during the follow-up (event time 11 ± 7 months after enrolment). Demographic characteristics of CAD patients with and without the new CV event are presented in Table 1. CAD patients with new CV event had higher E/E’ (13 ± 4 vs. 9 ± 3, p = 0.005) and they were more commonly treated with nitrates than their counterparts without the CV event (55 vs. 18%, p = 0.028). Other characteristics did not differ statistically significantly between the groups. However, patients with new event tended to have more often type 2 diabetes and more patients were alcohol users (Table 1). The mean number of R-R interval home measurements over 8 weeks was 12 ± 3 and 12 ± 3 for patients without and with new event, respectively (n = ns).
Table 1. Baseline characteristics of coronary artery disease patients with and without a new cardiovascular event.
Seated or supine mean R-R interval analyzed over 8 weeks did not differ between patients without and with new CV event (Table 2). The variability of mean R-R interval was greater in patients with new CV event compared to patient without the event when measured at supine (95 ± 34 vs. 59 ± 26 ms, p < 0.001, and p = 0.006 after adjustment) and sitting position (92 ± 28 vs. 62 ± 24 ms, p = 0.001, and p < 0.001 after adjustment) (Table 2). Representative examples of the variability of mean R-R interval during 8 weeks from patients with and without CV event are shown in Figure 1.
Table 2. Heart rate and heart rate variability in CAD patients with and without new cardiovascular event without and with adjusted for left ventricular diastolic function and use of nitrates, alcohol, type 2 diabetes, and sex.
Figure 1. Representative examples of the variability of mean R-R intervals in supine and sitting positions for 8 weeks preceding cardiovascular event (A,C) and without the event (B,D).
Any mean value of HR variability measures analyzed over 8 weeks did not differ between groups (Table 2). However, the variability of HF power over 8 weeks was greater in patients with new CV event when measured in supine (0.74 ± 0.23 vs. 0.56 ± 0.23 ms2, p = 0.012, and p = 0.003 after adjustment) and sitting positions (0.90 ± 0.35 vs. 0.61 ± 0.26 ms2, p = 0.007, and p = 0.038 after adjustment). On the contrary, the variability of LF power while seated was smaller in patients with new CV event (0.70 ± 0.27 vs. 0.83 ± 0.40 ms2, p = 0.009, and p = 0.010 after adjustment). The determinants of the variation in the mean R-R interval and HR variability are presented in Table 3. No significant sex × CV event interaction was observed in any HR variable.
Table 3. Determinants of measurement-to-measurement variation in the mean R-R interval and heart rate variability by stepwise linear regression.
When comparing the changes in the mean R-R interval and HR variability from the first two and last two measurements during 8-weeks period, significant interaction was observed in seated and supine mean R-R interval, which increased more in the patients with new CV event compared to patients without new CV event (Table 4). Also, sitting SDNN and HF power increased more in the patients with new CV event than those without.
Table 4. Heart rate and heart rate variability in CAD patients with and without a new cardiovascular event as average of the two first and two last measurements of the 8-weeks period without and with adjusted for left ventricular diastolic function, use of nitrates, alcohol, type 2 diabetes, and sex.
The main finding of the present study was that day-to-day variability of mean R-R interval was greater before new CV event in patients with initially stable CAD. The analysis of day-to-day variation of HR variability measures did not provide additional information to the variation of mean R-R interval. The present results suggest that a more unstable autonomic balance precedes the new CV events. Therefore, home monitoring of resting HR may become a simple and useful method for providing prognostic information on the risk of imminent CV event in patients with stable CAD.
The risk of CV events related to increased day-to-day variability blood pressure and HR have been well documented in general populations (Kikuya et al., 2008; Johansson et al., 2012). Temporal changes HR and HR variability before imminent CV event are evidently less studied. Adamson et al. (2004) demonstrated that HR variability decreased, and the mean HR increased before hospitalization among patients with heart failure. However, daily variation in HR variability or mean HR were not reported. Furthermore, home telemonitoring, using, e.g., HR, blood pressure, and weight control, have shown to reduce hospitalization in heart failure patients (Koulaouzidis et al., 2016; Sardu et al., 2016; Tse et al., 2018). In the present study, no significant differences were found in HR variability or mean R-R interval between the patients with and without new CV event, when averaged over 8 weeks preceding the CV event or corresponding period in event-free patients. Instead, we observed that day-to-day variation in the mean R-R interval and HF power of HR variability was greater among patients who encountered the new CV event compared to those who did not. One possible explanation for the present finding is fluctuating cardiac autonomic regulation with progressing myocardial ischemia, which was underlying most of the new CV events in the present study. The triggers, such as physical, emotional and other stressors, and, thus, symptoms of angina pectoris may vary day-to-day which may well explain the present findings (Task Force Members et al., 2013). Notably, higher variation of the mean R-R interval was partly explained by the history of coronary artery bypass grafting and the use of nitrates indicating a more advanced CAD. We do not have information about patients’ compliance and symptom-based need for nitrate medication that may possibly contribute to daily variation of angina pectoris and cardiac autonomic regulation. However, the difference in the variation of the mean R-R interval and HR variability between patients with and without new CV event remained significant even after adjustments for these confounders.
We observed increase in SDNN as well as in the mean R-R interval and HF power during the 8-weeks period preceding the new CV event compared to event-free control group. This is contradictory to our expectations since acute decrease in cardiac vagal activity has been typically observed at the onset of myocardial ischemia (Sroka et al., 1997). The present findings also oppose the previous results regarding hospitalization due to heart failure by Adamson et al. (2004) one possible explanation for these divergent findings is that assessed progression of HR variability before hospitalization due to heart failure, whereas the present study assessed HR dynamics before the new CV event that were mainly related to progression new coronary events. While increased sympathetic activity and decreased vagal activity with progressing heart failure is well documented (Nolan et al., 1992; van de Borne et al., 1997), myocardial ischemia may produce varying autonomic responses. For example, Lombardi et al. (1996) showed that autonomic response to acute coronary syndrome depends on the localization of myocardial ischemia, i.e., anterior wall ischemia may induce sympathetic activation and inferior wall ischemia vagal activation. Therefore, decreased HR variability with progressing CAD may not be constantly observed. Secondly, a co-activation of vagal and sympathetic outflow may result in either increased HR variability due to the sudden and large beat-to-beat R-R interval changes or decreased HR variability due to the sudden onset of fixed R-R interval dynamics (Tulppo et al., 1998; Tulppo et al., 2005). The episodes of co-activation of both autonomic arms during short-term R-R interval recordings, such as in the present study, may partly explain our inconsistent HR variability dynamics before CV events.
The present study may have important practical implications. It seems that monitoring of the variation of mean HR may suffice to evaluate risk for imminent CV event because analysis of HR variability did not provide additional information. In addition, HR measurements at seated position appeared to be related to the occurrence of CV event more than supine measurements. These findings suggest that HR monitoring may be easily accompanied with home monitoring of BP. Some caution may be needed since HR, acquired during BP measurement, may not be as robust as mean HR obtained from R-R interval measurements. Finally, it must be acknowledged that HR dynamics may complement a symptom-based follow-up of CAD patients as it may also capture autonomic alterations also related to silent myocardial ischemia (Goseki et al., 1994).
The major limitation of the present study is the lack of ECG data measured at the home condition. Unfortunately, at the time of our data collection (starting 2007) we did not have suitable date recorders available to collect ECG at home. Remote technologies to measure not only ECG but also respiratory frequency, temperature, physical activity, etc., are developing very rapidly. Today, ECG and other signals can be measured at home and send data easily to the health care centers via internet (Fung et al., 2015). The present study is also limited by the small number of patients having CV event during the follow-up. Therefore, more research will be needed in a larger sample size to confirm the definite value of day-to-day variation of HR and HR variability to provide additional information of CV risk among CAD patients. Also, it must be considered that the measurements at home were short and could not be strictly controlled for potential confounders. However, they were performed in real-life conditions that may have greater value for practical implications. Additionally, stationary data can be obtained with very short stabilization period (Flatt and Esco, 2016).
Day-to-day variability of mean R-R interval was greater before new CV event in CAD patients suggesting to a more unstable autonomic balance preceding this event. If confirmed in a large sample size, home HR measurement may become a simple and useful method for providing information on the risk of imminent CV event in patients with stable CAD.
All datasets generated for this study are included in the manuscript and/or the supplementary files.
MT, AK, and HH contributed conception and design of the study. MT and AK organized the database and performed the statistical analysis. MT wrote the first draft of the manuscript. MT, AK, MJJ, and HH wrote the sections of the manuscript. All authors contributed to manuscript revised, read, and approved the submitted version.
This work was funded by the Finnish Technology Development Center (TEKES) Helsinki, Finland, the Academy of Finland (#267435), Helsinki, Finland, and the Sigrid Jusélius Foundation, Helsinki, Finland.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors gratefully acknowledge the assistance of the registered nurses Pirkko Huikuri, Päivi Koski, Päivi Kastell, and Sari Kaarlenkaski.
Adamson, P. B., Smith, A. L., Abraham, W. T., Kleckner, K. J., Stadler, R. W., Shih, A., et al. (2004). Continuous autonomic assessment in patients with symptomatic heart failure: prognostic value of heart rate variability measured by an implanted cardiac resynchronization device. Circulation 110, 2389–2394. doi: 10.1161/01.CIR.0000139841.42454.78
Airaksinen, K. E., Ikaheimo, M. J., Linnaluoto, M. K., Niemela, M., and Takkunen, J. T. (1987). Impaired vagal heart rate control in coronary artery disease. Br. Heart J. 58, 592–597. doi: 10.1136/hrt.58.6.592
Bigger, J. T. Jr., Fleiss, J. L., Steinman, R. C., Rolnitzky, L. M., Schneider, W. J., Stein, P. K., et al. (1995). RR variability in healthy, middle-aged persons compared with patients with chronic coronary heart disease or recent acute myocardial infarction. Circulation 917, 1936–1943. doi: 10.1161/01.CIR.91.7.1936
Chowdhury, E. K., Owen, A., Krum, H., Wing, L. M., Nelson, M. R., Reid, C. M., et al. (2014). Systolic blood pressure variability is an important predictor of cardiovascular outcomes in elderly hypertensive patients. J. Hypertens. 32, 525–533. doi: 10.1097/HJH.0000000000000028
Fung, E., Jarvelin, M. R., Doshi, R. N., Shinbane, J. S., Carlson, S. K., Grazette, L. P., et al. (2015). Electrocardiographic patch devices and contemporary wireless cardiac monitoring. Front. Physiol. 6:149. doi: 10.3389/fphys.2015.00149
Goseki, Y., Matsubara, T., Takahashi, N., Takeuchi, T., and Ibukiyama, C. (1994). Heart rate variability before the occurrence of silent myocardial ischemia during ambulatory monitoring. Am. J. Cardiol. 73, 845–849. doi: 10.1016/0002-9149(94)90807-9
Huikuri, H. V., Niemelä, M. J., Ojala, S., Rantala, A., Ikaheimo, M. J., Airaksinen, K. E., et al. (1994). Circadian rhythms of frequency domain measures of heart rate variability in healthy subjects and patients with coronary artery disease. Effects of arousal and upright posture. Circulation 90, 121–126. doi: 10.1161/01.CIR.90.1.121
Johansson, J. K., Niiranen, T. J., Puukka, P. J., and Jula, A. M. (2012). Prognostic value of the variability in home-measured blood pressure and heart rate: the finn-home study. Hypertension 59, 212–218. doi: 10.1161/HYPERTENSIONAHA.111.178657
Junttila, M. J., Kiviniemi, A. M., Lepojarvi, E. S., Tulppo, M., Piira, O. P., Kentta, T., et al. (2018). Type 2 diabetes and coronary artery disease: preserved ejection fraction and sudden cardiac death. Heart Rhythm 15, 1450–1456. doi: 10.1016/j.hrthm.2018.06.017
Karjalainen, J. J., Kiviniemi, A. M., Hautala, A. J., Piira, O. P., Lepojarvi, E. S., Perkiomaki, J. S., et al. (2015). Effects of physical activity and exercise training on cardiovascular risk in coronary artery disease patients with and without type 2 diabetes. Diabetes Care 38, 706–715. doi: 10.2337/dc14-2216
Kikuya, M., Ohkubo, T., Metoki, H., Asayama, K., Hara, A., Obara, T., et al. (2008). Day-by-day variability of blood pressure and heart rate at home as a novel predictor of prognosis: the ohasama study. Hypertension 52, 1045–1050. doi: 10.1161/HYPERTENSIONAHA.107.104620
Lombardi, F., Sandrone, G., Spinnler, M. T., Torzillo, D., Lavezzaro, G. C., Brusca, A., et al. (1996). Heart rate variability in the early hours of an acute myocardial infarction. Am. J. Cardiol. 77, 1037–1044. doi: 10.1016/S0002-9149(96)00127-0
Nolan, J., Flapan, A. D., Capewell, S., MacDonald, T. M., Neilson, J. M., Ewing, D. J., et al. (1992). Decreased cardiac parasympathetic activity in chronic heart failure and its relation to left ventricular function. Br. Heart J. 67, 482–485. doi: 10.1136/hrt.67.6.482
Sardu, C., Santamaria, M., Rizzo, M. R., Barbieri, M., di Marino, M., Paolisso, G., et al. (2016). Telemonitoring in heart failure patients treated by cardiac resynchronisation therapy with defibrillator (CRT-D): the TELECART Study. Int. J. Clin. Pract. 70, 569–576. doi: 10.1111/ijcp.12823
Task Force Members, Montalescot, G., Sechtem, U., Achenbach, S., Andreotti, F., Arden, C., et al. (2013). 2013 ESC guidelines on the management of stable coronary artery disease: the Task Force on the management of stable coronary artery disease of the European Society of Cardiology. Eur. Heart J. 34, 2949–3003. doi: 10.1093/eurheartj/eht296
Tse, G., Chan, C., Gong, M., Meng, L., Zhang, J., Su, X. L., et al. (2018). Telemonitoring and hemodynamic monitoring to reduce hospitalization rates in heart failure: a systematic review and meta-analysis of randomized controlled trials and real-world studies. J. Geriatr. Cardiol. 15, 298–309.
Tulppo, M. P., Kiviniemi, A. M., Hautala, A. J., Kallio, M., Seppänen, T., Mäkikallio, T. H., et al. (2005). Physiological background of the loss of fractal heart rate dynamics. Circulation 112, 314–319. doi: 10.1161/CIRCULATIONAHA.104.523712
Tulppo, M. P., Mäkikallio, T. H., Seppänen, T., Airaksinen, J. K., and Huikuri, H. V. (1998). Heart rate dynamics during accentuated sympathovagal interaction. Am. J. Physiol. 274(3 Pt 2), H810–H816. doi: 10.1152/ajpheart.1998.274.3.H810
van de Borne, P., Montano, N., Pagani, M., Oren, R., and Somers, V. K. (1997). Absence of low-frequency variability of sympathetic nerve activity in severe heart failure. Circulation 95, 1449–1454. doi: 10.1161/01.CIR.95.6.1449
Webb, A. J. S., Mazzucco, S., Li, L., and Rothwell, P. M. (2018). Prognostic Significance of blood pressure variability on beat-to-beat monitoring after transient ischemic attack and stroke. Stroke 49, 62–67. doi: 10.1161/STROKEAHA.117.019107
Yu, J. M., Kong, Q. Y., Schoenhagen, P., Shen, T., He, Y. S., Wang, J. W., et al. (2014). The prognostic value of long-term visit-to-visit blood pressure variability on stroke in real-world practice: a dynamic cohort study in a large representative sample of Chinese hypertensive population. Int. J. Cardiol. 177, 995–1000. doi: 10.1016/j.ijcard.2014.09.149
Keywords: heart rate, coronary artery disease, cardiovascular event, home monitoring, heart rate variability
Citation: Tulppo MP, Kiviniemi AM, Junttila MJ and Huikuri HV (2019) Home Monitoring of Heart Rate as a Predictor of Imminent Cardiovascular Events. Front. Physiol. 10:341. doi: 10.3389/fphys.2019.00341
Received: 18 December 2018; Accepted: 13 March 2019;
Published: 27 March 2019.
Edited by:Gaetano Santulli, Columbia University, United States
Reviewed by:James Hodgkinson, University of Birmingham, United Kingdom
Amalia Peix, Instituto de Cardiología y Cirugía Cardiovascular, Cuba
Copyright © 2019 Tulppo, Kiviniemi, Junttila and Huikuri. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.