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(2016) 1367–1443] for the deterministic problem in the perturbative regime, and in [E.S. Daus, S. Jin and L. Liu,Kinet. Relat. Models12(2019) 909–922] for the single-species problem with uncertainty. The well-posedness result of the sensitivity system presented here has not been obtained so far neither in the single species case nor in the multi-species case.