In this paper the nonlinear multi-species Boltzmann equation with random uncertainty coming from the initial data and collision kernel is studied. Well-posedness and long-time behavior – exponential decay to the global equilibrium – of the analytical solution, and spectral gap estimate for the corresponding linearized gPC-based stochastic Galerkin system are obtained, by using and extending the analytical tools provided in [M. Briant and E. S. Daus, Arch. Ration. Mech. Anal., 3, 1367–1443, 2016] for the deterministic problem in the perturbative regime, and in [E. S. Daus, S. Jin and L. Liu, Kinet. Relat. Models, 12, 909–922, 2019] for the single-species problem with uncertainty. The well-posedness result of the sensitivity system presented here has not been obtained so far even for the single-species case.