We propose a novel numerical approach for the Boltzmann equation with uncertainties. The method combines the efficiency of classical direct simulation Monte Carlo (DSMC) schemes in the phase space together with the accuracy of stochastic Galerkin (sG) methods in the random space. This hybrid formulation makes it possible to construct methods that preserve the main physical properties of the solution along with spectral accuracy in the random space. The schemes are developed and analyzed in the case of space homogeneous problems as these contain the main numerical difficulties. Several test cases are reported, both in the Maxwell and in the variable hard sphere (VHS) framework, and confirm the properties and performance of the new methods.