Not yet reviewed by Pith; the record is open.

We obtain a model-independent prediction for the two-photon exchange contribution to the hyperfine splitting in muonic hydrogen. We use the relation of the Wilson coefficients of the spin-dependent dimension-six four-fermion operator of NRQED applied to the electron-proton and to the muon-proton sectors. Their difference can be reliably computed using chiral perturbation theory, whereas the Wilson coefficient of the electron-proton sector can be determined from the hyperfine splitting in hydrogen. This allows us to give a precise model-independent determination of the Wilson coefficient for the muon-proton sector, and consequently of the two-photon exchange contribution to the hyperfine splitting in muonic hydrogen, which reads $\delta \bar E_{p\mu,\rm HF}^{\rm TPE}(nS)=-\frac{1}{n^3}1.161(20)$ meV. Together with the associated QED analysis, we obtain a prediction for the hyperfine splitting in muonic hydrogen that reads $E^{\rm th}_{p\mu,\rm HF}(1S)=182.623(27)$ meV and $E^{\rm th}_{p\mu,\rm HF}(2S)=22.8123(33)$ meV. The error is dominated by the two-photon exchange contribution.