We overview the general status of higher order corrections to Bhabha scattering and review recent progress in the determination of the two-loop virtual corrections. Quite recently, they were derived from combining a massless calculation and contributions with electron sub-loops. For a massive calculation, the self-energy and vertex master integrals are known, while most of the two-loop boxes are not. We demonstrate with an example that a study of systems of differential equations, combined with Mellin-Barnes representations for single masters, might open a way for their systematic calculation.