We employ our Bayesian Machine Learning framework BINGO (Bayesian INference for Galactic archaeOlogy) to obtain high-quality stellar age estimates for 68 360 red giant and red clump stars present in the 17th data release of the Sloan Digital Sky Survey, the APOGEE-2 high-resolution spectroscopic survey. By examining the denoised age-metallicity relationship of the Galactic disc stars, we identify a drop in metallicity with an increase in [Mg/Fe] at an early epoch, followed by a chemical enrichment episode with increasing [Fe/H] and decreasing [Mg/Fe]. This result is congruent with the chemical evolution induced by an early-epoch gas-rich merger identified in the Milky Way-like zoom-in cosmological simulation Auriga. In the initial phase of the merger of Auriga 18 there is a drop in metallicity due to the merger diluting the metal content and an increase in the [Mg/Fe] of the primary galaxy. Our findings suggest that the last massive merger of our Galaxy, the Gaia-Sausage-Enceladus, was likely a significant gas-rich merger and induced a starburst, contributing to the chemical enrichment and building of the metal-rich part of the thick disc at an early epoch.