Route to non-Gaussian statistics in convective turbulence