A novel method for QRST cancellation during atrial fibrillation (AF) is introduced for use in recordings with two or more leads. The method is based on an echo state neural network (ESN) which estimates the time-varying, nonlinear transfer function between two leads, one lead with atrial activity and another lead without, for the purpose of canceling ventricular activity. The performance is evaluated on ECG signals, with simulated f-waves of low amplitude added, by determining the root mean square error P between the true f-wave signal and the estimated signal, as well as by evaluating the dominant AF frequency. When compared to average beat subtraction (ABS), being the most widely used method for QRST cancellation, the performance is found to be significantly better with equal to mean and standard deviation of PESN 24.8±7.3 and PABS 34.2±17.9 μV (p