We study Jacobi matrices on trees whose coefficients are generated by multiple orthogonal polynomials. Hilbert space decomposition into an orthogonal sum of cyclic subspaces is obtained. For each subspace, we find generators and the generalized eigenfunctions written in terms of the orthogonal polynomials. The spectrum and its spectral type are studied for large classes of orthogonality measures.