An analytical study of diffusiophoresis in a homogeneous suspension of identical spherical charge-regulating particles with an arbitrary thickness of the electric double layers in a solution of a symmetrically charged electrolyte with a uniform prescribed concentration gradient is presented. The charge regulation due to association/dissociation reactions of ionogenic functional groups on the particle surface is approximated by a linearized regulation model, which specifies a linear relationship between the surface charge density and the surface potential. The effects of particle-particle electrohydrodynamic interactions are taken into account by employing a unit cell model, and the overlap of the double layers of adjacent particles is allowed. The electrokinetic equations that govern the electric potential profile, the ionic concentration distributions, and the fluid flow field in the electrolyte solution surrounding the particle in a unit cell are linearized assuming that the system is only slightly distorted from equilibrium. Using a regular perturbation method, these linearized equations are solved with the equilibrium surface charge density (or zeta potential) of the particle as the small perturbation parameter. Closed-form formulas for the diffusiophoretic velocity of the charge-regulating sphere correct to the second order of its surface charge density or zeta potential are derived. Our results indicate that the charge regulation effect on the diffusiophoretic mobility is quite sensitive to the boundary condition for the electric potential specified at the outer surface of the unit cell. For the limiting cases of a very dilute suspension and a very thin or very thick electric double layer, the particle velocity is independent of the charge regulation parameter.