The notion of second-order relative spectrum of a self-adjoint operator acting on a Hilbert space has been studied recently in connection with the phenomenon of spectral pollution in the Galerkin method. In this paper we examine how the second-order spectrum encodes precise information about the multiplicity of the isolated eigenvalues of the underlying operator. Our theoretical findings are supported by various numerical experiments on the computation of guaranteed eigenvalue inclusions via finite element bases.