Mössbauer isomer shifts of 129 I and 127 I in the ICl, IBr and I 2 molecules are studied. Filatov's formulation is used, based on calculating the electronic energy change of the two systems involved in the Mössbauer γ transition, the source and absorber. The energy difference between the transitions in the two systems determines the shift. The effects of relativity and electron correlation on the shifts are investigated. The exact two-component (X2C) and the four-component relativistic schemes give virtually identical results; the non-relativistic approach yields about 50% of the relativistic shifts. Electron correlation is included by coupled-cluster singles-and-doubles with perturbative triples [CCSD(T)]; it reduces Hartree-Fock shifts by 15%-20%. Basis sets are increased until the isomer shifts converge. The final results, calculated with the converged basis in the framework of the X2C Hamiltonian and CCSD(T) correlation, give an agreement of 10% or better with experimental data.