We estimate in medium properties of axion i.e., its mass and self-coupling within a three flavor Polyakov loop extended Nambu–Jona-Lasinio (PNJL) model with Kobayashi-Maskawa-t’Hooft determinant interaction. We also estimate the topological susceptibility of strong interaction within the same model. It is observed that (statistical) confinement effects simulated by Polyakov loop potential play an important role in the estimation of all these quantities, particularly, near the critical temperature. Both the mass and the self-coupling of the axion get correlated with the chiral and deconfinement transition. The results for all these quantities obtained within the PNJL model is compared with chiral perturbation theory, Nambu–Jona-Lasinio (NJL) model and lattice QCD simulation results wherever available. The results for properties of axions at finite baryon densities are also presented.