In this paper, the high-order space-time discontinuous Galerkin cell vertex scheme (DG-CVS) developed by the authors for hyperbolic conservation laws is extended for time dependent diffusion equations. In the extension, the treatment of the diffusive flux is exactly the same as that for the advective flux. Thanks to the Riemann-solver-free and reconstruction-free features of DG-CVS, both the advective flux and the diffusive flux are evaluated using continuous information across the cell interface. As a result, the resulting formulation with diffusive fluxes present is still consistent and does not need any extra ad hoc techniques to cure the common "variational crime" problem when traditional DG methods are applied to diffusion problems. For this reason, DG-CVS is conceptually simpler than other existing DG-typed methods. The numerical tests demonstrate that the convergence order based on the L²-norm is optimal, i.e. O(h^p+1) for the solution and O(h^p) for the solution gradients, when the basis polynomials are of odd degrees. For even-degree polynomials, the convergence order is sub-optimal for the solution and optimal for the solution gradients. The same odd-even behaviour can also be seen in some other DG-typed methods.