  Radiative quarkonia decays offer an ideal setting for probing Axion-Like
Particle (ALP) interactions. This paper provides a comprehensive review of ALP
production mechanisms through the $e^+ e^- \to \gamma\,a$ process at B- and
Charm-factories, alongside an analysis of potential ALP decay channels. We
derive constraints on ALP couplings to Standard Model (SM) fields, based on
recent experimental results on quarkonia decays by the Belle II and BESIII
collaborations. The analysis distinguishes between "invisible" and "visible"
ALP decay scenarios. The "invisible" scenario, characterised by a mono-$\gamma$
plus missing-energy signature, enables stringent limits on ALP-photon and
ALP-quark ($b$ or $c$) couplings. Moreover, extensive research at flavour
factories has explored various "visible" ALP decays into SM final states, which
depend on a larger set of ALP-SM couplings. To streamline the "visible" ALP
scenario, we introduce additional theoretical assumptions, such as universal
ALP-fermion couplings, or we adopt specific benchmark ALP models, aiming to
minimise the number of independent variables in our analysis.
