In this work, we study the very low frequency (VLF) modulation (range 0.01-0.03 Hz) in QRS slopes, heart rate variability (HRV) and ECG-derived respiration in hemodialysis patients. First, the relation between QRS slopes and HRV in the VLF band is measured using ordinary coherence. Then, partial coherence is used to measure the former relationship once the effect related to respiration is removed. Ordinary coherence values above a statistical threshold revealed linear relationship between VLF modulation in QRS slopes and HRV in about 10% of analyzed segments, with mean ± SD values of 0.79 ± 0.07 for upward slope and 0.77 ± 0.06 for downward slope. For these segments, partial coherence values drop below the threshold for 64% of the cases for upward slope and 76% for downward slope, suggesting that the origin of the VLF modulation in QRS slopes is mainly driven by respiration or linearly related to it. In the rest of the cases, partial coherence values dropped with respect to ordinary coherence from 0.89 to 0.77 for upward slope and from 0.86 to 0.75 for downward slope, suggesting that other ANS effects non-linearly related to respiration also contribute to the VLF modulation in QRS slopes.