The most promising mechanism of generating PBHs is by the enhancement of power spectrum of the primordial curvature perturbation, which is usually accompanied by the the enhancement of non-Gaussianity that crucially changes the abundance of PBHs. In this chapter I will discuss how non-Gaussianity is generated in single field inflation as well as in the curvaton scenario, and then discuss how to calculate PBH mass function and induced gravitational waves (GWs) with such non-Gaussianities. When the PBH abundance is fixed, non-Gaussianity only has mild impact on the spectral shape of the induced GWs, which gives relatively robust predictions in mHz and nHz GW experiments.