  We present a systematic study of galaxy bias in the presence of primordial
non-Gaussianity in General Relativity (GR) at second order in perturbation
theory. The non-linearity of the Poisson equation in GR and primordial
non-Gaussianity are consistently included. We show that the inclusion of
non-local primordial non-Gaussianity in addition to local non-Gaussianity is
important to show the absence of the modulation of small scale clustering by
the long-wavelength mode in the single field slow-roll inflation. We study the
bispectrum of the relativistic galaxy density in several gauges and identify
the effect of primordial non-Gaussianity and GR corrections.
