We explore predictions of the wounded quark model for particle production and properties of the initial state formed in ultrarelativistic heavy-ion collisions. The approach is applied uniformly to A+A collisions in a wide collision energy range, as well as for p+A and p+p collisions at the CERN Large Hadron Collider (LHC). We find that generically the predictions from wounded quarks for such features as eccentricities or initial sizes are close (within 15\%) to predictions of the wounded nucleon model with an amended binary component. A larger difference is found for the size in p+Pb system, where the wounded quark model yields a smaller (more compact) initial fireball than the standard wounded nucleon model. The inclusion of subnucleonic degrees of freedom allows us to analyze p+p collisions in an analogous way, with predictions that can be used in further collective evolution. The approximate linear dependence of particle production in A+A collisions on the number of wounded quarks, as found in previous studies, makes the approach based on wounded quarks natural. Importantly, at the LHC energies we find approximate uniformity in particle production from wounded quarks, where at a given collision energy per nucleon pair similar production of initial entropy per source is needed to explain the particle production from p+p collisions up to A+A collisions. We also discuss the sensitivity of the wounded quark model predictions to distribution of quarks in nucleons, distribution of nucleons in nuclei, and to the quark-quark inelasticity profile in the impact parameter. In our procedure, the quark-quark inelasticity profile is chosen in such a way that the experiment-based parametrization of the proton-proton inelasticity profile is properly reproduced. The parameters of the overlaid multiplicity distribution is fixed from p+p and p+Pb data.