Optimal target spaces are exhibited in arbitrary-order Sobolev type embeddings for traces of n n -dimensional functions on lower dimensional subspaces. Sobolev spaces built upon any rearrangement-invariant norm are allowed. A key step in our approach consists of showing that any trace embedding can be reduced to a one-dimensional inequality for a Hardy type operator depending only on n n and on the dimension of the relevant subspace. This can be regarded as an analogue for trace embeddings of a well-known symmetrization principle for first-order Sobolev embeddings for compactly supported functions. The stability of the optimal target space under iterations of Sobolev trace embeddings is also established and is part of the proof of our reduction principle. As a consequence, we derive new trace embeddings, with improved (optimal) target spaces, for classical Sobolev, Lorentz-Sobolev and Orlicz-Sobolev spaces.