The ratios among the leading-order (LO) hadronic vacuum polarization (HVP) contributions to the anomalous magnetic moments of electron, muon and τ -lepton, a HVP,LO =e,µ,τ , are computed using lattice QCD+QED simulations. The results include the effects at order O(α 2 em ) as well as the electromagnetic and strong isospin-breaking corrections at orders O(α 3 em ) and O(α 2 em (m u − m d )), respectively, where (m u − m d ) is the u-and d-quark mass difference. We employ the gauge configurations generated by the Extended Twisted Mass Collaboration with N f = 2 + 1 + 1 dynamical quarks at three values of the lattice spacing (a 0.062, 0.082, 0.089 fm) with pion masses in the range 210−450 MeV. The calculations are based on the quark-connected contributions to the HVP in the quenched-QED approximation, which neglects the charges of the sea quarks. The quark-disconnected terms are estimated from results available in the literature. We show that in the case of the electronmuon ratio the hadronic uncertainties in the numerator and in the denominator largely cancel out, while in the cases of the electron-τ and muon-τ ratios such a cancellation does not occur. For the electron-muon ratio we get R e/µ ≡ (m µ /m e ) 2 (a HVP,LO e /a HVP,LO µ ) = 1.1478 (70) with an uncertainty of 0.6%. Our result, which represents an accurate Standard Model (SM) prediction, agrees very well with the estimate obtained using the results of dispersive analyses of the experimental e + e − → hadrons data. Instead, it differs by 2.7 standard deviations from the value expected from present electron and muon (g − 2) experiments