We compute the skewness of the matter distribution arising from nonlinear evolution and from non-Gaussian initial perturbations. We apply our result to a very generic class of models with non-Gaussian initial conditions and we estimate analytically the ratio between the skewness due to nonlinear clustering and the part due to the intrinsic non-Gaussianity of the models. We finally extend our estimates to higher moments.