We present the closed analytic expression of the form factors of the
 two-loop QED vertex amplitude for on-shell electrons of finite mass $m$
 and arbitrary momentum transfer $S=-Q^2$. The calculation is carried out
 within the continuous $D$-dimensional regularization scheme, with a single
 continuous parameter $D$, the dimension of the space-time, which
 regularizes at the same time UltraViolet (UV) and InfraRed (IR)
 divergences. The results are expressed in terms of 1-dimensional harmonic
 polylogarithms of maximum weight 4.