Given the important role that the galaxy bispectrum has recently acquired in cosmology and the scale and precision of forthcoming galaxy clustering observations, it is timely to derive the full expression of the large-scale bispectrum going beyond approximated treatments which neglect integrated terms or higher-order bias terms or use the Limber approximation. On cosmological scales, relativistic effects that arise from observing on the past light-cone alter the observed galaxy number counts, therefore leaving their imprints on N-point correlators at all orders. Working in spherical Bessel coordinates, in this talk I will show that it is possible to derive a compact expression for the power spectrum and bispectrum that encompasses all the physical effects at first and second order, including integrated (along the line of sight) terms.