Symbolic Regression

7+ years ago I posted a question here regarding symbolic regression in Mathematica. At that time there did not seem to be much built-in support for this. I’m interested in multi-input models of the form y=f(x1,x2,x3,x4,x5) where I need to find the best fit "f" given a data set.

Has any meaningful progress been made during the last 7 years in this regard in Mathematica ? I’m aware of TuringBot. Is that considered the state of the art in symbolic regression ?

Thanks for any knowledge on this subject.