We proved direct and inverse theorems on B-spline quasi-interpolation sampling representation with a Littlewood-Paley-type norm equivalence in Sobolev spaces $W_p^r$ of mixed smoothness $r$. Based on this representation, we established estimates of the approximation error of recovery in $L_q$-norm of functions from the unit ball $U_p^r$ in the spaces $W_p^r$ by linear sampling algorithms and the asymptotic optimality of these sampling algorithms in terms of Smolyak sampling width $r_n^s(U_p^r,L_q)$ and sampling width $r_n(U_p^r,L_q)$.