  When a TeV-scale leptoquark has a sizeable Yukawa coupling, its dominant
production mechanism at hadron colliders is the partonic-level lepton-quark
fusion. Even though the parton distribution functions for leptons inside the
proton are minuscule, they get compensated by the resonant enhancement. We
present the first computation of higher order radiative corrections to the
resonant leptoquark production cross section at the Large Hadron Collider
(LHC). Next-to-leading (NLO) QCD and QED corrections are similar in size but
come with the opposite sign. We compute NLO $K$-factors for a wide range of
scalar leptoquark masses, as well as, all possible combinations of quark and
lepton flavors and leptoquark charges. Theoretical uncertainties due to the
renormalisation and factorisation scale variations and the limited knowledge of
parton distribution functions are quantified. We finally discuss how to
disentangle the flavor structure of leptoquark interactions by exploiting the
interplay between different production channels.
