Next Article in Journal
A Double-Blind, Randomized, Placebo-Controlled Trial to Evaluate the Efficacy of a Hydrolyzed Chicken Collagen Type II Supplement in Alleviating Joint Discomfort
Next Article in Special Issue
Fatty Acids Quality in Middle Eastern Traditional Dishes, Arabic Sweets and Market Foods Frequently Consumed in Lebanon
Previous Article in Journal
Nutrition Strategy and Life Style in Polycystic Ovary Syndrome—Narrative Review
Previous Article in Special Issue
Beneficial Outcomes of Omega-6 and Omega-3 Polyunsaturated Fatty Acids on Human Health: An Update for 2021
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Heart Rate Variability and Long Chain n-3 Polyunsaturated Fatty Acids in Chronic Kidney Disease Patients on Haemodialysis: A Cross-Sectional Pilot Study

1
Department of Nutritional Sciences, School of Life Course Sciences, King’s College London, 150 Stamford Street, London SE1 9NH, UK
2
School of Exercise & Nutrition Sciences, Queensland University of Technology, Brisbane, QLD 4059, Australia
*
Author to whom correspondence should be addressed.
Submission received: 25 June 2021 / Revised: 12 July 2021 / Accepted: 13 July 2021 / Published: 18 July 2021
(This article belongs to the Special Issue Dietary Fat and Human Health Outcomes)

Abstract

:
Low heart rate variability (HRV) is independently associated with increased risk of sudden cardiac death (SCD) and all cardiac death in haemodialysis patients. Long chain n-3 polyunsaturated fatty acids (LC n-3 PUFA) may exert anti-arrhythmic effects. This study aimed to investigate relationships between dialysis, sleep and 24 h HRV and LC n-3 PUFA status in patients who have recently commenced haemodialysis. A cross-sectional study was conducted in adults aged 40–80 with chronic kidney disease (CKD) stage 5 (n = 45, mean age 58, SD 9, 20 females and 25 males, 39% with type 2 diabetes). Pre-dialysis blood samples were taken to measure erythrocyte and plasma fatty acid composition (wt % fatty acids). Mean erythrocyte omega-3 index was not associated with HRV following adjustment for age, BMI and use of β-blocker medication. Higher ratios of erythrocyte eicosapentaenoic acid (EPA) to docosahexaenoic acid (DHA) were associated with lower 24 h vagally-mediated beat-to-beat HRV parameters. Higher plasma EPA and docosapentaenoic acid (DPAn-3) were also associated with lower sleep-time and 24 h beat-to-beat variability. In contrast, higher plasma EPA was significantly related to higher overall and longer phase components of 24 h HRV. Further investigation is required to investigate whether patients commencing haemodialysis may have compromised conversion of EPA to DHA, which may impair vagally-mediated regulation of cardiac autonomic function, increasing risk of SCD.

1. Introduction

Cardiovascular disease is particularly prevalent in chronic kidney disease (CKD) [1], and accounts for 43% of all-cause mortality among dialysis patients [2], compared to an estimated 11% of deaths in the general population [3]. Risk of sudden cardiac death (SCD) is doubled when a patient with CKD stage 5 (kidney failure) starts dialysis [4] and in patients on haemodialysis, SCD accounts for two thirds of all cardiac deaths and one fourth of all-cause mortality [2]. Patients with kidney failure are a population with low heart rate variability (HRV) indicating autonomic dysfunction [5]. Low HRV has been independently associated with increased risk of all-cause mortality, SCD and all cardiac death in haemodialysis patients [6,7,8]. Patients on haemodialysis have lower serum long chain n-3 polyunsaturated fatty acids (LC n-3 PUFA) levels compared to control populations without CKD [9,10], and lower serum LC n-3 PUFA: arachidonic acid (AA) ratios were predictive of CVD events [10]. Fish consumption has been associated with reduced mortality in an observational study with incident dialysis patients [11], and Friedman et al. reported an inverse association between erythrocyte LC n-3 PUFA levels and mortality in haemodialysis patients [12]. Supplementation with LC n-3 PUFA is a promising candidate for SCD prevention in this population due to their potential anti-inflammatory, anti-oxidative and anti-arrhythmic effects on the cardiac myocytes [13,14,15]. Mixed results have been obtained in previous investigations into the effect of LC n-3 PUFA on HRV in healthy subjects but there is stronger evidence in myocardial infarction survivors of benefits from fish oil supplementation [16]. Few studies have investigated the relationship between LC n-3 PUFA and HRV in CKD patients [17,18]. Identifying positive associations between tissue LC n-3 PUFA status and HRV in patients with kidney failure would provide important evidence to support the potential application of LC n-3 PUFA supplementation to improve cardiac autonomic function and help prevent SCD.
It is hypothesised that the proportion of eicosapaentanoic acid + docosahexaenoic acid (EPA + DHA) in the erythrocyte will be independently positively associated with 24 h HRV in patients with CKD stage 5 who have recently commenced haemodialysis. This study aimed to describe LC n-3 PUFA status in patients who recently commenced haemodialysis, and to investigate relationships between erythrocyte membrane LC n-3 PUFA composition and HRV in these patients, including comparison to a reference healthy population. Associations between plasma LC n-3 PUFA composition and HRV will also be explored.

2. Materials and Methods

This was a cross-sectional pilot study with ethical approval obtained from NRES Committee London–Camberwell St Giles (REC ref: 14/LO/0186). The study was registered on clinicaltrials.gov (NCT02014792) and conducted at King’s College Hospital (KCH) and Guy’s and St Thomas’ Hospital (GSTT) and their 10 respective satellite dialysis units in the London area. All subjects gave their informed consent for inclusion before they participated in the study. The study was conducted between July 2014 and January 2016, in accordance with the rules of the Declaration of Helsinki of 1975, revised in 2013.
The primary dependent variables were measures of 24 h overall HRV during haemodialysis: SDNN and triangular index (defined below). The primary independent variable was the proportion of EPA + DHA (Omega-3 Index, O3I; %) in the erythrocyte membrane. Secondary dependent outcome variables included other time- and frequency-domain, and non-linear HRV parameters, during 2 h dialysis, 24 h and sleep-time. Secondary independent outcome variables included erythrocyte and plasma proportions (%) of EPA, docosapentaenoic acid (DPAn-3) and DHA. Other descriptive variables included background dietary intakes (using the EPIC-Norfolk food frequency questionnaire; data was analysed using an in-house customized Microsoft Access database and exported to Microsoft Excel for further processing. The individual LC n-3 PUFA intake calculations were determined using the nutrient composition values from the McCance and Widdowson’s ‘The Composition of Foods Integrated Dataset’ with additional data obtained from manufacturers and USDA values when required.) and clinical risk factors, including history of sleep apnoea or suspected sleep apnoea (using Epworth Sleepiness Scale and Berlin Questionnaire), waist circumference, 12 month average SBP and DBP and BP variability prior to starting haemodialysis, hs-C reactive protein (CRP) at start of dialysis session, serum creatinine, albumin, sodium, potassium, phosphate and calcium before and after the dialysis session.
The inclusion criteria were: women and men with CKD stage 5 commencing haemodialysis (6 to 10 weeks after the first session), aged 40–80 years and able to provide written informed consent. The exclusion criteria included: unwilling to participate, history of chronic liver disease or neuropathy, infection or antibiotics within the last month and history of prior renal replacement therapy including transplantation.
Participants were patients under the care of a nephrologist. They were identified through the course of usual clinical care, or via electronic database search conducted by the project dietitian and nephrologist at King’s College Hospital NHS Foundation Trust and Guy’s and St Thomas’ NHS Foundation Trust, respectively. Patients had haemodialysis treatment on three non-consecutive days a week. Eligible patients enrolled in the study 6 to 10 weeks after commencing haemodialysis treatment and two study days, at least 1 week apart, were scheduled to yield duplicate data. Study days were not scheduled right after the two consecutive days of intradialytic period, also known as the “long gap”, to avoid even greater elevations in uraemia and to maintain consistency across duplicate study measurements.
In both study visits, pre- and post-dialysis blood sample were taken by a nurse (Appendix A). Waist circumference was taken on patient’s arrival on the first study visit before going on the dialysis machine. Immediately after the nurse started the patient on haemodialysis treatment, the researcher fitted the monitor to record HR and HRV for a 24 h period. Additional information was obtained from the medical records. All patients continued with the best medical care, as appropriate, at the nephrologist’s discretion.
Heart rate variability was measured using the Actiheart monitoring equipment (CamNtech Ltd., Cambridge, UK). Data processing of the HRV recordings was carried out using the Actiheart software (version 4.0.91, CamNtech Ltd., Cambridge, UK) and Kubios HRV analysis software (Biosignal Analysis and Medical Imaging Group, Department of Physics, University of Kuopio, Finland) where all artefacts were excluded from the analysis. Patients were given an activity and event diary to complete with information about the duration of dialysis, type of activities done, and nap and sleep times at night. HRV and HR/interbeat interval (IBI) data were analysed for a standardised 2 h of haemodialysis (from the moment the monitors were fitted), for 24 h (in order to be considered valid and used for analysis a minimum of 18 h of the recording had to be analysable) and sleep time. HRV outcomes included time and frequency-domain parameters; time-domain parameters are based on the time intervals between adjacent QRS (Q, R and S being points on the R wave seen on an ECG during ventricular depolarisation, and R being the peak upward deflection) complexes (normal-to-normal (NN) intervals) whereas frequency-domain parameters employ power spectral analysis of NN intervals to determine the power (variance) within frequency bands [19]. Time-domain parameters included SDNN, standard deviation of the average normal-to-normal intervals in 5 min segments of the whole recording (SDANN), square root of the mean of the sum of the squares of differences between adjacent NN intervals (RMSSD), the percentage of adjacent normal-to-normal intervals that differed by >50 ms (pNN50) and triangular index (TI), the integral of the density distribution (the number of all NN intervals) divided by the maximum of the density distribution. Frequency-domain parameters included high-frequency (HF), low-frequency (LF) and very-low-frequency (VLF) power, and the ratio of the LF and HF band powers (LF:HF). A non-linear parameter using Poincaré plots of short-term variability (SD1) against long-term variability (SD2) was also calculated as a measure of complexity of HRV distribution over the duration of the recording. SDNN, LF and TI represent overall variability. Short-term (beat-to-beat) components of HRV include RMSSD, pNN50 and HF. SDANN and VLF reflect longer-phase components of variability.
Pre- and post-dialysis blood samples collected in SSTTM tubes (BD VacutainerTM, Franklin Lakes, NJ, USA) were used to determine serum sodium, potassium, phosphate and calcium as well as albumin and creatinine before and after dialysis and hs-CRP before dialysis only. These analyses were performed by a clinical pathology accredited biochemistry laboratory (ViaPath, Kings College Hospital). Pre-dialysis blood samples collected in EDTA tubes were used to determine erythrocyte membrane fatty acid composition (% weight) and plasma fatty acid composition (% weight) at King’s College London. Proportions of fatty acids in whole plasma and erythrocyte membranes were analysed by GC (Agilent 7890A GC; Agilent Technologies, Santa Clara, CA, USA) with a BPX70 GC column (length 25 m, internal diameter 0·32 mm, film thickness 0·25 μm) custom designed for separation of fatty acid methyl esters (SGE Analytic Science, Ringwood Victoria, Australia) following transesterification, as previously described [20], but substituting toluene for benzene and using pentadecanoic acid as an internal standard. The O3I was defined as the sum of % weight EPA + DHA in erythrocytes.
Statistical analyses were carried out using IBM SPSS Statistics 21.0 (Statistical Product and Service Solutions; IBM Corp., Armonk, NY, USA) and statistically significance was considered at p-value < 0.05. Values from first and second visit were averaged in all subjects. Simple linear regression was used to analyse relationships between LC n-3 PUFA status (O3I, erythrocyte and plasma EPA, DHA and DPA, and erythrocyte EPA:DHA ratio) and HRV variables. Multiple linear regression was used to control for confounding factors, age, BMI and use of β-blockers, which are all related to HRV. Non-normally distributed data were normalised by natural logarithm (LN) before performing the analysis. Standardised residuals were checked for normality. Chi-square (χ 2) test and independent t-tests were used to compare categorical and continuous variables, respectively, between haemodialysis patients and a healthy cohort. Non-normally distributed data were normalised by LN transformation (results shown as geometric means and 95% CI) before analysis by independent t-test. If LN transformation failed to yield a normal distribution, a Mann–Whitney U test was applied to compare groups (results shown as medians with lower and upper quartiles).

3. Results

3.1. Patient Characteristics

A total of 49 patients were enrolled, 45 patients completed the first study visit and 31 completed both study visits (Figure 1). Of the 45 patients, 39 patients provided HRV recordings of sufficient quality to be included in the 2 h dialysis period, 40 provided sleep-time HRV recordings of sufficient quality (minimum 3 h period), and 35 provided 24 h recordings of sufficient quality for analysis (minimum 18 h period). There was only one patient where it was not possible to obtain erythrocyte membrane fatty acid composition due to technical issues.
Demographic and clinical characteristics for the total population enrolled in the study is presented in Table 1, including 2 h HRV parameters during dialysis. Most patients presented hypertension (90.8%) and the second most common co-morbidity was type 2 diabetes (38.6%). Nearly half the patients were taking statins, β-blockers, erythropoietin therapy and iron injections, and more than half were taking vitamin D supplements. Protein, fat and carbohydrate made up 18 (5), 36 (8) and 45 (7)% (medians with IQR) of energy intakes respectively. Median total sugar intake was 75 g/d (18% (7) energy intake) and mean saturated fatty acids (SFA) intakes were 13% (SD 3).

3.2. LC n-3 PUFA and 2 h Dialysis, Sleep-Time and 24 h HRV

For participants with complete duplicate visits there were significant moderate to strong correlations between visits 1 and 2 for all HRV parameters: dialysis, n = 23, 0.578 < r < 0.880, p < 0.01; sleep time, n = 24, 0.468 < r < 0.913, p < 0.05; and 24h, n = 19, 0.630 < r < 0.812, p < 0.01. Correlations between study visits 1 and 2 for LC n-3 PUFA were also moderate to strong: erythrocyte EPA, r = 0.687, p < 0.001; erythrocyte DHA, r = 0.497, p = 0.010; and erythrocyte O3I, r = 0.538, p = 0.005 (all n = 26); plasma EPA, r = 0.729, p < 0.001; plasma DPA(n-3), r = 0.907, p < 0.001; and DHA, r = 0.885, p < 0.001 (all n = 31).
The median [lower and upper quartile] erythrocyte O3I of this population was very low (3.3%; 2.5, 3.8); only 2 individuals had O3I > 5%. Median erythrocyte EPA [lower and upper quartile] was 0.7% [0.5, 0.9] and DHA was 2.5% [2.0, 3.2]. Median plasma EPA was 0.8% [0.5, 1.2], DPA(n-3) was 0.4% [0.4, 0.5], and DHA was 1.7% [1.3, 2.2]. There was a strong significant correlation between plasma and erythrocyte EPA (r = 0.704; p < 0.001), and between plasma and erythrocyte DHA (r = 0.608; p < 0.001). No correlations were found between plasma or erythrocyte LC n-3 PUFA proportions and estimated PUFA, EPA, DHA or oily fish intake assessed by FFQ.
Mean 24 h HR was 83 bpm (SD 10.9) and IBI was 747 ms (SD 131). Twenty-four hour overall HRV was low: 84 ms (SD 26) for SDNN and median Ti was 17 [lower and upper quartiles 15 and 27]. Mean and median parameters of long-phase and short-phase variability are presented in Appendix B.

3.3. LC n-3 PUFA and HRV in Haemodialysis vs. Reference Healthy Population

A paired sub-group of a healthy cohort from a previously reported randomized controlled trial, the MARINA study [21], were used as a reference population. Individuals from the current haemodialysis cohort were pair-matched by age, sex and BMI with individuals from the healthy reference population. Table 2 displays the results comparing LC n-3 PUFA intake and blood composition, hs-CRP, and HRV parameters during sleep time. Compared to the healthy cohort, patients in the haemodialysis cohort had significantly higher sleep-time HR and lower HRV, higher hsCRP, as well as significantly lower proportions of erythrocyte membrane and plasma LC n-3 PUFA, and higher erythrocyte EPA:DHA ratio, whilst EPA and DHA dietary intakes were not significantly different. There were other notable differences in erythrocyte fatty acid composition, with significantly higher proportions of MUFAs, cis-vaccenic acid (18:1n-7), oleic acid (18:1n-9), and palmitoleic acid (16:1n-7); SFAs, stearic acid (18:0) and palmitic acid (16:0), and the erythrocyte palmitic acid (PA):linoleic acid (LA, 18:2n-6) ratio in patients in the haemodialysis cohort compared to the healthy reference group. Furthermore, there were significantly lower proportions of erythrocyte n-6 PUFAs, docosapentaenoic acid (DPA; 22:5n-6), arachidonic acid (20:4n-6), dihomo-γ-linolenic acid (DGLA; 20:3n-6), and linoleic acid (18:2n-6) in patients in the haemodialysis cohort compared to the healthy reference group.

3.4. Associations between LC n-3 PUFA and HRV in Haemodialysis Patients

3.4.1. Erythrocyte Membrane LC n-3 PUFA

No significant associations were observed between erythrocyte O3I, EPA or DHA (%) and HRV parameters during dialysis, sleep nor 24 h following adjustment (Table 3). EPA:DHA ratios were inversely associated with RMSSD (adjusted standardised β = −0.377, p = 0.044), LF power (adjusted standardised β = −0.289, p = 0.043), HF power (adjusted standardised β = −0.436, p = 0.008), and SD1:SD2 (adjusted standardised β = −0.364, p = 0.047).

3.4.2. Plasma LC n-3 PUFA

Plasma DHA (%) was positively associated with Ti during dialysis (adjusted standardised β = 0.331, p = 0.022) (Table 4) but no significant associations were observed between plasma EPA or DPA and any of the HRV parameters during dialysis.
During sleep, DPA(n-3) was negatively associated with pNN50 (adjusted standardised β = −0.321, p = 0.041) and SD1:SD2 (adjusted standardised β = −0.330, p = 0.045), both indicators of short-term variability, and EPA was also negatively associated with SD1:SD2 (adjusted standardised β = −0.354, p = 0.036). DHA was not related to sleep-time HRV, nor 24 h HRV.
Significant moderate positive associations were observed between plasma EPA and longer phase components of 24 h HRV (SDANN adjusted standardised β = 0.405, p = 0.016; VLF power adjusted standardised β = 0.481, p = 0.004) as well as overall 24 h HRV (SDNN adjusted standardised β = 0.441, p = 0.009). A moderate significant negative association between EPA and 24 h SD1:SD2 was observed (adjusted standardised β = −0.482, p = 0.004), and between plasma DPA(n-3) and 24 h short-term variability indicators RMSSD (adjusted standardised β = −0.387, p = 0.024), and pNN50 (adjusted standardised β = −0.394, p = 0.021).

4. Discussion

Arrhythmias and cardiac arrests account for over a quarter of deaths among haemodialysis patients [22] and the first three months following commencement of haemodialysis treatment are particularly high-risk for SCD due to factors relating to the parameters of the dialysis itself that may trigger arrhythmia, but also the high prevalence of pre-existing cardiovascular disease and diabetes [23]. This study set out to characterise erythrocyte membrane and plasma LC n-3 PUFA in a UK patient population who have recently commenced haemodialysis, as well as exploring relationships with HRV parameters during this critical period in their treatment. Previous reports on haemodialysis patients in other countries (Italy, Denmark, USA, Germany and Serbia) have shown lower EPA, DHA, or both EPA + DHA, in erythrocyte membranes and/or plasma compared to healthy controls or healthy populations [9,24,25,26,27,28]. Shoji et al. showed greater levels of EPA and DHA in plasma in haemodialysis patients compared to healthy controls in Japan [10] and Friedman et al. reported a greater erythrocyte DHA and O3I, but lower plasma DHA, in haemodialysis patients (n = 75) compared to a healthy control group (n = 25) in the USA [29]. Compared to other haemodialysis cohorts, the mean levels of DHA and O3I in the present study population were amongst the lowest, being 61% and 57% lower, respectively, than a matched healthy cohort.
Erythrocyte membrane fatty acid composition is considered to be a better marker of cardiovascular disease risk than plasma, as erythrocyte membrane phospholipids are more reflective of the cardiomyocytes fatty acid profiles [30]. None of the patients in this study presented an O3I above 8% which has been associated with the greatest cardio-protection; mean O3I was below 4%, which has been associated with the least cardio-protection [31]. Background diet may play a role, but this was not supported by food frequency questionnaire data on self-reported fish intake. It is most likely that the major contributing factor to low erythrocyte LC n-3 PUFA content was related to the pathology of kidney disease. Anaemia is prevalent in stage 5 CKD due to reduced production of erythropoietin, a hormone produced by the kidney that regulates erythrocyte production from bone marrow (necessitating erythropoietin therapy in most haemodialysis patients), and also as a result of uraemic inhibitors of erythrocyte production [32]. Other contributors to anaemia include reduced lifespan of erythrocytes due to oxidative stress, inflammation and haemolysis, which may be a limitation in using erythrocyte membrane fatty acid composition as a proxy marker of tissues of interest such as cardiac myocytes or neurons. Thus, although the utility of erythrocyte fatty acid composition as a biomarker of moderate-long term dietary intake has been established in other populations [33], the same relationships with dietary intake and tissue status cannot be assumed in populations in advanced stages of CKD. These relationships are likely to be further disrupted in patients with kidney failure due to haemodialysis-related factors such as chronic inflammation, oxidative stress and altered lipid metabolism/uptake into membranes. Our data show that variability in erythrocyte LC n-3 PUFA contents is indeed limited to a small range, and it is yet to be established whether this is an accurate representation of general tissue status in haemodialysis patients compared to other sample types such as leucocytes and adipose tissue.
The primary HRV dependent outcome variables, 24 h SDNN and Ti, did not correlate with O3I neither did any of the other HRV parameters during dialysis, sleep-time or 24 h, following adjustment for age, BMI and use of β-blockers. Few published studies to date have assessed the relationship between LC n-3 PUFA status and HRV in stage 5 CKD patients on dialysis. Christensen et al. found positive correlations between granulocyte EPA, DHA and total LC n-3 PUFA content and 24 h SDNN in 17 patients on haemodialysis or peritoneal dialysis after a fish-oil vs. olive oil dietary intervention [17]. Svensson et al. measured LC n-3 PUFA in serum phospholipids in 30 patients and found no correlations with 24 h HRV measurements at baseline [18]. Rantanen et al. measured LC n-3 PUFA in plasma phospholipids in 135 patients and found no significant associations with 24 h time-domain or frequency-domain HRV [34].
We report that plasma EPA was significantly positively correlated with overall and longer phase components of HRV and negatively correlated with vagally-modulated beat-to-beat variability over 24 h, in contrast to the complete lack of associations between 24 h HRV and plasma DHA. The fact that a greater proportions of plasma EPA (and DPAn-3, the next fatty acid in the interconversion cascade from EPA to DHA) were associated with reduced beat-to-beat variability might suggest that cell membrane uptake and/or conversion of EPA to DHA is inhibited in patients on haemodialysis with lower parasympathetic tone. In fact, erythrocyte EPA to DHA ratios were significantly higher in this haemodialysis cohort compared to an age-, sex- and BMI-matched cohort of healthy participants and within the haemodialysis cohort, erythrocyte EPA:DHA ratio was inversely associated (fully adjusted) with 24 h RMSSD, HF power, LF power and SD1:SD2 ratio, all markers of vagal regulation of heartbeat [19]. Reduced availability of DHA is likely to be directly involved in impaired autonomic function by a number of mechanisms, including adverse effects on neuronal membrane properties, increasing neuroinflammation, reducing neural regeneration, and failure to inhibit oxidative stress, as well as a likely direct adverse effect on pacemaker activity [35]. These results indicate that a greater EPA proportion in plasma, is associated with better overall and longer-phase autonomic regulation but not short-term vagally-driven variability (alongside low DHA). It might be speculated that greater availability of EPA-derived lipid mediators may be involved in increased responsivity of HR to factors such as circadian fluctuations in hormone secretion and baroreflex (reflected by long-phase HRV). Reduced parasympathetic tone could be due to reduced availability of DHA in the membranes of cardiomyocytes and/or a reduction in DHA-derived lipid mediators involved in neuroprotective processes in this patient population.
This study confirms previous reports that HRV is substantially reduced in haemodialysis patients compared to healthy controls [36,37] and extends these findings by demonstrating consistently that time-domain and frequency-domain HRV parameters were significantly lower in haemodialysis patients compared to age-, sex- and BMI-matched healthy controls. Patients in this study also had SDNN measurements amongst the lowest reported in comparison to other haemodialysis cohorts [6,8,18,38]. This might be related to the fact that these were patients that recently started haemodialysis (measurements taken between 6 and 10 weeks of commencing treatment). Commencing haemodialysis is associated with autonomic dysfunction; lower SDNN values were reported in patients on haemodialysis for less than 30 months compared to those on haemodialysis for more than 30 months [37]. Moreover, after the start of haemodialysis treatments autonomic dysfunction of patients has been reported to improve or remain unaltered compared to patients with chronic renal failure before starting haemodialysis [39]. Decreases in Ti independently predict cardiac death [6], and Hayano et al. reported an independent prognostic value of Ti (<22) to predict both all-cause and sudden death (hazard ratio (95% CI); 8.1 (1.3, 48.6) and 12.6 (1.3, 126.4), respectively) [7]. In the present study 22 patients presented Ti < 22 ms which corresponds to 65% of those with valid 24 h HRV readings. Hayano et al. also reported a 24 h SDNN below 50 ms to predict a greater risk of sudden death in 31 chronic haemodialysis patients with stage 5 CKD after 5 years of follow-up [7]. Although only 2 patients presented SDNN under 50 ms in the present study, the mean 24 h SDNN of this study population was 84 ms; considered to be moderately depressed. This finding is consistent with other published studies reporting depressed HRV in haemodialysis patients [6,8,18,38,40,41,42,43]. Effects of medication on HRV should be considered in this poly-medicated patient population. For example, those taking β-blockers (n = 16) had a significantly lower HR than non-users (78 vs. 87 respectively) and higher HF (161 vs. 81), as might be expected as evidence shows that treatment with β-blockers increases HRV [44,45,46].
The inflammatory burden in CKD is amplified upon commencing haemodialysis treatment [47]. CRP has been previously positively associated with mortality in haemodialysis patients [48]. The haemodialysis patients in this study presented significantly higher hs-CRP values compared to the healthy cohort, raising the question of whether very high chronic inflammation is a possible underlying cause of low LC n-3 PUFA and impaired HRV. One study reported an inverse relationship between IL-6, and SDNN, SDANN and VLF in CKD patients, but in dialysis patients, no associations between HRV and inflammatory markers were found [49]. Studies with a larger sample size showed independent associations between HRV and CRP in healthy participants and in patients with CHD [50] and a meta-analysis showed that CRP predicted CHD independently of traditional cardiovascular risk factors [51].
Where erythrocyte LC n-3 PUFA and n-6 PUFA contents as a proportion of total fatty acids were lower in haemodialysis patients compared to the reference group, this was met by a reciprocal increase in the proportions of MUFA and SFA in red blood cell membranes. Membrane SFA and MUFA are poor indicators of dietary fatty acid composition [52]. The fatty acids that were higher in erythrocyte membranes of the haemodialysis group–16:0, 16:1n-7, 18:0, 18:1n-7, and 18:1n-9–are all endogenously synthesised and some may be associated with metabolic dysregulation [53]. Higher erythrocyte membrane levels of 16:0, 16:1n-7, and 18:1n-7 were associated with greater risk of SCD in a US population-based case-control study [54]. By altering membrane fluidity and impairing potassium ion channel function, it is possible that palmitic, palmitoleic, and vaccenic acid-enrichment of cardiomyocyte membranes might also increase risk of SCD in haemodialysis patients independently of depletion in tissue LC n-3 PUFA and n-6 PUFA status.
This study presents several limitations, including multiple testing. However, the relationships observed were consistent across HRV parameters measuring either short-term, long-term or overall HRV making type 1 errors less likely to be a concern. There were challenges in the participant recruitment, compliance and retention of patients for the second visit. In addition, the patients enrolled in this study had different treatment schedules with patients dialysing at different times of the day (morning, afternoon or evening) which may have induced circadian differences in the readings. Another limitation of this study is the number of comorbidities presented by most patients which could have affected independently the HRV measurements not allowing the detection of potential associations with LC n-3 PUFA status.
Patients with stage 5 CKD are at high risk of SCD and there are numerous mechanisms that could be involved in the electrophysiological instability that predispose these patients to malignant cardiac arrhythmias. These patients have a moderate to severe cardiac autonomic dysfunction and low tissue LC n-3 PUFA status may contribute to increased risk of arrhythmic cardiac events and death. Research is needed to investigate the potential efficacy of EPA+DHA or DHA supplementation in this patient population for prevention of sudden cardiac death. The main finding from this feasibility study is that erythrocyte LC n-3 PUFA are uniformly very low in patients who have recently commenced haemodialysis, with %DHA being disproportionately lower relative to EPA. It can be concluded that other tissue biomarkers of LC n-3 PUFA status should be explored for their utility in assessing risk in patients undergoing haemodialysis treatment, for example, in platelets and mononuclear cells [55].

5. Conclusions

In summary, haemodialysis patients presented low proportions of erythrocyte membrane EPA + DHA. There were no relationships between erythrocyte membrane LC n-3 PUFA status and HRV, nor with estimated dietary intake of LC n-3 PUFA. However, negative associations between plasma EPA and erythrocyte EPA:DHA ratio, and 24 h parasympathetically driven parameters of HRV, suggest that further research is needed to investigate compromised conversion of EPA/DPA to DHA in patients that recently started haemodialysis as this may contribute to impaired cardiac autonomic function, conferring a greater risk of SCD.

Author Contributions

Conceptualization, W.L.H. and H.L.M.; methodology, W.L.H., H.L.M. and A.M.P.; formal analysis, W.L.H. and A.M.P.; investigation, A.M.P. and H.L.M.; writing—original draft preparation, W.L.H. and A.M.P.; writing—review and editing, W.L.H., A.M.P. and H.L.M.; supervision, W.L.H. and H.L.M.; project administration, A.M.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Committee of NRES Committee London–Camberwell St Giles (REC ref: 14/LO/0186, approved on the 14 March 2014).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

Robert Gray for technical support with the plasma and erythrocyte membrane fatty acid analysis. Renal consultants Nicola Kumar and Jasmine Lee for support with recruitment of patients.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Outline of study days.
Figure A1. Outline of study days.
Nutrients 13 02453 g0a1

Appendix B

Table A1. 24 h HRV measurements in the study participants with analysable records (n = 35).
Table A1. 24 h HRV measurements in the study participants with analysable records (n = 35).
24 h HRV ParametersMeanSD
HR (bpm)83.010.9
IBI (ms)747131
Ti
Median16.8
Lower and upper quartiles15.0, 27.1
SDNN (ms)84.125.6
SDANN (ms)76.225.6
pNN50 (%)4.84.0
RMSSD (ms)29.011.7
LF (ms2)
Median92.7
Lower and upper quartiles40.7, 194.5
HF (ms2)
Median128.7
Lower and upper quartiles68.4, 216.8
LF:HF
Median0.79
Lower and upper quartiles0.57, 1.45
LogLF:HF−0.220.79
VLF (ms2)48443218
SD1:SD2
Median0.15
Lower and upper quartiles0.12, 0.24
IBI, interbeat interval (also known as RR interval); HR, heart rate; bpm, beats per minute; HRV, heart rate variability; NN, normal-to-normal; Ti, triangular index; SDNN, standard deviation of all NN intervals; SDANN, standard deviation of the averaged NN intervals, calculated from 5 min epochs; RMSSD, the square root of the mean of the sum of squares of differences between adjacent NN intervals; PNN50, percentage of adjacent NN intervals that differed by >50 ms; LF, low-frequency power; HF, high-frequency power; VLF, very-low-frequency power; SD1:SD2, the ratio of the SD of beat-to-beat IBI variability (SD1) against the SD of long-term IBI variability (SD2).

References

  1. Foley, R.N.; Murray, A.M.; Li, S.; Herzog, C.A.; McBean, A.M.; Eggers, P.W.; Collins, A.J. Chronic Kidney Disease and the Risk for Cardiovascular Disease, Renal Replacement, and Death in the United States Medicare Population, 1998 to 1999. J. Am. Soc. Nephrol. 2004, 16, 489–495. [Google Scholar] [CrossRef]
  2. Herzog, C.A.; Mangrum, J.M.; Passman, R. Non-coronary heart disease in dialysis patients: Sudden cardiac death and dialysis patients. Semin. Dial. 2008, 21, 300–307. [Google Scholar] [CrossRef]
  3. Myerburg, R.J.; Junttila, M.J. Sudden Cardiac Death Caused by Coronary Heart Disease. Circulation 2012, 125, 1043–1052. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Pun, P.H.; Smarz, T.R.; Honeycutt, E.F.; Shaw, L.K.; Al-Khatib, S.M.; Middleton, J.P. Chronic kidney disease is associated with increased risk of sudden cardiac death among patients with coronary artery disease. Kidney Int. 2009, 76, 652–658. [Google Scholar] [CrossRef] [Green Version]
  5. Ranpuria, R.; Hall, M.; Chan, C.T.; Unruh, M. Heart rate variability (HRV) in kidney failure: Measurement and consequences of reduced HRV. Nephrol. Dial. Transplant. 2008, 23, 444–449. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Fukuta, H.; Hayano, J.; Ishihara, S.; Sakata, S.; Mukai, S.; Ohte, N.; Ojika, K.; Yagi, K.; Matsumoto, H.; Sohmiya, S.; et al. Prognostic value of heart rate variability in patients with end-stage renal disease on chronic haemodialysis. Nephrol. Dial. Transplant. 2003, 18, 318–325. [Google Scholar] [CrossRef]
  7. Hayano, J. Prognostic value of heart rate variability during long-term follow-up in chronic haemodialysis patients with end-stage renal disease. Nephrol. Dial. Transplant. 1999, 14, 1480–1488. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Oikawa, K.; Ishihara, R.; Maeda, T.; Yamaguchi, K.; Koike, A.; Kawaguchi, H.; Tabata, Y.; Murotani, N.; Itoh, H. Prognostic value of heart rate variability in patients with renal failure on hemodialysis. Int. J. Cardiol. 2009, 131, 370–377. [Google Scholar] [CrossRef]
  9. Madsen, T.; Christensen, J.H.; Svensson, M. Marine n-3 Polyunsaturated Fatty Acids in Patients With End-stage Renal Failure and in Subjects Without Kidney Disease: A Compararative Study. J. Ren. Nutr. 2011, 21, 169–175. [Google Scholar] [CrossRef]
  10. Shoji, T.; Kakiya, R.; Hayashi, T.; Tsujimoto, Y.; Sonoda, M.; Shima, H.; Mori, K.; Fukumoto, S.; Tahara, H.; Shioi, A.; et al. Serum n-3 and n-6 Polyunsaturated Fatty Acid Profile as an Independent Predictor of Cardiovascular Events in Hemodialysis Patients. Am. J. Kidney Dis. 2013, 62, 568–576. [Google Scholar] [CrossRef] [PubMed]
  11. Kutner, N.G.; Clow, P.W.; Zhang, R.; Aviles, X. Association of fish intake and survival in a cohort of incident dialysis patients. Am. J. Kidney Dis. 2002, 39, 1018–1024. [Google Scholar] [CrossRef] [PubMed]
  12. Friedman, A.N.; Saha, C.; Watkins, B.A. Feasibility Study of Erythrocyte Long-Chain Omega-3 Polyunsaturated Fatty Acid Content and Mortality Risk in Hemodialysis Patients. J. Ren. Nutr. 2008, 18, 509–512. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Calder, P.C. Marine omega-3 fatty acids and inflammatory processes: Effects, mechanisms and clinical relevance. Biochim. Biophys. Acta Mol. Cell Biol. Lipids 2015, 1851, 469–484. [Google Scholar] [CrossRef]
  14. Billman, G.E. The effects of omega-3 polyunsaturated fatty acids on cardiac rhythm: A critical reassessment. Parmacology Ther. 2013, 140, 53–80. [Google Scholar] [CrossRef] [PubMed]
  15. Giordano, E.; Visioli, F. Long-chain omega 3 fatty acids: Molecular bases of potential antioxidant actions. Prostaglandins Leukot. Essent. Fat. Acids 2014, 90, 1–4. [Google Scholar] [CrossRef]
  16. Christensen, J.H. Omega-3 Polyunsaturated Fatty Acids and Heart Rate Variability. Front. Physiol. 2011, 2, 84. [Google Scholar] [CrossRef] [Green Version]
  17. Christensen, J.H.; Aarøe, J.; Knudsen, N.; Dideriksen, K.; Kornerup, H.J.; Dyerberg, J.; Schmidt, E.B. Heart rate variability and n-3 fatty acids in patients with chronic renal failure-a pilot study. Clin. Nephrol. 1998, 49, 102–106. [Google Scholar]
  18. Svensson, M.; Schmidt, E.B.; Jørgensen, K.A.; Christensen, J.H. The effect of n-3 fatty acids on heart rate variability in patients treated with chronic hemodialysis. J. Ren. Nutr. 2007, 17, 243–249. [Google Scholar] [CrossRef]
  19. Malik, M.; Bigger, J.T.; Camm, A.J.; Kleiger, R.E.; Malliani, A.; Moss, A.J.; Schwartz, P.J. Heart rate variability: Standards of measurement, physiological interpretation and clinical use. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Eur. Heart J. 1996, 17, 354–381. [Google Scholar] [CrossRef] [Green Version]
  20. Lepage, G.; Roy, C.C. Direct transesterification of all classes of lipids in a one-step reaction. J. Lipid Res. 1986, 27, 114–120. [Google Scholar] [CrossRef]
  21. Hall, W.L.; Hay, G.; Maniou, Z.; Seed, P.T.; Chowienczyk, P.; Sanders, T. Effect of low doses of long chain n-3 polyunsaturated fatty acids on sleep-time heart rate variability: A randomized, controlled trial. Int. J. Cardiol. 2013, 168, 4439–4442. [Google Scholar] [CrossRef] [Green Version]
  22. Saran, R.; Li, Y.; Robinson, B.; Ayanian, J.; Balkrishnan, R.; Bragg-Gresham, J.; Chen, J.T.; Cope, E.; Gipson, D.; He, K.; et al. US Renal Data System 2014 Annual Data Report: Epidemiology of Kidney Disease in the United States. Am. J. Kidney Dis. 2015, 66, A7-305. [Google Scholar] [CrossRef] [Green Version]
  23. Timofte, D.; Tanasescu, M.-D.; Balan, D.G.; Tulin, A.; Stiru, O.; Vacaroiu, I.A.; Mihai, A.; Constantin, P.C.; Cosconel, C.-I.; Enyedi, M.; et al. Management of acute intradialytic cardiovascular complications: Updated overview (Review). Exp. Ther. Med. 2021, 21, 1. [Google Scholar] [CrossRef] [PubMed]
  24. Dessì, M.; Noce, A.; Bertucci, P.; Noce, G.; Rizza, S.; De Stefano, A.; di Villahermosa, S.M.; Bernardini, S.; de Lorenzo, A.; di Daniele, N. Plasma and erythrocyte membrane phospholipids and fatty acids in Italian general population and hemodialysis patients. Lipids Health Dis. 2014, 21, 13–54. [Google Scholar] [CrossRef] [Green Version]
  25. Friedman, A.N.; Yu, Z.; Tabbey, R.; Denski, C.; Tamez, H.; Wenger, J.; Thadhani, R.; Li, Y.; Watkins, B.A. Low Blood Levels of Long-Chain n–3 Polyunsaturated Fatty Acids in US Hemodialysis Patients: Clinical Implications. Am. J. Nephrol. 2012, 36, 451–458. [Google Scholar] [CrossRef] [Green Version]
  26. Ristić, V.; Tepsić, V.; Ristić-Medié, D.; Perunicić, G.; Rasić, Z.; Postić, M.; Arsić, A.; Blazencić-Mladenović, V.; Ristić, G. Plasma and erythrocyte phospholipid fatty acids composition in Serbian hemodialyzed patients. Ren. Fail. 2006, 28, 211–216. [Google Scholar] [CrossRef] [PubMed]
  27. Peck, W.; Monsen, R. Effect of three sources of long-chain fatty acids on the plasma fatty acid profile, plasma prostaglandin concentrations, and pruritus symptoms in hemodialysis. Am. J. Clin. Nutr. 1996, 64, 210–214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Gollasch, B.; Dogan, I.; Rothe, M.; Gollasch, M.; Luft, F.C. Effects of hemodialysis on blood fatty acids. Physiol. Rep. 2020, 8, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Friedman, A.N.; Moe, S.M.; Perkins, S.M.; Li, Y.; Watkins, B.A. Fish Consumption and Omega-3 Fatty Acid Status and Determinants in Long-Term Hemodialysis. Am. J. Kidney Dis. 2006, 47, 1064–1071. [Google Scholar] [CrossRef] [PubMed]
  30. Harris, W.S.; Sands, S.A.; Windsor, S.L.; Ali, H.A.; Stevens, T.L.; Magalski, A.; Porter, C.B.; Borkon, A.M. Omega-3 fatty acids in cardiac biopsies from heart transplantation patients: Correlation with erythrocytes and response to supplementation. Circulation 2004, 110, 1645–1649. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Harris, W.S.; von Schacky, C. The Omega-3 Index: A new risk factor for death from coronary heart disease? Prev. Med. 2004, 39, 212–220. [Google Scholar] [CrossRef]
  32. Babitt, J.L.; Lin, H.Y. Mechanisms of Anemia in CKD. J. Am. Soc. Nephrol. 2012, 23, 1631–1634. [Google Scholar] [CrossRef] [Green Version]
  33. Serra-Majem, L.; Nissensohn, M.; Øverby, N.C.; Fekete, K. Dietary methods and biomarkers of omega 3 fatty acids: A systematic review. Br. J. Nutr. 2012, 107 (Suppl. S2), S64–S76. [Google Scholar] [CrossRef] [Green Version]
  34. Rantanen, J.M.; Schmidt, E.B.; Riahi, S.; Lundbye-Christensen, S.; Christensen, J.H. Marine n-3 PUFA, heart rate variability and ventricular arrhythmias in patients on chronic dialysis: A cross-sectional study. Br. J. Nutr. 2018, 120, 317–325. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Hall, W.L. The future for long chain n-3 PUFA in the prevention of coronary heart disease: Do we need to target non-fish-eaters? Proc. Nutr. Soc. 2017, 76, 408–418. [Google Scholar] [CrossRef]
  36. Steinberg, A.A.; Mars, R.L.; Goldman, D.S.; Percy, R.F. Effect of end-stage renal disease on decreased heart rate variability. Am. J. Cardiol. 1998, 82, 1156–1158. [Google Scholar] [CrossRef]
  37. Tamura, K.; Tsuji, H.; Nishiue, T.; Yajima, I.; Higashi, T.; Iwasaka, T. Determinants of heart rate variability in chronic hemodialysis patients. Am. J. Kidney Dis. 1998, 31, 602–606. [Google Scholar] [CrossRef]
  38. Mann, M.C.; Exner, D.V.; Hemmelgarn, B.R.; Hanley, D.A.; Turin, T.C.; MacRae, J.M.; Ahmed, S.B. The VITAH Trial—Vitamin D Supplementation and Cardiac Autonomic Tone in Patients with End-Stage Kidney Disease on Haemodialysis: A Blinded, Randomized Controlled Trial. Nutrients 2016, 8, 608. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Agarwal, A.; Anand, I.S.; Sakhuja, V.; Chugh, K.S. Effect of dialysis and renal transplantation on autonomic dysfunction in chronic renal failure. Kidney Int. 1991, 40, 489–495. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Tong, Y.-Q.; Hou, H.-M. Alteration of Heart Rate Variability Parameters in Nondiabetic Hemodialysis Patients. Am. J. Nephrol. 2007, 27, 63–69. [Google Scholar] [CrossRef] [PubMed]
  41. Poulikakos, D.; Malik, M.; Banerjee, D. Sex-Dependent Association between Heart Rate Variability and Pulse Pressure in Haemodialysis Patients. Nephron Clin. Pr. 2014, 128, 361–366. [Google Scholar] [CrossRef]
  42. Coquet, I.; Mousson, C.; Rifle, G.; Laurent, G.; Moreau, D.; Cottin, Y.; Zeller, M.; Touzery, C.; Wolf, J.E. Influence of Ischemia on Heart-Rate Variability in Chronic Hemodialysis Patients. Ren. Fail. 2005, 27, 7–12. [Google Scholar] [CrossRef] [PubMed]
  43. Chan, C.T.; Levin, N.W.; Chertow, G.M.; Larive, B.; Schulman, G.; Kotanko, P. Determinants of cardiac autonomic dysfunction in ESRD. Clin. J. Am. Soc. Nephrol. 2010, 5, 1821–1827. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Sandrone, G.; Mortara, A.; Torzillo, D.; La Rovere, M.T.; Malliani, A.; Lombardi, F. Effects of beta blockers (atenolol or metoprolol) on heart rate variability after acute myocardial infarction. Am. J. Cardiol. 1994, 74, 340–345. [Google Scholar] [CrossRef]
  45. Wennerblom, B.; Lurje, L.; Karlsson, T.; Tygesen, H.; Vahisalo, R.; Hjalmarson, A. Circadian variation of heart rate variability and the rate of autonomic change in the morning hours in healthy subjects and angina patients. Int. J. Cardiol. 2001, 79, 61–69. [Google Scholar] [CrossRef]
  46. Tory, K.; Horváth, E.; Süveges, Z.; Fekete, A.; Sallay, P.; Berta, K.; Szabó, T.; Szabó, A.J.; Tulassay, T.; Reusz, G.S. Effect of propranolol on heart rate variability in patients with end-stage renal disease: A crossover pilot trial. Clin. Nephrol. 2004, 61, 316–323. [Google Scholar] [CrossRef]
  47. Yong, K.; Dogra, G.; Boudville, N.C.; Lim, W.H. Increased Inflammatory Response in Association with the Initiation of Hemodialysis Compared with Peritoneal Dialysis in a Prospective Study of End-Stage Kidney Disease Patients. Perit. Dial. Int. 2018, 38, 18–23. [Google Scholar] [CrossRef] [PubMed]
  48. Bazeley, J.; Bieber, B.; Li, Y.; Morgenstern, H.; De Sequera, P.; Combe, C.; Yamamoto, H.; Gallagher, M.; Port, F.K.; Robinson, B.M. C-Reactive Protein and Prediction of 1-Year Mortality in Prevalent Hemodialysis Patients. Clin. J. Am. Soc. Nephrol. 2011, 6, 2452–2461. [Google Scholar] [CrossRef] [Green Version]
  49. Psychari, S.N.; Sinos, L.; Iatrou, C.; Liakos, G.; Apostolou, T.S. Relations of inflammatory markers to lipid levels and autonomic tone in patients with moderate and severe chronic kidney disease and in patients under maintenance hemodialysis. Clin. Nephrol. 2005, 64, 419–427. [Google Scholar] [CrossRef]
  50. Von Känel, R.; Carney, R.M.; Zhao, S.; Whooley, M.A. Heart rate variability and biomarkers of systemic inflammation in patients with stable coronary heart disease: Findings from the Heart and Soul Study. Clin. Res. Cardiol. 2011, 100, 241–247. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Danesh, J.; Wheeler, J.G.; Hirschfield, G.; Eda, S.; Eiriksdottir, G.; Rumley, A.; Lowe, G.D.; Pepys, M.B.; Gudnason, V. C-Reactive Protein and Other Circulating Markers of Inflammation in the Prediction of Coronary Heart Disease. N. Engl. J. Med. 2004, 350, 1387–1397. [Google Scholar] [CrossRef]
  52. Poppitt, S.D.; Kilmartin, P.; Butler, P.; Keogh, G.F. Assessment of erythrocyte phospholipid fatty acid composition as a biomarker for dietary MUFA, PUFA or saturated fatty acid intake in a controlled cross-over intervention trial. Lipids Health Dis. 2005, 4, 30. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Zong, G.; Zhu, J.; Sun, L.; Ye, X.; Lu, L.; Jin, Q.; Zheng, H.; Yu, Z.; Zhu, Z.; Li, H.; et al. Associations of erythrocyte fatty acids in the de novo lipogenesis pathway with risk of metabolic syndrome in a cohort study of middle-aged and older Chinese. Am. J. Clin. Nutr. 2013, 98, 319–326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Lemaitre, R.N.; King, I.B.; Sotoodehnia, N.; Knopp, R.H.; Mozaffarian, D.; McKnight, B.; Rea, T.D.; Rice, K.; Friedlander, Y.; Lumley, T.S.; et al. Endogenous red blood cell membrane fatty acids and sudden cardiac arrest. Metabolism 2010, 59, 1029–1034. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Browning, L.M.; Walker, C.G.; Mander, A.P.; West, A.L.; Madden, J.; Gambell, J.M.; Young, S.; Wang, L.; Jebb, S.A.; Calder, P.C. Incorporation of eicosapentaenoic and docosahexaenoic acids into lipid pools when given as supplements providing doses equivalent to typical intakes of oily fish. Am. J. Clin. Nutr. 2012, 96, 748–758. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Study flow chart.
Figure 1. Study flow chart.
Nutrients 13 02453 g001
Table 1. Characteristics of the study population (n = 45).
Table 1. Characteristics of the study population (n = 45).
N%
Sex
Male2556
Female2044
Ethnicity
White2044
Black1942
Other613
Co-morbidities
Hypertension4091
Coronary artery disease511
Peripheral vascular disease25
Heart failure37
Cerebrovascular disease37
Type 2 Diabetes1739
Hyperlipidaemia921
Anaemia511
Medications
ACE inhibitor1125
ARB1023
β-blocker2046
α-blocker1227
Diuretic1023
Statin2148
Insulin1125
Anti-diabetic614
Anti-depressant614
Vitamin D2557
Calcium-based phosphate binders1841
Iron injections2148
Erythropoietin therapy2045
Subjective daytime sleepiness
Normal2374
Mild516
Moderate13
Severe27
Risk of sleep apnoea
High2065
Low1135
MeanSD
Age (years)589
Waist Circumference (cm)
Male105.613.2
Female97.119.4
BMI (kg/m2)27.87.0
Average SBP (mmHg) 114116
Average DBP (mmHg) 1769
Fluid removal (L)2.00.8
Pre- and post-dialysis measurements
Pre-dialysis creatinine (μmol/L)585192
Post-dialysis creatinine (μmol/L)228100
Pre-dialysis albumin (g/L)
Median38
Lower and upper quartiles36, 40
Post-dialysis albumin (g/L)
Median40
Lower and upper quartiles37, 44
Pre-dialysis sodium (mmol/L)
Median139.3
Lower and upper quartiles137.6, 141.0
Post-dialysis sodium (mmol/L)
Median138.0
Lower and upper quartiles137.0, 139.0
Pre-dialysis calcium (mmol/L)2.170.15
Post-dialysis calcium (mmol/L)2.250.19
Pre-dialysis potassium (mmol/L)
Geometric Mean5.1
95% CI4.9, 5.4
Post-dialysis potassium (mmol/L)
Geometric Mean3.8
95% CI3.6, 3.9
Dialysis 2 h HR and HRV 2
HR (bpm)77.215.0
IBI (ms)787132
Ti11.13.6
SDNN (ms)45.813.9
pNN50 (%)
Median1.86
Lower and upper quartilesRMSSD (ms)0.68, 5.41
Median21.4
Lower and upper quartiles17.4, 36.5
LF (ms2)
Median101
Lower and upper quartiles46, 197
HF (ms2)
Median114
Lower and upper quartiles55, 195
logLF:HF−0.020.43
SD1:SD2
Median0.24
Lower and upper quartiles0.18, 0.42
1 12 months prior to commencing haemodialysis. ACE, angiotensin converting enzyme; ARB, angiotensin II receptor antagonists; BMI, body mass index; DBP, diastolic blood pressure; SBP, systolic blood pressure. Results expressed as percentage and numbers; mean values and standard deviations (SD); medians and interquartile rage (IQR); or geometric mean and 95% confidence intervals (CI). 2 Average of V1 and V2 n = 23, V1 or V2 n = 39, missing HRV values n = 6.
Table 2. Comparison of LC n-3 PUFA status and intake, sleep-time HRV and inflammatory marker in haemodialysis patients and age-, sex- and BMI-matched healthy controls from the MARINA study.
Table 2. Comparison of LC n-3 PUFA status and intake, sleep-time HRV and inflammatory marker in haemodialysis patients and age-, sex- and BMI-matched healthy controls from the MARINA study.
Haemodialysis Patients (n = 42)Healthy Reference (n = 42)p-Value
Sex (M:F)22:2022:20
Age (y)57 ± 857 ± 8
Ethnicity (White:Black:Other)19:18:532:7:30.013 §
BMI (kg/m2)27.6 ± 7.026.1 ± 3.40.231
hs-CRP (mg/L) a5.0 (3.3, 7.6)0.6 (0.4, 0.9)<0.001
EPA dietary intake (mg/day)0.09 (0.06, 0.14)0.10 (0.07, 0.13)0.732
DHA dietary intake (mg/day)0.27 (0.16, 0.45)0.18 (0.13, 0.25)0.176
Erythrocyte EPA (%) a0.6 (0.54, 0.77)1.2 (1.09, 1.35)<0.001
Erythrocyte DPA (%) a1.5 (1.37, 1.63)3.3 (3.1, 3.5)<0.001
Erythrocyte DHA (%)2.6 ± 0.946.6 ± 1.42<0.001
Omega-3 index (%)3.4 ± 1.177.9 ± 1.67<0.001
Erythrocyte DPAn-6 (%) a0.23 (0.17, 0.42)0.45 (0.38, 0.52)<0.001
Erythrocyte AA (%)9.9 ± 1.2917.2 ± 1.73<0.001
Erythrocyte DGLA (%) b1.10 (0.88, 1.23) 1.81 (1.63, 2.07)<0.001
Erythrocyte LA (%)7.6 ± 1.2811.1 ± 1.61<0.001
Erythrocyte cis-VA (%) b1.41 (1.36, 1.72)1.23 (1.11, 1.38)<0.001
Erythrocyte OA (%)15.7 ± 1.714.9 ± 1.04<0.01
Erythrocyte SA (%)23.1 ± 1.5915.7 ± 1.71<0.001
Erythrocyte POA (%) a0.42 (0.32, 0.53)0.27 (0.23, 0.32)<0.005
Erythrocyte PA (%) b30.4 (29.0, 31.9)20.0 (18.8, 20.7)<0.001
Erythrocyte EPA:DHA ratio b0.31 (0.22, 0.39)0.18 (0.16, 0.22)<0.001
Erythrocyte PA:LA ratio b4.00 (3.48, 4.47)1.80 (1.65, 1.95)<0.001
Erythrocyte AA:EPA ratio b14.0 (10.4, 19.3)14.2 (10.7, 17.5)0.741
Plasma EPA (%) a0.8 (0.63, 0.91)1.0 (0.91, 1.18)0.007
Plasma DPA (%) a0.44 (0.40, 0.47)0.64 (0.60, 0.68)<0.001
Plasma DHA (%)1.8 ± 0.652.4 ± 0.63<0.001
HR and HRV parameters (Sleep time)N = 40N = 40
HR (bpm)77.4 ± 12.164.6 ± 9.7<0.001
IBI (ms)800 ± 133962 ± 115<0.001
Ti12.2 ± 3.822.6 ± 6.3<0.001
SDNN (ms)53.1 ± 16.794.3 ± 25.9<0.001
pNN50 (%) a2.6 (1.7, 4.1)6.2 (4.1, 9.5)0.006
RMSSD (ms) a23.5 (19.9, 27.8)34.4 (29.4, 40.4)0.001
LF (ms2) a104 (79, 137)655 (504, 852)<0.001
HF (ms2) a106 (76, 148)316 (231, 432)<0.001
logLF:HF−0.01 ± 0.822.64 ± 1.91<0.001
VLF (ms2) a1514 (1257, 1824)4932 (4132, 5888)<0.001
SD1:SD2 a0.24 (0.21, 0.29)0.20 (0.17, 0.23)0.073
Results expressed as a ratio or mean ± SD, except a geometric means (95%CI) and p-value obtained by independent t-test, except § χ2 test, and b median (lower and upper quartiles) and p-value obtained by independent samples Mann-Whitney U-test. AA, arachidonic acid (20:4n-6); BMI, body mass index; DBP, diastolic blood pressure; DGLA, dihomo-γ-linolenic acid(20:3n-6); DHA, docosahexaenoic acid; DPA, docosapentaenoic acid; EPA, eicosapentaenoic acid; HR, heart rate; MA, myristic acid (14:0); OA, oleic acid (18:1n-9); PA, palmitic acid (16:0); POA, palmitoleic acid (16:1n-7); PUFA, polyunsaturated fatty acids; SA, stearic acid (18:0); SBP, systolic blood pressure; cis-VA, cis-vaccenic acid (18:1n-7); IBI, interbeat interval (also known as RR interval), the time interval between R spikes of the QRS complex of the electrocardiogram; HR, heart rate; bpm, beats per minute; HRV, heart rate variability; NN, normal-to-normal; Ti, triangular index (total number of all NN intervals divided by the height of the histogram of all NN intervals); SDNN, standard deviation of all NN intervals (NN intervals, similar to R-R, but on normalised IBI data); SDANN, standard deviation of the averaged NN intervals, calculated from 5min epochs; RMSSD, the square root of the mean of the sum of squares of differences between adjacent NN intervals; PNN50, percentage of adjacent NN intervals that differed by >50 ms; LF, low-frequency power; HF, high-frequency power; VLF, very-low-frequency power; SD1:SD2, the ratio of the SD of beat-to-beat IBI variability (SD1) against the SD of long-term IBI variability (SD2).
Table 3. Simple and multivariable linear regression analyses between erythrocyte membrane LC n-3 PUFA and HRV during dialysis (2 h), sleep time and over a 24 h period.
Table 3. Simple and multivariable linear regression analyses between erythrocyte membrane LC n-3 PUFA and HRV during dialysis (2 h), sleep time and over a 24 h period.
EPA (wt %) aDHA (wt %)Omega-3 Index
β1p1β2p2β1p1β2p2β1p1β2p2
Dialysis (n = 36)
 HR (bpm)−0.3790.017−0.2530.051−0.2510.123−0.1340.298−0.3440.032−0.1980.124
 IBI (ms)0.1900.2480.8800.5630.1870.2530.1530.3020.2300.1600.1700.255
 Ti0.1650.3160.0470.7590.2310.1570.2800.0560.2630.1060.2790.060
 SDNN (ms)0.1360.4080.0230.8820.1010.5420.1420.3470.1330.4190.1410.354
 pNN50 (%) a0.0530.7470.0120.9400.0900.5840.1270.4300.1000.5460.1240.445
 RMSSD (ms) a−0.0990.549−0.1430.3930.0420.7990.0880.5920.0110.9490.0430.794
 LF (ms2) a0.2250.1690.0880.5280.0950.5650.1510.2630.1540.3500.1670.217
 HF (ms2) a−0.0010.994−0.0740.6270.1290.4330.1950.1860.1200.4660.1640.271
 logLF:HF0.2500.1240.1950.230−0.0920.579−0.1010.528−0.0130.937−0.0410.800
 SDANN (ms)0.1140.4890.0460.7910.1430.3870.1620.3390.1660.3130.1660.329
 SD1:SD2 a−0.2010.220−0.1790.305−0.0290.863−0.0020.993−0.0850.608−0.0520.766
Sleep time (n = 39)
 HR (bpm)−0.3060.062−0.2740.110−0.0860.608−0.0340.840−0.1690.304−0.1100.522
 IBI (ms)0.3090.0590.2700.1140.0790.6370.0320.8500.1540.3490.1060.532
 Ti0.1840.2680.1280.469−0.0180.915−0.0660.7010.0370.825−0.0270.879
 SDNN (ms)0.2200.1850.1560.378−0.0350.836−0.0660.7030.0320.850−0.0180.916
 pNN50 (%) a−0.0430.798−0.1350.4290.1800.2800.1620.3410.1140.4960.0330.846
 RMSSD (ms) a−0.0430.797−0.1270.4660.1730.3000.1320.4350.1480.3760.0890.602
 LF (ms2) a0.1560.3490.0180.9090.1640.3270.1570.3090.1970.2350.1530.323
 HF (ms2) a0.0960.568−0.0130.9380.2300.1660.2020.2070.2410.1450.1870.247
 logLF:HF0.0550.7430.0220.902−0.0990.556−0.0730.674−0.0760.652−0.0630.719
 VLF (ms2) a0.2680.1040.1940.266−0.1010.548−0.1370.4210.0180.913−0.0350.844
 SDANN (ms) a0.1450.3770.2030.221−0.1270.446−0.1610.361−0.0590.725−0.0990.580
 SD1:SD2 a−0.2000.229−0.2390.1770.2140.1970.1930.2670.1410.3990.1150.515
24 h (n = 35)
 HR (bpm)−0.2550.145−0.1520.451−0.2160.220−0.1460.430−0.2790.110−0.2000.308
 IBI (ms)0.3150.0690.2130.2820.2170.2170.1480.4140.2970.0880.2190.256
 Ti a0.4300.0110.3340.094−0.0080.965−0.0780.6770.1090.5380.0010.998
 SDNN (ms)0.3920.0220.2100.289−0.0350.842−0.1730.3420.0730.684−0.1360.484
 pNN50 (%)−0.0510.775−0.1450.5080.1640.3550.1860.3510.1450.4140.1680.432
 RMSSD (ms)−0.2300.190−0.3300.1300.1910.2790.2560.2010.1230.4890.1970.359
 LF (ms2) a0.1730.329−0.1750.2960.2650.1300.2210.1450.3040.0810.1980.222
 HF (ms2) a−0.0900.611−0.3110.1120.2860.1020.3230.0680.2520.1500.2760.149
 logLF:HF0.3000.0850.1270.5030.0330.851−0.0470.7860.1140.521−0.0190.918
 VLF (ms2)0.4120.0160.2930.150−0.0500.777−0.1680.3730.0630.722−0.1090.588
 SDANN (ms)0.3950.0210.2390.239−0.0670.707−0.2110.2550.0430.810−0.1710.390
 SD1:SD2 a−0.3490.043−0.3100.1460.1960.2670.3160.1020.0950.5940.2690.197
β1unadjusted standardised beta coefficient; p1—unadjusted p-value; β2adjusted standardised beta coefficient; p2 value, adjusted for age, BMI and β-blockers (yes or no). 24 h HRV analysis has activity level from accelerometery data (counts per minute, cpm) as an additional adjustment for p2. VLF during dialysis not presented due to shorter duration of recording (2 h). a Regression analysis based on LN variables. IBI, interbeat interval (also known as RR interval), the time interval between R spikes of the QRS complex of the electrocardiogram; HR, heart rate; bpm, beats per minute; HRV, heart rate variability; NN, normal-to-normal; Ti, triangular index (total number of all NN intervals divided by the height of the histogram of all NN intervals); SDNN, standard deviation of all NN intervals (NN intervals, similar to R-R, but on normalised IBI data); SDANN, standard deviation of the averaged NN intervals, calculated from 5 min epochs; RMSSD, the square root of the mean of the sum of squares of differences between adjacent NN intervals; PNN50, percentage of adjacent NN intervals that differed by >50 ms; LF, low-frequency power; HF, high-frequency power; VLF, very-low-frequency power; SD1:SD2, the ratio of the SD of beat-to-beat IBI variability (SD1) against the SD of long-term IBI variability (SD2).
Table 4. Simple and multivariable linear regression analyses between plasma LC n-3 PUFA and HRV during dialysis (2 h), sleep time and over a 24 h period.
Table 4. Simple and multivariable linear regression analyses between plasma LC n-3 PUFA and HRV during dialysis (2 h), sleep time and over a 24 h period.
EPA aDPA(n-3) aDHA
β1p1β2p2β1p1β2p2β1p1β2p2
Dialysis (n = 39)
 HR (bpm)−0.1290.435−0.1190.3470.2070.2070.1450.253−0.2980.065−0.2080.101
 IBI (ms)0.1510.3590.1590.273−0.1780.279−0.1110.4480.2390.1430.2230.128
 Ti0.1430.3860.1340.362−0.0510.759−0.0470.7520.3280.0420.3310.022
 SDNN (ms)0.1230.4540.1180.422−0.1000.546−0.0870.5580.2520.1220.2580.081
 pNN50 (%) a−0.0670.686−0.0520.743−0.1590.334−0.1180.4590.0560.7370.0920.567
 RMSSD (ms)−0.1960.231−0.1840.249−0.2790.085−0.2530.112−0.0230.8920.0150.925
 LF (ms2) a0.0950.5640.0810.539−0.1470.372−0.1530.2480.2390.1420.2350.076
 HF (ms2) a−0.0860.602−0.0730.617−0.2210.176−0.1850.2020.1200.4680.1680.254
 logLF:HF0.1960.2310.1660.2850.0950.5660.0420.7910.1000.5460.0440.781
 SDANN (ms)0.1820.2660.1770.282−0.0130.939−0.0070.9650.2700.0960.2700.105
 SD1:SD2 a−0.3080.057−0.2910.076−0.2310.158−0.2110.209−0.2020.218−0.1630.337
Sleep time (n = 40)
 HR (bpm)−0.2630.106−0.2790.0800.1310.4270.0800.626−0.2250.168−0.2050.210
 IBI (ms)0.2800.0850.2960.060−0.1000.543−0.0460.7800.2180.1830.2010.218
 Ti0.1910.2440.1860.264−0.1460.375−0.1330.4310.1270.4410.0900.599
 SDNN (ms)0.2540.1190.2460.1370.0400.8080.0520.7580.0460.7810.0060.970
 pNN50 (%) a−0.1800.274−0.1740.272−0.3620.024−0.3210.041−0.1130.495−0.1640.309
 RMSSD (ms)−0.1830.264−0.1760.277−0.3310.040−0.2900.071−0.1160.482−0.1530.353
 LF (ms2) a0.0310.8510.0250.869−0.1420.390−0.1100.4790.0680.6800.0350.823
 HF (ms2) a−0.0480.773−0.0420.789−0.2220.174−0.1730.268−0.0150.929−0.0470.765
 logLF:HF0.0880.5930.0750.6580.0980.5530.0790.6440.0860.6040.0870.616
 VLF (ms2) a0.2850.0790.2770.0940.0780.6390.0850.6100.0990.5510.0530.754
 SDANN (ms) a0.3200.0470.3080.0640.1810.2700.1720.3160.0390.8120.0040.980
 SD1:SD2 a−0.3540.027−0.3400.036−0.3640.023−0.3300.045−0.1590.333−0.1720.311
24 h (n = 35)
 HR (bpm)−0.2800.109−0.2230.1950.1960.2660.1540.369−0.2980.087−0.2560.170
 IBI (ms)0.3520.0410.2960.077−0.1440.416−0.0970.5680.3350.0530.2970.103
 Ti a0.3470.0440.2810.112−0.2220.206−0.2110.2210.2820.1070.2210.244
 SDNN (ms)0.4410.0090.3780.022−0.0520.769−0.0260.8790.2540.1470.1240.506
 pNN50 (%)−0.3020.083−0.3310.073−0.3940.021−0.3760.038−0.0750.675−0.0970.636
 RMSSD (ms)−0.3370.052−0.3530.057−0.3870.024−0.3860.034−0.1200.498−0.1170.573
 LF (ms2) a0.0620.729−0.0170.909−0.1340.449−0.1080.4510.2190.2130.0890.574
 HF (ms2) a−0.2110.232−0.2590.125−0.3430.047−0.3100.0610.0310.8630.0030.986
 logLF:HF0.2980.0870.2580.1060.2100.2340.2050.1990.2400.1710.1250.481
 SDANN (ms)0.4640.0060.4050.016−0.0160.929−0.0100.9540.2470.1590.1210.525
 VLF (ms2)0.4810.0040.4340.010−0.0380.829−0.0100.9560.2490.1550.1550.422
 SD1:SD2 a−0.4820.004−0.4610.009−0.3310.056−0.3370.060−0.2270.197−0.1520.453
β1—unadjusted standardised beta coefficient; p1—unadjusted p-value; β2—adjusted standardised beta coefficient; p2 value, adjusted for age, BMI and β-blockers (yes or no). 24 h HRV analysis has activity level from accelerometery data (counts per minute, cpm) as an additional adjustment for p2. VLF during dialysis not presented due to shorter duration of recording (2 h). a Regression analysis based on LN variables. IBI, interbeat interval (also known as RR interval), the time interval between R spikes of the QRS complex of the electrocardiogram; HR, heart rate; bpm, beats per minute; HRV, heart rate variability; NN, normal-to-normal; Ti, triangular index (total number of all NN intervals divided by the height of the histogram of all NN intervals); SDNN, standard deviation of all NN intervals (NN intervals, similar to R-R, but on normalised IBI data); SDANN, standard deviation of the averaged NN intervals, calculated from 5 min epochs; RMSSD, the square root of the mean of the sum of squares of differences between adjacent NN intervals; PNN50, percentage of adjacent NN intervals that differed by >50 ms; LF, low-frequency power; HF, high-frequency power; VLF, very-low-frequency power; SD1:SD2, the ratio of the SD of beat-to-beat IBI variability (SD1) against the SD of long-term IBI variability (SD2).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pinto, A.M.; MacLaughlin, H.L.; Hall, W.L. Heart Rate Variability and Long Chain n-3 Polyunsaturated Fatty Acids in Chronic Kidney Disease Patients on Haemodialysis: A Cross-Sectional Pilot Study. Nutrients 2021, 13, 2453. https://0-doi-org.brum.beds.ac.uk/10.3390/nu13072453

AMA Style

Pinto AM, MacLaughlin HL, Hall WL. Heart Rate Variability and Long Chain n-3 Polyunsaturated Fatty Acids in Chronic Kidney Disease Patients on Haemodialysis: A Cross-Sectional Pilot Study. Nutrients. 2021; 13(7):2453. https://0-doi-org.brum.beds.ac.uk/10.3390/nu13072453

Chicago/Turabian Style

Pinto, Ana M, Helen L MacLaughlin, and Wendy L Hall. 2021. "Heart Rate Variability and Long Chain n-3 Polyunsaturated Fatty Acids in Chronic Kidney Disease Patients on Haemodialysis: A Cross-Sectional Pilot Study" Nutrients 13, no. 7: 2453. https://0-doi-org.brum.beds.ac.uk/10.3390/nu13072453

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop