Autonomic Dysfunction and Risk Stratification Assessed from Heart Rate Pattern
A Günther*, O.W Witte, D Hoyer
Identifiers and Pagination:Year: 2010
First Page: 39
Last Page: 49
Publisher Id: TONEUJ-4-39
Article History:Received Date: 11/5/2009
Revision Received Date: 22/12/2009
Acceptance Date: 4/2/2010
Electronic publication date: 15/6/2010
Collection year: 2010
open-access license: This is an open access article licensed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted, non-commercial use, distribution and reproduction in any medium, provided the work is properly cited.
The modulation of the autonomic nervous system (ANS) under physiological and pathophysiological conditions is in focus of recent research. Many patients with cardio- and cerebrovascular diseases display features of sympathovagal dysregulation. Measuring specific ANS parameters could improve risk stratification. Thus, the early diagnosis of ANS dysfunction in these patients poses a great challenge with high prognostic relevance.
The most relevant methods and measures of Heart Rate Variability (HRV) analysis and HRV monitoring will be described in detail in this chapter. The grown importance of these easily obtainable heart rate patterns in stratifying the risk of patients with myocardial infarction and heart failure as well as ischemic stroke will be demonstrated based on recent clinical studies. In order to perspectively improve clinical management of these patients further large scale clinical investigations on the role of ANS dysfunction will be useful.
Recent research has made important progress in characterizing the role of the autonomic nervous system (ANS) under physiological and pathophysiological conditions as an integrative component in the interorgan interplay.
Early diagnosis of ANS dysfunction within intensive care patients poses a challenge as the majority will go on to display features of sympatho-vagal dysregulation during the course of the disease (e.g. ).
Through the measurement of specific ANS parameters clinicians and researchers are enabled to improve cardiovascular and cerebrovascular risk stratification.
The following chapter briefly outlines ANS functions and interactions and summarizes the most relevant methods of its measurement through Heart Rate Variability (HRV) analysis. Prognostic evaluation of critical care patients using such a cost and time-efficient tool as a Holter-ECG in conjunction with clinical scores have been assessed in a number of clinical studies. Here, we will review the present evidence from these studies for such highly relevant pathophysiological conditions of heart failure and myocardial infarction with focus on prognosis of sudden cardiac death, and stroke with a special emphasis neurological outcome related to ANS dysfunction.
2. THE AUTONOMIC NERVOUS SYSTEM (ANS)
Autonomic nervous system (ANS) responses to (patho-) physiological stimuli have been characterized by a complex interaction between the different components of the ANS on the one hand, and visceral organs and (involuntary) innervated glands on the other hand. In conjunction with the sensory and motor nervous system, the autonomic nervous system (ANS) is responsible for fast, continuous and normally subconscious control of basic physiological functions, such as heart rate, blood pressure, respiratory rate, body temperature and gastrointestinal motility. Forming a functional network of peripheral and central mechanisms throughout the body the autonomic nervous system, however, has two basic divisions: the parasympathetic and the sympathetic component. The classical thinking that parasympathetic and sympathetic nervous system are usually acting in opposition holds true for many ANS functions as in cardiovascular, respiratory as well as gastrointestinal and urinary tract innervations and reactions. But in many other situations the two divisions work synergistically in response to a simultaneous sympathetic and vagus nerve stimulus.
The ANS consists of efferent neurons running in series to innervate smooth musculature of all viscerals and vessels, the heart and glands. Afferent neurons relay visceral input to different parts of the brain. The sympathetic postganglionic transmitter is norepinephrine, whereas the parasympathetic postganglionic transmitter is acetylcholine.
ANS, however, also interacts with different brain structures. The central and peripheral autonomic information is processed by the central nervous system, involving forebrain structures (frontal premotor cortex), hypothalamus, pituitary, amygdala, and the insular cortex. Included in this central network are also several brainstem structures, such as nucleus tractus solitarius (NTS), dorsal motor nucleus of the vagus (DMN), and nucleus ambiguous (NA).
To generate a homeostatic signal level three pathways are relevant: the hypothalamo-pituitary-adrenal (HPA) axis, the parasympathetic and the sympathetic nervous system . Hypothalamic signal input results in the release of corticotrophin-releasing factor from the paraventricular nucleus, by which the HPA axis and the sympathetic nervous system as well as the parasympathetic system are modulated. The HPA axis activation finally leads to production of glucocorticoids, resulting in immunomodulation among other hormone related effects.
Peripherally, afferent vagus nerve fibers forward visceral (inflammatory) signals to the brain. Inhibition of inflammation by suppressing pro-inflammatory cytokine release through efferent vagus nerve signals (mediated by acetylcholine, the principle vagal neurotransmitter) has been proposed as the cholinergic anti-inflammatory pathway . Synergistic with the parasympathetic nervous system, the sympathetic activation and secretion of epinephrine and norepinephrine from the adrenal medulla and release from sympathetic nerve terminals also results in beta adrenergic mediated immunological and autonomic modulation [4, 5].
These mechanisms aim to keep sympatho-vagal homeostasis. The role of the ANS, however, remains to be further elucidated as a possible key player in the bodies’ response to pathophysiological conditions in order to consolidate cardio- and cerebrovascular and immune function as well as to develop implications for early diagnosis and treatment of respective dysfunctions.
3. METHODS OF HRV ANALYSIS
3.a. Introduction of Linear and Nonlinear Methods
The fluctuations of heart rate essentially reflect modulations mediated by the autonomic nervous system. They provide insights into different parts of autonomic control as well as risk indices in various pathological entities. In several studies traditional time and frequency domain measures of heart rate variability (HRV) have been investigated in conjunction with cardiovascular outcomes as an important compartment of risk stratification. Furthermore, non-linear HRV measures, such as fractal scaling exponents, complexity indices and indices based on Poincarè plots have been applied in this context.
All of these measures are based on electrocardiographic (ECG) recordings. By means of mobile so called “Holter ECG” systems, the patient’s HRV can be assessed under “normal” daily activity. In clinical settings the ECG can be recorded both by means of Holter systems and by stationary bedside monitoring systems. In order to obtain a reasonable precision of the heart rate fluctuations the ECG should be sampled with rate of at least 1000 Hz. The heart beats have to be exactly identified for example by template matching or a steepest ascent criterion followed by an annotation of normal and arrhythmic beats. In order to focus on the autonomic rhythms and discriminate from those beats reflecting cardiogenic arrhythmias the convention of analyzing only intervals between normal heart beats (so called NN intervals) was set more than a decade ago . Clean NN-interval series can be constituted after treating arrhythmias by appropriate interpolation or deletion. The relative number of such removed arrhythmic beats is an important measure of the data quality and can itself provide additional prognostic value beyond the “Task Force HRV” analysis .
Indices of HRV Under Spontaneous Conditions
HRV indices can be calculated from the overall 24 hour recording or from shorter periods which reflect particular situations. Also the circadian rhythm of autonomic activity should be taken into consideration.
Linear HRV assessment of those type of data was first organized by a Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology in 1996 , and has recently been successfully validated in several clinical studies as demonstrated in the studies referred below. These guidelines are still the established HRV standard, even though an update has been recommended [1, 7]. This mainly “linear” HRV analysis is outlined in more detail below (see 3b).
Nonlinear HRV assessment was proposed by using different mathematical approaches:
Fractal approaches such as the beta-slope in the log-log power spectrum and indices of detrended fluctuation analysis (DFA, [8, 9] assess the self-similar qualities of the HR time series between different time scales. Their success as predictors of outcome supports the potential relevance of measuring complex cardiovascular coordination over a wide range of scales and frequencies, respectively. Normal behaviour is associated with fractal-like correlated properties. Pathophysiological conditions can change the fractal scaling exponent or even distort the self-similarity in favour of single scale rhythms at one extreme or random noise at the other.
Complexity approaches assess the randomness (unpredictability) of the NN interval series. This can be done based on the original data leading to approximate entropy , sample entropy (, conditional entropy , or information loss [13-15]
Those approaches are furthermore a fundamental way of analyzing the information flow that characterizes communication. Cardiovascular control operates by communication between different parts of the autonomic nervous system at different time scales [13,15] Complexity measures as mentioned above are calculated using one selected time scale [10,18] Costa et al.  introduced “Multiscale entropy”, a complexity function over multiple scales, and demonstrated the relevance of heart rate complexity pattern on different time scales.
Autonomic Information Flow (AIF) functions reflect the information flow over different prediction time horizons and time scales, respectively, of NN interval series. AIF functions allow the assessment of complexity of both short and long term heart rate patterns at scales analogous to those of traditional HRV measures . In that way the complex rhythms and interactions, essentially mediated by the autonomic nervous system were introduced into cardiovascular risk assessment. The capability of AIF analysis opens a new window into the development of physiologically relevant complexity measures by considering underlying autonomic mechanisms and the extensive knowledge base already gained using the standard HRV approaches of task force monitoring .
These methods focus on autonomic rhythms, reflected by normal-to-normal (NN) beat intervals after removing and interpolating cardiogenic arrhythmic parts from the measured heart beat interval series. On the other hand, cardiogenic arrhythmias such as ectopic beats are of particular diagnostic and prognostic value using different mathematical tools.
Heart Rate Turbulence (HRT) is an alternative method where the heart rate response to premature ventricular complexes (VPC) is evaluated . This response reflects autonomic blood pressure regulation associated activity induced by a decreased heart ejection volume period. Reduced HRT is significantly associated with mortality after myocardial infarction .
Heart Rate Responses Under Test Conditions
Orthostatic tilt test, baroreflex sensitivity by phenylephrine administration, chemoreflex sensitivity by breathing air with increased or decreased oxygen, cardiovascular effects of other drugs, and controlled ventilation are prominent examples of selectively challenging the autonomic system and consequently stimulation of well defined heart rate responses [24,25]. Those tests address particular control loops more specifically than the indices of HRV under resting conditions. Otherwise they are laborious and do not allow a long time monitoring.
The spontaneous and instantaneous heart rate response to blood pressure fluctuations is an index of “spontaneous baroreflex sensitivity”. The so-called “alpha-index” reflects the coupling within particular frequency bands of heart rate and blood pressure rhythms [26,27]. The “sequence technique” reflects the coupling between sequences of increasing or decreasing, respectively, heart rate and blood pressure values [26,28].
3.b. “Task Force” HRV Monitoring
Besides the heart rate fluctuations assessed as HRV measures already the mean heart rate provides valuable information about the autonomic control. The main time and frequency domain HRV measures according to  will be introduced as followed.
Selected Prominent Measures of Linear HRV Analysis
In the time domain different measures of the dispersion of the NN intervals are assessed from the entire data set:
SDNN: Standard deviation of all normal to normal (NN) intervals. This measure reflects the overall variability whereby the slower rhythms contribute stronger due to their larger amplitudes.
RMSSD: Square root of the mean squared differences of successive NN intervals. This measure reflects faster rhythms which are usually assigned to vagal rhythms.
Further related indices are:
SDANN: Standard deviation of the average of NN intervals in all 5 min segments of the entire recording.
SDNNindex: Mean of the standard deviations of all NN intervals for all 5 min segments.
HRVtriangular index: Total Number of all NN intervals divided by the height of the histogram of all NN intervals.
TINN: Baseline width of the maximum square difference triangular interpolation of the highest peak of the histogram of all NN intervals.
pNN50: Number of pairs of adjacent NN intervals differing by more than 50 ms divided by the total number of NN intervals in the entire recording.
In the frequency domain the power spectral density reflects the distribution of independent sine waves. Usually the NN interval series is equidistantly re-sampled and power spectra are calculated from the resulting time series. This can be done either by a Fast Fourier Transform algorithm using appropriate windowing or by statistically reliable autoregressive models. Different components of the power spectrum (frequency band power) are assigned corresponding to typical rhythms of the autonomic modulations of the heart period:
The short term rhythms are accentuated by using the mean of the measures obtained by analyzing short windows (5min: HF, 5min: LF, 30min: VLF) in which the entire data set is subdivided:
TP (Total Power): lowest signal frequency-0.4 Hz (overall HRV).
LF (Low Frequency): 0.04-0.15 Hz, reflecting vagal and sympathetical mediated control.
LFnorm: LF/(TP-VLF): emphasizes the relative contribution of LF.
HF (High Frequency): 0.15-0.4 Hz, predominantly reflecting vagal mediated control.
HFnorm: HF/(TP-VLF): emphasizes the relative contribution of HF.
LF/HF: ratio interpreted as sympatho-vagal balance.
In Fig. (1) the change of a short term power spectrum due to tilt load is shown. The shift of the sympatho-vagal balance is reflected in the relative increase of LF and decrease of HF band power.
Sketch of short term HRV frequency domain characteristics according to : Power spectral density during rest (dotted line) and tilt load (dashed line). Note the shift of the sympathovagal balance represented by the LF and HF peaks.
The long term rhythms are analyzed from the entire data set as a whole. Besides frequency bands corresponding to the short term measures the Ultra Low Frequency (ULF) band is additionally considered. The resulting 24 hours power spectrum mainly reflects the circadian rhythm. The frequency bands of these long term power spectra were found to be predominant predictive. The square root of the total power (TP) is equivalent to SDNN of the same data set. Moreover, the decay of a double logarithmic plot of this power spectrum (beta-slope), which quantifies its fractal behaviour, turned out as important measure of complex multi-scale cardiovascular behaviour.
Beta slope: Slope of a line fitted to the log-log-power spectrum between 10-4 and 10-2 Hz (Fig. 2).
Sketch of a long term power spectrum. Note the slope of the fitted line between 10-4 and 10-2 Hz, respectively -4 and -2 in the logarithmic scale.
Each of the indices in time and frequency domain has different advantages, and robustness when applied to measured clinical data. Hence the most appropriate of redundant indices (correlations see Table 1) have to be chosen.
Fundamental Correlations Between Time Domain Measures (Grouped Left) and Frequency Domain Measures (from [6,30])
|Time Domain||Frequency Domain|
|SDNN, SDANN, HRVtri.ind., TINN||x|
The classical physiological interpretation of the mediators underlying the heart rate fluctuations reflects consensus on a simplified level, but also taking into account complex inter-relationships between multiple autonomic, humoral and other processes. Under particular pathological conditions these relationships can even be changed . As a basic rule of dynamic processes can be applied that faster processes influence slower ones, such as discussed for vagal influences reflected in VLF band ). In critical care risk stratification parasympathetic indices (RMSSD and HFshort term) as well as VLFshort term provided prognostic value. Some respective implications such as those related to autonomic regulations of inflammation, the activity of the renin-angiotensin system, or respiratory dysfunctions are discussed in [1,32]. Furthermore, reduced SDNN which reflects the overall HRV and hence the effect of multiple different physiological controllers and influences up to the circadian rhythm turned out as risk predictor in several studies, a result that confirms the comprehensive systemic approach of analyzing autonomic rhythms by HRV analysis. The same holds true for Beta-slope as a comprehensive long term HRV measure.
The general procedure of (linear) HRV analysis is outlined in Fig. (3) according to . The later developed (non-linear) complexity and AIF indices are not yet organized following a respective standard scheme. Nevertheless, they are usually assessed following related considerations, however, with focus at complex interactions within autonomic control.
Flow chart of the established standard HRV analysis according to . Data sets with more than 5% arrhythmic beats excluded from the analysis. This scheme holds also true for novel complexity and AIF indices.
Even though easily applicable, the method of HRV analysis absolutely relies on good signal quality and its major limitation is the preclusion of atrial fibrillation, sinoatrial dysfunction and >20% ectopic complexes.
4. HEART FAILURE AND MYOCARDIAL INFARCTION
Myocardial infarction (MI) and heart failure are the leading causes of death and morbidity in industrial countries. Investigations of patients after MI revealed ANS dysfunction as a predictor for cardiac mortality in association with increased incidence of sudden cardiac death . For the general clinical use HRV measures are of particular interest in cases where they add prognostic value independent of the established clinical variables or where they can substitute complicated clinical scores by easily obtainable HRV indices during monitoring.
Therefore, since more than one decade HRV analysis is established in clinical cardiology and accepted as a part of risk assessment concerning sudden cardiac death. Currently, the development of HRV based ICD (implantable cardioverter defibrillator) control parameters is one of the predominant questions.
The autonomic rhythms represented by heart rate patterns reflect certain pathophysiological aspects of several cardiac disorders, like myocardial infarction, congestive heart failure, cardiomyopathy, or ischemic heart disease. The comorbidity of cardiac diseases with for example psychiatric or mental disorders leads to even more complicated identification of the origin of a measured HRV pattern. For example, the very low frequency (VLF) band power in depressive in comparison to non-depressive patients (both groups with recent acute myocardial infarction (MI)) is reduced due to depression and MI, whereby both disorders constitute independent factors of mortality .
Representative results of larger clinical prognosis studies of patients with cardiac diseases using 24 hours HRV are summarized in Table 2. A comprehensive review of HRV based predictors of mortality after myocardial infarction by Kleiger et al.  includes the selected key results presented here. The strongest clinical risk predictor is a reduced left ventricular ejection fraction (LVEF). A main objective of clinical studies is to evaluate whether HRV indices provide independent additional prognostic value.
Studies of Heart Rate Related ANS Dysfunction in Heart Failure and Myocardial Infarction Patients
|Prognostic Heart Failure Trials||Data Sample Set||Results|
|Prospective study of heart rate variability and mortality in chronic heart failure: results of the United Kingdom heart failure evaluation and assessment of risk trial (UK-heart) ||482+/-161 days follow-up.
433 patients, 54 deaths.
24 hour Holter records .
|Reduced SDNN is an independent predictor besides cardiothoracic ratio, left ventricular end-systolic diameter, and serum sodium.|
|Depressed low frequency power of heart rate variability as an independent predictor of sudden death in chronic heart failure ||22+/-18 month follow-up.
190 patients, 55 deaths.
24 hour / 10-19h daytime Holter records.
|Reduced SDNN is an independent predictor besides coronary artery disease and cardiothoracic ratio.
Reduced LF during daytime is an independent predictor besides coronary artery disease in a sudden death subgroup.
|Short-term heart rate variability strongly predicts sudden cardiac death in chronic heart failure patients ||Three years follow up.
202 patients (total mortality: 37%, sudden death mortality: 9.4%) in derivation sample.
242 patients (total mortality: 22%, sudden death mortality: 8%) in the validation sample.
5 min data sets during controlled breathing.
|In the derivation sample reduced LF was an independent sudden death predictor besides left ventricular end-diastolic diameter.
In the validation sample reduced LF was an independent sudden death predictor besides ventricular premature contractions.
|Different spectral components of 24 h heart rate variability are related to different modes of death in chronic heart failure ||Three years follow up.
330 patients, mortality due to pump failure / urgent transplantation: 24%, due to sudden death: 9%.
24 hour, night-time 0-5h, daytime 9-19h Holter records.
|Reduced VLF (night-time), high pulmonary wedge pressure and low left ventricular ejection fraction independently predic death for progressive pump failure / urgent transplantation.
Reduced LF (night-time) and increased left ventricular end-systolic diameter independently predict sudden death.
|Nonlinear indices of heart rate variability in chronic heart failure patients: redundancy and comparative clinical value ||31 (4-36) months follow up.
200 patients, cardiac death n=60, heart transplantation n=37, non-cardiac death n=2, drop out n=2.
24 hour Holter records.
|Only 1VP (symbolic pattern) and IMAI2 (LF-HF related empirical mode decomposition) added independent prognostic value to major clinical and functional predictors: symptom severity, left ventricular ejection fraction, peak VO2 at cardiopulmonary exercise testing and systolic arterial pressure, however other candidates such as sample entropy, DFA, and alpha slope did not.|
|Autonomic response to cardiac dysfunction in chronic heart failure - a risk predictor based on autonomic information flow ||31 (4-36) months follow up.
200 patients, 100 cardiac death or heart transplantation.
24 hour Holter records.
|BDNN and PDmVLF added independent prognostic value to NYHA class, LVEF, SAP, peak_VO2 (see example below).|
|Prognostic Post Myocardial Infarction Trails|
|Baroreflex sensitivity and heart-rate variability in prediction of total cardiac mortality after myocardial infarction. ATRAMI (Autonomic Tone and Reflexes After Myocardial Infarction) Investigators ||21 (8) months follow up
1284 patients, 44 cardiac deaths and 5 non-fatal cardiac arrests
24 hour Holter records
|Low values of either BRS or SDNN are independent predictors besides left ventricular ejection fraction and the rate of ventricular premature contractions. The association of low BRS and SDNN further increases risk.|
|Fractal analysis of heart rate dynamics as a predictor of mortality in patients with depressed left ventricular function after acute myocardial infarction. TRACE Investigators (TRAndolapril Cardiac Evaluation) ||4-years follow-up
159 patients with depressed left ventricular ejection fraction after myocardial infarction, 72 patients died.
24 hour Holter records
|DFA based scaling exponent (<11 beats) was independent predictor in a multivariable Cox proportional regression model besides age, history of angina pectoris, diuretic medication, and left ventricular ejection fraction.
None of the linear HRV indices had prognostic value.
|Mortality in patients with unstable angina|
|Prognostic value of ventricular arrhythmias and heart rate variability in patients with unstable angina ||6 months follow up
543 patients, 40 death
24 hour Holter records
|Ventricular arrhythmias add independent predictive value to clinical and laboratory variables.
SDNNindex, LF, LF/HF are predictive in univariate analysis. But LF only adds independent predictive value in multivariate Cox survival analysis.
Reflecting the briefly summarized results presented in Table 2 obviously the reduced overall variability (SDNN, HRV index) or slow rhythms (VLF band power) provide consistent HRV based risk predictors.
In chronic heart failure patients and those with unstable angina also reduced LF or LF/HF showed predictive value. A reduced DFA index (reduced fractal correlations) was found of superior predictive value than SDNN in a study of post-myocardial infarction patients .
Due to the interindividual heterogeneity and the complex contributions such as autonomic, mental, psychiatric, physical activity, circadian rhythm, etc. besides the cardiac dysfunction, and furthermore the often unclear causality and comorbidity, it is unlikely to find the “optimal” HRV index. However, selected indices may be sensitive to particular aspects of pathophysiological autonomic control and hence particular patient groups.
Maestri et al.  compared 20 nonlinear indices and SDNN as the best discriminating linear HRV index in a group of 200 patients with chronic heart failure (CHF). They found two nonlinear indices (one from empirical mode decomposition, one from symbolic dynamics) as significant predictors independent of the clinical parameters and SDNN. A third respective set of significant nonlinear predictors of this CHF patient group was found by autonomic information flow (AIF) characteristics . These provided additional predictive value to the best clinical indices NYHA class and LVEF. The Kaplan Mayer survival curves (Fig. 4) show their discriminatory power. The result of this patient group indicates the predictive value of the complex overall autonomic adjustment as superior compared to particular rhythms, such as HF or LF, which did not significantly contribute in that study.
Kaplan-Meier survival curves for total cardiac mortality, AIF index according to .
In a scientific statement from the American Heart Association Council on Clinical Cardiology Committee on Electrocardiography and Arrhythmias and Council on Epidemiology and Prevention noninvasive risk stratification techniques for identifying patients at risk for Sudden Cardiac Death (SCD) were reported . LVEF, QRS duration, QT interval, QT dispersion, signal-averaged electrocardiogram, as well as, short-term HRV were concluded as with limited predictive value. Long-term ECG recordings provide information about both, ventricular arrhythmias and long-term HRV. Ventricular arrhythmias were found as an abundant and recommended SCD predictor in post-MI patients. Depressed HRV was confirmed to be a clear predictor of overall mortality in several studies such as also reported in Table 2. In particular in non-arrhythmic mortality long-term HRV may be a better predictor. Finally, HRT, which assesses the autonomic response modulating the heart rate, turned out as promising candidate for SCD prediction. The latter two results, namely long-term HRV and HRT, confirm the relevance of considering the autonomic control and dysfunction, respectively, as independent predictors of SCD.
In conclusion, the key results of the studies summarized in Table 2 clearly show the independent prognostic value of the autonomic nervous system related indices in extension to clinical and functional variables. Nevertheless, other influences such as age, gender, history, NYHA class, treatment (thrombolysis, beta-blocker, ACE inhibitors), major depression, etc. were partly considered as confounders of HRV indices.
The consistent prognostic value of long term variability and correlations (SDNN, SDNNindex, VLF, beta-slope) may reflect the importance of the complex adjustment of several underlying mechanisms such as vagal and sympathetic rhythms (Table 2). Those results confirm the importance of complex autonomic adjustments and thus strongly motivate applying AIF analysis in the risk assessment in CHF and MI patients . Furthermore, these results support the importance of a complex systemic analysis, reflecting both the particular cardiac disorder and compensation/adaptation by autonomic and other systemic mechanisms. In conclusion, the established linear HRV analysis provides sensitive markers of the cardiovascular-systemic response to cardiac dysfunctions, which can significantly be extended by novel nonlinear methods.
In stroke as the third most important cause of death in industrial countries a variety of ANS disturbances can occur. They are attributed to the damage of the different compartments of the central autonomic network, i.e. frontoparietal cortical areas and brain stem (see above, and Fig. 5) or may result from disrupted autonomic pathways (hypothalamus-mesencephalon-pons-medulla oblongata-spinal cord). Among the ANS dysfunctions in stroke victims cardiovascular sympathetic/ parasympathetic dysbalance is most evident in the acute phase of stroke. This cardiovascular as well as infection related morbidity within the first days after stroke is very relevant for survival and neurological outcome after stroke .
Stroke related autonomic disturbances and their clinical relevance.
Since the early 1990th studies in ischemic brain infarction revealed frequent ANS dysfunctions in association with pathological activation of the sympathetic system. Such a sympathetic hyperactivation and autonomic imbalance in acute ischemic stroke patients resulted in increased risk and incidence of severe cardiac arrhythmias and sudden cardiac death [53-56].
Based on NN interval analysis Colivicchi et al. showed in 208 stroke patients the best predictive value for 1 year survival using SDNN, however the occurrence of ventricular tachycardia had an even higher predictive value . Interestingly, stroke localization seems to play an important role in this context: HRV dysregulation (i.e. decreased LF, HF and SDNN) was more frequent in association with right hemispheric and insular stroke independent of the use of either time frame or spectral analysis. Right hemispheric brain infarction with insular involvement is associated with the highest degree of ANS dysfunction (measured as decreased HRV and increased heart rate) . This ANS dysfunction was associated with higher incidence of sudden death .
The prognostic importance of a sympathetic activation post stroke was demonstrated in more than 100 patients studied by . An elevated norepinephrine concentration was associated with insular involvement, nocturnal elevated blood pressure and prolonged QTc interval in the ECG. One year post stroke these patients showed a higher incidence of cardio- and cerebrovascular events, a lower degree of neurological restitution of stroke symptoms in conjunction with lower rate of independence .
Relevant studies analyzing stroke associated changes in ANS activity and of those assessing the prognostic value of HRV indices after stroke based on 24 hours ECG recordings are reviewed in Table 3.
Clinical Studies of Heart Rate Related ANS Dysfunction in Stroke Patients
|Author||Data Sample Set||Results|
|Barron et al., 1994 ||40 stroke patients 4-11 days after event vs. 40 matched controls, power spectrum analysis of HRV, all pts. in sinus rhythm.||Cardiac parasympathetic innervations reduced after stroke; more pronounced after right sided infarct.|
|Naver et al. 1996 ||23 stroke / 11 TIA patients vs. 21 matched controls, no medication except aspirin, follow up 8-48 days (stroke) up to 2 years (TIA), respiratory HRV, orthostatic test, isometric handgrip.||Right sided stroke associated with reduced respiratory HRV, peripheral reflexes equal between right an left lesions.|
|Korpelainen et al. 1997 ||32 stroke patients vs. 32 matched controls; 24h Holter, time and frequency domain acute and after 6 months.||LF/HF ratio unchanged during night-time in stroke patients reversible abolition of circadian HRV rhythm, loss of relative vagal nocturnal dominance.|
|Giubilei et al., 1998 ||10 patients with MCA infarct vs. controls, 1 day after stroke, power spectrum analysis of HRV during sleep.||Sympathetic/parasympathetic balance (VLF+HF/HF), higher in patients. Sympathetic predominance during sleep.|
|Korpelainen et al., 1999 ||46 stroke patients (early: 1-7days poststroke; follow up at 6 month) vs. 30 controls. 24h Holter.||Abnormal HRV measures involved in cardiac complications. SDNN, VLF, LF impaired in acute and late state post stroke.|
|Tokgözoglu et al., 1999 ||62 stroke patients vs. 62 matched controls. 72h ECG recording post stroke. Power spectrum analysis.||Decreased LF, HF, SDNN. Decreased sympathetic and parasympathetic parameters, most pronounced in right insular cortex infarcts.|
|Colivicchi et al., 2004 ||103 stroke patients, 24h Holter, time and frequency domain HRV.||Lower SDNN, rMSSD; higher LF/HF ratio in right insular infarct patients in association with more frequent and complex arrhythmias.|
|Mäkikallio et al., 2004 ||84 stroke patients, 24h Holter, follow up 7 years, power-law slope beta (beta<-1.5).||Conventional HRV measures without prognostic power, but abnormal long-term HR dynamics predict poststroke mortality.|
|Meyer et al., 2004 ||29 stroke patients within 24h after event, HR and BP.||Pathological sympathetic activation more pronounced after right sided stroke.|
|Colivicchi et al., 2005 ||208 stroke patients, 24h Holter, 12 month follow up.||Decreased SDNN, right insular infarct and non sustained ventricular tachycardia associated with higher risk of early mortality.|
|McLaren et al., 2005 ||Cross sectional case-control study in 76 stroke patients (9 months poststroke) vs. 70 matched controls; long term effects; power spectrum analysis, baro reflex sensitivity, cardiovascular reflex autonomic test.||Impaired autonomic function may increase mortality in older stroke survivors.|
|Bassi et al., 2007 ||85 stroke patients, 24 h Holter, time domain measure.||Lower SDNN associated with unfavourable outcome after rehabilitation.|
|Dütsch et al., 2007 ||28 stroke patients vs. 21 matched controls, RR intervals, BP means oscillations.||Post stroke parasympathetic cardiac deficit and increased sympathetic cardiovascular modulation. LF/HF ratio of RR intervals elevated in right hemispheric vs. left sided stroke.|
|Rufa et al., 2007 ||23 CADASIL patients vs. 22 controls, frequency domain measures of HRV, QTc Interval.||LF/HF ratio higher in CADASIL patients. No correlation with MRI lesion volume.|
Predictors for a lethal outcome are diminished heart rate variability (HRV), insular cortex involvement, and the occurrence of complex arrhythmias within 72 h after admission [55,57]. A decreased nocturnal vagal dominance leading to hypertension and decreased HRV during night-time resulted in a long lasting autonomic dysfunction and thus chronic risk . The standard deviation of all normal to normal R-R (NN) intervals (SDNN) is a very strong prognostic predictor, together with the cardiac arrhythmia rate [55,57,60,61].
In summary, a number of experimental and clinical studies have revealed that ANS dysfunction represents a relevant pathogenetic factor and indicator for post stroke morbidity and mortality. Thus, using HRV related indices in acute and postacute stroke victims may improve prognostic evaluation and also counseling of patients, and may thus eventually lead to new preventive strategies recognizing the important relationship of ANS dysfunction and central nervous system injury.
Heart rate patterns which can easily be obtained from Holter ECG recordings reflect important characteristics of the complex autonomic control. They are not only associated with cardiac arrhythmias, but implicate sympathetic and parasympathetic modulations.
In stroke and cardiovascular disease as well as in other very relevant clinical constellations, such as in sepsis and MODS (multiple organ dysfunction syndrome; e.g. [14,25] a selective pathological disturbance leads to a complex systemic compensation and adaptation process. Consequently, final decompensation and respective risk assessment might be related rather to an overall systemic response than to a particular mechanism, unless this is fatally altered itself. Therefore, it is not surprising, that for heart failure / MI patients the comprehensive HRV measures like SDNN, VLF, beta slope, complexity, and AIF measures were found to have the best predictive value (as reviewed above). On the other hand, specific indices, such as the shifted adjustment between sympathetic and vagal activity (sympatho-vagal (im)-balance) due to stroke and cardiovascular diseases provided predictive value with more specific physiological implications.
From the clinical point of view the main question is, whether HRV indices provide additional pathophysiological and predictive value going beyond the conventional risk factors and clinical parameters which are so far not completely satisfying in particular in severely and acutely ill patients. This objective becomes even more challenging since multiple “cardiovascular” confounders, such as diabetes, age, daytime, as well as pharmacological interventions like beta-blockers, angiotensin-converting enzyme (ACE) inhibitors, and statins which modulate autonomic tone and systemic proinflammatory responses  have to be considered.
Besides the autonomic modulations, which are assessed following the established HRV guidelines , cardiac arrhythmias such as atrial fibrillation, sinoatrial dysfunction and ectopic complexes are of essential prognostic value and their patterns are subject of current research. Especially after stroke the HRV-index SDNN and the presence of non-sustained ventricular tachycardia were independent predictors of the 1-year mortality . Cardiogenic arrhythmias and autonomic modulations cannot be completely separated, because of the complex autonomic-cardiovascular system, but both aspects have to be considered in a way appropriate for the respective situation, such as sinus rhythm or cardiac arrhythmia. Furthermore, stroke related immunodepression  has to be considered as another important player in this relationship between ANS, CNS and cardiovascular system.
Due to the multiple influences on the autonomic function and the rapid development of confounding or targeted therapies large scale clinical studies are immediately necessary to further evaluate the clinical value of the above mentioned HRV indices in cardio- and cerebrovascular risk assessment und management.
We thank Dr. Hischam Taha, Frenchay Hospital, Bristol, for helpful comments.