The large scale structure bispectrum in the squeezed limit couples large with small scales. Since relativity is important at large scales and non-linear loop corrections are important at small scales, the proper calculation of the observed bispectrum in this limit requires a non-linear relativistic calculation. We compute the matter bispectrum in general relativity in the weak field approximation. The calculation is as involved as existing second-order results. We find several differences with the Newtonian calculation such as the non-cancellation of IR divergences, the need to renormalize the background, and the fact that initial conditions must be set at second order in perturbation theory. For the bispectrum, we find relativistic corrections to be as large as the newtonian result in the squeezed limit. In that limit relativistic one-loop contributions, which we compute for the first time, can be as large as tree level results and have the same 1/k2 dependence as a primordial local non-Gaussianity signal where k is the momentum approaching zero. Moreover, we find the time dependence of the relativistic corrections to the bispectrum to be the same as that of a primordial non-Gaussianity signal.