Neutron stars are ideal to probe, not only nuclear physics, but also strong-field gravity. Approximate universal relations insensitive to the star's internal structure exist among certain observables and are useful in testing General Relativity, as they project out the uncertainties in the equation of state. One such set of universal relations between the moment of inertia $(I)$, the tidal Love number and the quadrupole moment $(Q)$ has been studied both in General Relativity and in modified theories. In this paper, we study the relations in dynamical Chern-Simons gravity, a well-motivated, parity-violating effective field theory, extending previous work in various ways. First, we study how projected constraints on the theory using the I-Love relation depend on the measurement accuracy of $I$ with radio observations and that of the Love number with gravitational-wave observations. Provided these quantities can be measured with future observations, we find that the latter could place bounds on dynamical Chern-Simons gravity that are six orders of magnitude stronger than current bounds. Second, we study the I-Q and Q-Love relations in this theory by constructing slowly-rotating neutron star solutions to quadratic order in spin. We find that the approximate universality continues to hold in dynamical Chern-Simons gravity, and in fact, it becomes stronger than in General Relativity, although its existence depends on the normalization of the dimensional coupling constant of the theory. Finally, we study the variation of the eccentricity of isodensity contours inside a star and its relation to the degree of universality. We find that, in most cases, the eccentricity variation is smaller in dynamical Chern-Simons gravity than in General Relativity, providing further support to the idea that the approximate self-similarity of isodensity contours is responsible for universality.