In this work, we study the approximation properties of multipatch dG-IgA methods, that apply the multipatch IsogeometricAnalysis discretization concept and the discontinuous Galerkin technique on the interfaces between the patches, for solving linear diffusion problems with diffusion coefficients that may be discontinuous across the patch interfaces. The computational domain is divided into nonoverlapping subdomains, called patches in IgA, where Bsplines, or NURBS approximations spaces are constructed. The solution of the problem is approximated in every subdomain without imposing any matching grid conditions and without any continuity requirements for the discrete solution across the interfaces. Numerical fluxes with interior penalty jump terms are applied in order to treat the discontinuities of the discrete solution on the interfaces. We provide a rigorous a priori discretization error analysis for diffusion problems in two- and three-dimensional domains, where solutions patchwise belong to Wl,p, with some l ≥ 2 and p ∈ (2d/(d + 2(l - 1)), 2]. In any case, we show optimal convergence rates of the discretization with respect to the dG - norm. [ABSTRACT FROM AUTHOR]