We examine the influence of quadrupole moment of a slowly rotating neutron star (NS) on the oscillations of a fluid accretion disk (torus) orbiting a compact object the spacetime around which is described by the Hartle-Thorne geometry. Explicit formulae for non-geodesic orbital epicyclic and precession frequencies, as well as their simplified practical versions that allow for an expeditious application of the universal relations determining the NS properties, are obtained and examined. We demonstrate that the difference in the accretion disk precession frequencies for NSs of the same mass and angular momentum, but different oblateness, can reach up to tens of percent. Even higher differences can arise when NSs with the same mass and rotational frequency, but different equations of state (EoS), are considered. In particular, the Lense-Thirring precession frequency defined in the innermost parts of the accretion region can differ by more than one order of magnitude across NSs with different EoS. Our results have clear implications for models of the LMXBs variability.