Even before the introduction of Conformal Field Theory by Belavin, Polyakov and Zamolodchikov, it appeared indirectly in the work of den Nijs and Nienhuis using Coulomb gas techniques. The latter postulate (unrigorously) that height functions of lattice models of statistical mechanics (like percolation, Ising, 6-vertex models etc) converge to the Gaussian Free Field, allowing to derive many exponents and dimensions.This convergence remains in many ways mysterious, in particular it was never formulated in the presence of a boundary, but rather on a torus or a cylinder. We will discuss the original arguments as well as some recent progress, including possible formulations on general domains or Riemann surfaces and their relations to CFT, SLE and conformal invariance of critical lattice models. Interestingly, new objects in complex geometry and potential theory seem to arise.