The dipole moment of a charge distribution $\rho(\mathbf r)$, which is always defined as
$$
\mathbf d=\int\mathbf r\rho(\mathbf r) \: \mathrm d\mathbf r
$$
with respect to some coordinate origin, is independent of said coordinate origin if and only if the system is neutral. Otherwise, as you note, it transforms to a new origin $\mathbf r_0$ as
$$
\mathbf d'
=\int(\mathbf r-\mathbf r_0)\rho(\mathbf r) \: \mathrm d\mathbf r
=\mathbf d-\mathbf r_0\int\rho(\mathbf r) \: \mathrm d\mathbf r
=\mathbf d-Q\mathbf r_0
$$
for $Q$ the total charge of the system.
(This transformation law comes from the fact that in our original equation, $\mathbf r$ actually plays two roles, that of a position and that of an arrow. In this calculation of $\mathbf d'$ I'm keeping this calculation in the old reference frame, which keeps the $\rho(\mathbf r)$ intact, but I need to change the arrow to $\mathbf r-\mathbf r_0$ so that it now points from the new coordinate origin, $\mathbf r_0$, to the position $\mathbf r$.)
For a charged system, the dipole moment takes a backseat role, but it is still a useful quantity - for example, it is still the coefficient in the $1/r^2$ term of the multipole expansion for the electrostatic potential far away from the system,
$$
V(\mathbf r) = \frac{Q}{r} + \frac{\mathbf d\cdot\mathbf r}{r^3}+\cdots,
$$
with quadrupole terms after that. However, it is always important to keep in mind this dependence on the coordinate origin, and to clearly state this origin whenever the dipole moment is reported for a charged system.