We introduce novel relations between derivatives of the Dirac eigenvalue spectrum with respect to the light sea quark mass and the $(n+1)$-point correlations among the eigenvalues of the massless Dirac operator. Using these relations we present LQCD results for the derivatives at light pion masses and at a temperature of about 1.6 times the chiral phase transition temperature. We find that eigenvalue density develops a peaked structure. We demonstrate that this phenomena is responsible for the manifestations of axial anomaly in two-point correlation functions of light scalar and pseudoscalar mesons. After continuum and chiral extrapolations we find that axial anomaly remains manifested in two-point correlation functions of scalar and pseudoscalar mesons in the chiral limit.