We propose a viable model based on the SU (3)C × SU (3)L × U (1)X gauge group supplemented by the S4 family symmetry and other auxiliary cyclic symmetries, whose spontaneous breaking gives rise to the observed pattern of SM fermion masses and mixing angles. In the proposed model the small light active neutrino masses are generated from a linear seesaw mechanism mediated by three Majorana neutrinos. The model is capable of reproducing the experimental values of the physical observables of both quark and lepton sectors. Our model is predictive in the quark sector having 9 effective parameters that allow to successfully reproduce the four CKM parameters and the six Standard Model (SM) quark masses. In the SM quark sector, there is particular scenario, motivated by naturalness arguments, which allows a a good fit for its ten observables, with only six effective parameters. We also study the single heavy scalar production via gluon fusion mechanism at proton-proton collider. Our model is also consistent with the experimental constraints arising from the Higgs diphoton decay rate.