We study the discrete spectrum of fiber Hamiltonians \(H_{\mu _1\mu _2}(K)\) for a two-boson lattice model on \(\mathbb {Z}^3\) with on-site and nearest-neighbour interactions. The operators depend on the total quasi-momentum \(K\in \mathbb {T}^3\) and the real coupling constants \(\mu _1,\mu _2,\) and act on the Hilbert space \(L^{2,\textrm{e}}(\mathbb {T}^3)\) of even square-integrable functions. Each fiber Hamiltonian is a self-adjoint rank-four perturbation of the multiplication operator by the two-particle fiber dispersion relation, so its essential spectrum is explicitly determined. For \(K=0,\) the space \(L^{2,\textrm{e}}(\mathbb {T}^3)\) decomposes into three orthogonal invariant subspaces. In one subspace, the \((\mu _1,\mu _2)\) -plane is divided by smooth hyperbolic curves into regions where the number of eigenvalues below the lower edge of the essential spectrum is constant and equal to 0, 1 or 2. On these curves, eigenvalues reach the lower threshold of the essential spectrum and generate threshold resonances. In each of the remaining two subspaces, a unique critical point on the \(\mu _2\) -axis separates the parameter regimes with 0 and 1 eigenvalue below the essential spectrum; at this point, the lower threshold becomes an eigenvalue. An analogous description holds for the number of eigenvalues above the upper edge of the essential spectrum, together with the corresponding threshold phenomena at the upper spectral edge. In addition, we obtain a complete classification of the discrete spectrum of \(H_{\mu _1\mu _2}(0)\) describing the exact number of eigenvalues lying below and above the essential spectrum in each region of the parameter plane. Moreover, we establish stability of these eigenvalues under perturbations of the total quasi-momentum \(K\in \mathbb {T}^3\) : for \(K\ne 0\) , the operator \(H_{\mu _1\mu _2}(K)\) has at least as many eigenvalues below and above the essential spectrum as \(H_{\mu _1\mu _2}(0).\) If the lower edge of the essential spectrum is a threshold resonance of \(H_{\mu _1\mu _2}(0)\) , then \(H_{\mu _1\mu _2}(K)\) has at least one more eigenvalue below the essential spectrum than \(H_{\mu _1\mu _2}(0).\) If the lower edge is a threshold eigenvalue of multiplicity two then \(H_{\mu _1\mu _2}(K)\) has at least two more such eigenvalues. Analogous statement holds for the upper edge of the essential spectrum.