We consider the transverse-momentum (qT) distribution of Z bosons produced in hadronic collisions. At small values of qT, we perform the analytic resummation of the logarithmically enhanced QED contributions up to next-to-leading logarithmic accuracy, including the mixed QCD-QED contributions at leading logarithmic accuracy. Resummed results are consistently matched with the next-to-leading fixed-order results at small, intermediate and large values of qT and combined the known QCD results. We show numerical results at LHC and Tevatron energies, studying the impact of the QED corrections and providing an estimate of the corresponding perturbative uncertainty. We show preliminary results about the extension of our analytic results, valid for the production of generic neutral and colourless high-mass systems in hadronic collision, to the case of charged W boson production.