Two-loop QED corrections to the Altarelli-Parisi splitting functions

We compute the two-loop QED corrections to the Altarelli-Parisi (AP) splitting functions by using a deconstructive algorithmic Abelianization of the well-known NLO QCD corrections. We present explicit results for the full set of splitting kernels in a basis that includes the leptonic distribution functions that, starting from this order in the QED coupling, couple to the partonic densities. Finally, we perform a phenomenological analysis of the impact of these corrections in the splitting functions.