Specific matrix elements of exchange and correlation kernels in time-dependent density-functional theory are computed. The knowledge of these matrix elements not only constrains approximate timedependent functionals. It also allows one to link different practical approaches to excited states, based either on density-functional theory or on many-body perturbation theory, despite the approximations that have been performed to derive them. [S0031-9007(99) [10,11,15], due to Görling and Levy, builds a perturbation theory (GLPT) in the difference between the many-body and the second-quantized Kohn-Sham (KS) Hamiltonians, where the parameter of the perturbation is the coupling constant of the particle interaction, in such a way that, at each order, the exact density is recovered.As an alternative to these DFT-based efforts, one may start from many-body perturbation theory and perform partial resummations of diagrams such as to build a screened interaction between dressed particles (quasiparticles) [16]. At the lowest order in the screened interaction, one obtains Hedin's GW approximation [17] to quasiparticle energies (one particle is added or subtracted to the system), while the energy of excited states for which the number of particles is conserved can be deduced, in a subsequent step, from a Bethe-Salpeter (BS) equation describing the interaction between quasiparticles (e.g., electron-hole pairs). The application of these techniques to real materials is very demanding [18].In the present Letter, we explore the relationships between excitation energies derived from time-dependent density-functional theory, the Görling-Levy perturbation theory, and the screened-interaction many-body perturbation theory. If excitation energies were derived without approximations, in the three formalisms, the results should be identical. Because of the specific theoretical developments, the most obvious simplifications are different, so that the practical schemes derived from these formalisms also differ. We find that some approximations used for practical calculations leave a connection between the approaches, at variance with the usual adiabatic local-density approximation (ALDA) in TDDFT [4][5][6], that leads to a different physical picture.The expressions from different approaches are linked thanks to a new technique for computing selected elements of TD functional kernels at resonance, that is, at the frequencies corresponding to differences in KS eigenvalues, for which the independent-particle susceptibility of the KS system is resonant. It is first applied to the exact exchange kernel, whose matrix elements appear in a simplified TDDFT treatment of excitation energies based on a Laurent expansion, and are found identical to the first-order corrections to KS eigenvalues differences, in the GLPT. The knowledge of these matrix elements imposes a new constraint on approximate functionals. The same technique is then applied to an explicit exchangecorrelation (XC) functional that includes an approximate correlation contribution. It i...