Certain subclasses \(\mathcal{S}\mathcal{H}^{m,\ell }_{p,q}[\Phi _{i}, \Psi _{j}; \alpha , A, B]\) and \(\overline{\mathcal{S}\mathcal{H}}^{m,\ell }_{p,q}[\Phi _{i}, \Psi _{j}; \alpha , A, B]\) of the function class \(\mathcal{S}\mathcal{H}(m)\) of multivalent harmonic functions associated with the (p, q)-derivative operator and subordination function are introduced and investigated. Further, on the one hand, we obtain the necessary and sufficient conditions and the coefficient characterization for the class of (p, q)-Sălăgean multivalent harmonic functions for the subclasses. On the other hand, we also provide the bound estimates of the coefficient \(a_{n}\) and \(b_{n}\) and distortion theorems and covering theorems as well as the closeness of convex combination for these harmonic function classes.