We investigate time-decay rates of strong solutions to the 3D compressible magnetohydrodynamics equations with large initial data. The main novelty of this paper is two-fold: first, we prove the upper optimal decay rates of the higher-order spatial derivatives of the solution, which are the same as those of the heat equation, and faster than the decay rates in the previous related works. Second, if the initial data satisfy some additional low frequency assumption, we also show the lower optimal decay rates of the solution as well as its all-order spatial derivatives. Therefore, our decay rates are optimal in this sense. Our methods mainly involve the Fourier splitting method, low-frequency and high-frequency decomposition and delicate energy estimates.