We consider the renormalization of the matrix elements of the bilinear quark operators (psi) over bar psi, (psi) over bar gamma(mu)psi, and (psi) over bar sigma(mu nu)psi at next-to-next-to-next-to-leading order in QCD perturbation theory at the symmetric subtraction point. This allows us to obtain conversion factors between the (MS) over bar scheme and the regularization invariant symmetric momentum subtraction (RI/SMOM) scheme. The obtained results can be used to reduce the errors in determinations of quark masses from lattice QCD simulations. The results are given in Landau gauge. (c) 2020 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).