The thermal restoration of chiral symmetry in QCD is known to proceed by an analytic crossover, which is widely expected to turn into a phase transition with a critical endpoint as the baryon density is increased. In the absence of a genuine solution to the sign problem of lattice QCD, simulations at zero and imaginary baryon chemical potential in a parameter space enlarged by a variable number of quark flavours and quark masses constitute a viable way to constrain the location of a possible non-analytic phase transition and its critical endpoint. In this article I review recent progress towards an understanding of the nature of the transition in the massless limit, and its critical temperature at zero density. Combined with increasingly detailed studies of the physical crossover region, current data bound a possible critical point to μB ≳ 3T.