We investigate the behavior of sulfur vacancy defects, the most abundant type of intrinsic defect in monolayer MoS2, using first-principles calculations based on density functional theory. We consider the dependence of the isolated defect formation energy on the charge state and on uniaxial tensile and compressive strain up to 5%. We also consider the possibility of defect clustering by examining the formation energies of pairs of vacancies at various relative positions, and their dependence on charge state and strain. We find that there is no driving force for vacancy clustering, independent of strain in the material. The barrier for diffusion of S vacancies is larger than 1.9 eV in both charged and neutral states regardless of the presence of other nearby vacancies. We conclude that the formation of extended defects from S vacancies in planar monolayer MoS2 is hindered both thermodynamically and kinetically.