Abstract: The wide measured range of thermal conductivities (k) for monolayer MoS2 and the corresponding incongruent calculated values in the literature all suggest that extrinsic defect thermal resistance is significant and varied in synthesized samples of this material. Here we present defect-mediated thermal transport calculations of MoS2 using interatomic forces derived from density functional theory combined with Green's function methods to describe phonon-point-defect interactions and a Peierls-Boltzmann formalism for transport. Conductivity calculations for bulk and monolayer MoS2 using different density functional formalisms are compared. Nonperturbative first-principles methods are used to describe defect-mediated spectral functions, scattering rates, and phonon k, particularly from sulfur vacancies (VS), and in the context of the plethora of measured and calculated literature values. We find that k of monolayer MoS2 is sensitive to phonon-VS scattering in the range of experimentally observed densities, and that first-principles k calculations using these densities can explain the range of measured values found in the literature. Furthermore, measured k values for bulk MoS2 are more consistent because VS defects are not as prevalent.