Since 
$$Z=X^\top (XX^\top + \mathrm{Id})^{-1} X=(Y + \mathrm{Id})^{-1} Y,$$ with $Y=X^\top X$, we can find the eigenvalue density of $Z$ from the eigenvalue density of $Y$, which for large matrix dimensions is given by the Marchenko-Pastur density $\rho_{\rm MP}(y)$.

In particular, the average eigenvalue is given by
$$\langle z\rangle=\int_0^\infty \frac{y}{1+y}\,\rho_{\rm MP}(y)\,dy=\frac{\sqrt{\gamma^2+4}+\gamma+2}{2 \gamma}.$$