We use self-consistent mean-field methods in combination with the interacting boson model (IBM) of nuclei, to establish a linkage between universal energy density functionals (EDFs) and partial dynamical symmetry (PDS). An application $^{168}$Er shows that IBM Hamiltonians derived microscopically from known nonrelativistic relativistic EDFs this region, conform SU(3)-PDS.