The leading-twist parton distribution amplitudes (PDAs) of ground-state and 131b4f9fde60977011b8cb83c041">- and -quarkonia are calculated using a symmetry-preserving continuum treatment of the meson bound-state problem which unifies the properties of these heavy-quark systems with those of light-quark bound-states, including QCD's Goldstone modes. Analysing the evolution of and 131b4f9fde60977011b8cb83c041"> PDAs with current-quark mass, , increasing away from the chiral limit, it is found that in all cases there is a value of for which the PDA matches the asymptotic form appropriate to QCD's conformal limit and hence is insensitive to changes in renormalisation scale, ζ. This mass lies just above that associated with the s-quark. At current-quark masses associated with heavy-quarkonia, on the other hand, the PDAs are piecewise convex–concave–convex. They are much narrower than the asymptotic distribution on a large domain of ζ ; but nonetheless deviate noticeably from , which is the result in the static-quark limit. There are also material differences between and 131b4f9fde60977011b8cb83c041"> PDAs, and between the PDAs for different vector-meson polarisations, which vanish slowly with increasing ζ . An analysis of moments of the root-mean-square relative-velocity, 〈v2m〉, in and 131b4f9fde60977011b8cb83c041"> systems reveals that 〈v4〉-contributions may be needed in order to obtain a reliable estimate of matrix elements using such an expansion, especially for processes involving heavy pseudoscalar quarkonia.