The aim of this work is to provide a mathematical model to describe the early stages of the embryonic development of zebrafish posterior lateral line (PLL). In particular, we focus on evolution of PLL proto-organ (said primordium), from its formation to the beginning of the cyclical behavior that amounts in the assembly of immature proto-neuromasts towards its caudal edge accompanied by the deposition of mature proto-neuromasts at its rostral region. Our approach has an hybrid integro-differential nature, since it integrates a microscopic/discrete particle-based description for cell dynamics and a continuous description for the evolution of the spatial distribution of chemical substances (i.e., the stromal-derived factor SDF1a and the fibroblast growth factor FGF10). Boolean variables instead implement the expression of molecular receptors (i.e., Cxcr4/Cxcr7 and fgfr1). Cell phenotypic transitions and proliferation are included as well. The resulting numerical simulations show that the model is able to qualitatively and quantitatively capture the evolution of the wild-type (i.e., normal) embryos as well as the effect of known experimental manipulations. In particular, it is shown that cell proliferation, intercellular adhesion, FGF10-driven dynamics, and a polarized expression of SDF1a receptors are all fundamental for the correct development of the zebrafish posterior lateral line.