We study the interaction of the pseudoscalar glueball with scalar and
pseudoscalar quark-antiquark meson fields as well as their first excited states. We introduce a corresponding chiral Lagrangian that describes the two- and three-body decays of a pseudoscalar glueball, $J^{PC}=0^{-+}$, with a mass of $2.6$ GeV as predicted by lattice QCD simulations. We compute the decay of the pseudoscalar glueball into (pseudo)scalar and their excited states. The various branching ratios are parameters-free predictions. These states and channels are in reach of the ongoing BESIII, Belle II, LHCb, and NICA experiments and the upcoming PANDA experiment at the FAIR facility.