  The correlation function of meson interpolators in $N_f=2+1+1$ lattice QCD
with optimal domain-wall quarks at the physical point are studied for six
temperatures in the range $ T \sim$ 190-770 MeV. The meson interpolators
include a complete set of Dirac bilinears, and each for six combinations of
quark flavors. In this paper, we focus on the meson correlators of $u$ and $d$
quarks, and discuss their implications for the effective restoration of
$U(1)_A$ and $SU(2)_L \times SU(2)_R$ chiral symmetries, as well as the
emergence of the approximate $SU(2)_{CS}$ chiral spin symmetry.
