The T2K collaboration has recently reported their results, which gives the best-fit values of the atmospheric mixing angle $sin^2 theta_{23} = 0.53$ and the Dirac CP-violating phase $delta = -1.89$ for normal neutrino mass ordering. We give a possible theoretical origin of such values based on the $mu$-$tau$ reflection symmetry. It has been found that the breaking of such symmetry using one-loop renormalization-group equations (RGEs) in the framework of minimal supersymmetric standard model can fit well with the latest T2K results and the recent global fits. To make a quantitative analysis, we have included for the first time the complete dataset of oscillation, beta decay, neutrinoless double-beta decay and cosmological observations in comparing the theory to experiments. We also further examine the importance of such breaking patterns for neutrinoless double-beta decay experiments.