We consider the effect of scalarization on static and slowly rotating neutron stars for a wide variety of realistic equations of state, including pure nuclear matter, nuclear matter with hyperons, hybrid nuclear and quark matter, and pure quark matter. We analyze the onset of scalarization, presenting a universal relation for the critical coupling parameter versus compactness. We find that the onset and the magnitude of the scalarization are strongly correlated with the value of the gravitational potential (the metric component gtt) at the center of the star. We also consider the moment-of-inertia-compactness relations and confirm universality for the nuclear matter, hyperon and hybrid equations of state.