Universal relations for quasinormal modes of neutron stars in $R^2$ gravity