  In this paper we discuss the effect of the complete leading-order QED
corrections to the DGLAP equations in the perturbative evolution of parton
distribution functions (PDFs). This requires the extension of the purely QCD
DGLAP evolution, including a PDF for the photons and, consistently, also for
the charged leptons $e^{\pm}$, $\mu^\pm$ and $\tau^\pm$. We present the
implementation of the QED-corrected DGLAP evolution in the presence of photon
and lepton PDFs in the APFEL program and, by means of different assumptions for
the initial scale PDFs, we produce for the first time PDF sets containing
charged lepton distributions. We also present phenomenological studies that aim
to assess the impact of the presence of lepton PDFs in the proton for some
relevant SM (and BSM) processes at the LHC at 13 TeV and the FCC-hh at 100 TeV.
The impact of the photon PDF is also outlined for those processes.
