Since this is a very active research topic the answer here might change regularly. Just within the last few months we've seen updatedupdates to the SCAN functional in the form of r2SCAN and the de-orbitalized r2SCAN-L. These are both functionals that attempt to satisfy all known physical constraints that can be in-principle satisfied by a semi-local functional, although it should be noted that r2SCAN relaxes one of these constraints (the fourth-order gradient expansion) for the benefit of numerical accuracy. To achieve a "universal" functional, these exact constraints have to be satisfied, and SCAN is the closest we have at least in terms of a semi-local functional.
However, there are limits to how good a semi-local functional can be even in principle, and ultimately non-local information is required for a true "universal functional." This is what has made the hybrid functionals so popular. In the solid-state community, HSE06 has become a de facto standard, but its universality is limited by including a fixed amount of Hartree-Fock exchange, when the true value of mixing is expected to vary according to the material's dielectric properties. There have been some "dielectric-dependent" functionals developed, as well as schemes to use HSE but to optimize the amount of mixing, and this is also an active research topic.