Two-loop operator matrix elements for massive fermionic local twist-2 operators in QED