  We report on our latest results in the calculation of the three-loop heavy
flavor contributions to the Wilson coefficients in deep-inelastic scattering in
the asymptotic region $Q^2 \gg m^2$. We discuss the different methods used to
compute the required operator matrix elements and the corresponding Feynman
integrals. These methods very recently allowed us to obtain a series of new
operator matrix elements and Wilson coefficients like the flavor non-singlet
and pure singlet Wilson coefficients.
