Divergence-form operators with stationary random coefficients homogenize over large scales. We investigate the effect of certain perturbations of the medium on the homogenized coefficients. The perturbations that we consider are rare at the local level, but when occurring, have an effect of the same order of magnitude as the initial medium itself. The main result of the paper is a first-order expansion of the homogenized coefficients, as a function of the perturbation parameter.