We discuss a family of multispecies ferromagnetic Ising models on multiregular random graphs. In the large volume limit, thermodynamic quantities are related to the solution of a belief propagation (BP) fixed point equation. A phase transition is identified and the critical region is determined by the spectral radius of a finite-dimensional matrix.