The breakup of viscous compound threads in the presence of insoluble surfactant at both interfaces is investigated. We use asymptotic methods in the limit of long axisymmetric waves to derive a coupled system of five one-dimensional ͑1-D͒ partial differential equations governing the evolution of the outer and inner interfaces, the surfactant concentrations there, and the leading order axial velocity component in the jet. The linear, and nonlinear, stability of these equations is then investigated for a wide range of outer to inner viscosity ratio, m, outer to inner surface tension ratio, ␥, the ratio of initial outer to inner radii, ␣, initial surfactant concentrations at the outer and inner interfaces, ⌫ 1 0 and ⌫ 2 0 , surfactant activities,  1 and  2 , and the Schmidt numbers, Sc 1 and Sc 2 , defined as the ratio of the kinematic viscosity to the surfactant surface diffusion coefficient. We also show that if Sc 1 ϭSc 2 , these results are recovered via solution of 1-D evolution equations governing the dynamics of an effective single surfactant covered thread, which are obtained through appropriate rescalings; these rescalings are detailed herein.