The formation of the first galaxies in the Universe is the new frontier of both galaxy formation and reionization studies. In fact, we will soon directly observe primeval galaxies thanks to the James Webb Space Telescope, and witness the reionization process through 21cm intensity mapping experiments. This unique moment in human history creates a fierce new challenge, i.e. to simultaneously understand in a unique and coherent picture the processes of galaxy formation and reionization, and – crucially – their connection. The latter, in particular, has escaped past numerical efforts. In this talk I will present the first results on this front from an years-long effort geared toward achieving such comprehensive picture, culminated in the Thesan suite of cosmological radiation-magneto-hydrodynamical simulations. I will briefly introduce the features of Thesan, highlighting the successes and failures of its physical model. Thesan produces realistic galaxy populations thanks to state-of-the-art physics, including self-consistent dust production+destruction and radiation transport. I will then show how Thesan can, for the first time, reproduce the connection between IGM and galaxies, as measured from the modulation of the Lyman-alpha flux around galaxies. Finally, I will chart the way forward towards and even deeper understanding of the emergence of the first structures in the Universe.