If the primordial curvature perturbation followed a Gaussian distribution, primordial black holes (PBHs) will be Poisson distributed with no additional clustering. We consider local non-Gaussianity and its impact on the initial PBH clustering and mass function due to mode coupling between long and short wavelength modes. We show that even a small amount of non-Gaussianity results in a significant enhancement on the PBH initial clustering and subsequent merger rate and that the PBH mass function shifts to higher mass PBHs.However, as the clustering becomes strong, the local number density of PBHs becomes large, leading to a large theoretical uncertainty in the merger rate.