An extension of the two Higgs doublet model including inverse seesaw neutrinos and neutral Higgs bosons was constructed based on the $A_4$ symmetry in order to explain the recent neutrino oscillation data. This model can distinguish two well-known normal and inverted order schemes of neutrino data once both the effective masses $m_{beta}$ in tritium beta decays and $langle mrangle$ in the neutrinoless double beta decay are observed. The lepton flavor violating decays of the charged leptons $e_brightarrow e_agamma$, $murightarrow3e$, the Standard model-like Higgs boson decays $hrightarrow e_be_a$, and the $mu$-e