The significance of new physics appearing in the loop-induced decays of neutral Higgs bosons into pairs of di-bosons γγ and Zγ will be discussed in the framework of the 3-3-1 models based on a recent work [1], where the Higgs sector becomes effectively the same as that in the Two Higgs Doublet models (2HDM) after the first symmetry breaking from SU (3) L scale into the electroweak scale. For large SU (3) L scale v 3 10 TeV, dominant one-loop contributions to the two decay amplitudes arise from only the single charged Higgs boson predicted by the 2HDM, leading to that experimental constraint on the signal strength µ 331 γγ of the Standard model like Higgs boson decay h → γγ will result in a strict upper bound on the signal strength µ 331 Zγ of the decay h → Zγ. For a particular model with lower v 3 around 3 TeV, contributions from heavy charged gauge and Higgs bosons may have the same order, therefore may give strong destructive or constructive correlations.As a by product, deviations from the SM prediction |µ 331 γγ − 1| ≤ 0.04 still allows |µ 331 Zγ − 1| to reach values near 0.1. We also show that there exists an CP-even neutral Higgs boson h 0 3 predicted by the 3-3-1 models, but beyond the 2HDM, has an interesting property that the branching ratio Br(h 0 3 → γγ) is very sensitive to the parameter β used to distinguish different 3-3-1 models. ¶ Corresponding author * Electronic address: hathanhhung@hpu2.edu.vn 1 arXiv:1907.06735v1 [hep-ph]