The innermost stable cicular orbit (ISCO) of an accretion disc that orbits a neutron star (NS) is often assumed to be a unique prediction of general relativity. However, it has been argued that ISCO also appears around highly elliptic bodies described by Newtonian theory. In this sense, the behaviour of an ISCO around a rotating oblate neutron star is formed by the interplay between relativistic and Newtonian effects. Here we briefly explore the consequences of this interplay using a straightforward analytic approach as well as numerical models that involve modern NS equations of state. We examine the ratio K between the ISCO radius and the radius of the neutron star. We find that, with growing NS spin, the ratio K first decreases, but then starts to increase. This non-monotonic behaviour of K can give rise to a neutron star spin interval in which ISCO appears for two very different ranges of NS mass. This may strongly affect the distribution of neutron stars that have an ISCO (ISCO-NS). When (all) neutron stars are distributed around a high mass M 0 , the ISCO-NS spin distribution is roughly the same as the spin distribution corresponding to all neutron stars. In contrast, if M 0 is low, the ISCO-NS distribution can only have a peak around a high value of spin. Finally, an intermediate value of M 0 can imply an ISCO-NS distribution divided into two distinct groups of slow and fast rotators. Our findings have immediate astrophysical applications. They can be used for example to distinguish between different models of high-frequency quasiperiodic oscillations observed in low-mass NS X-ray binaries.