We study the Cauchy problem of a two-species chemotactic model. Using the Fourier frequency localization and the Bony paraproduct decomposition, we establish a unique local solution and blow-up criterion of the solution, when the initial data (u 0, v 0, w 0) belongs to homogeneous Besov spaces \(\dot B_{p,1}^{ - 2 + 3/p}\left( {{\mathbb{R}^3}} \right) \times \dot B_{r,1}^{ - 2 + 3/r}\left( {{\mathbb{R}^3}} \right) \times \dot B_{q,1}^{3/q}\left( {{\mathbb{R}^3}} \right)\) for p, q and r satisfying some technical assumptions. Furthermore, we prove that if the initial data is sufficiently small, then the solution is global. Meanwhile, based on the so-called Gevrey estimates, we particularly prove that the solution is analytic in the spatial variable. In addition, we analyze the long time behavior of the solution and obtain some decay estimates for higher derivatives in Besov and Lebesgue spaces.