Reconstructed Discontinuous Galerkin (rDG) methods are presented for solving diffusion equations based on a first-order hyperbolic system (FOHS) formulation. The idea is to combine the advantages of the FOHS formulation and the rDG methods in an effort to develop a more reliable, accurate, efficient, and robust method for solving the diffusion equations. The developed hyperbolic rDG methods can be made to have higher-order accuracy than conventional DG methods with fewer degrees of freedom. A number of test cases for different diffusion equations are presented to assess accuracy and performance of the newly developed hyperbolic rDG methods in comparison with the standard BR2 DG method. Numerical experiments demonstrate that the hyperbolic rDG methods are able to achieve the designed optimal order of accuracy for both solutions and their derivatives on regular, irregular, and heterogeneous girds, and outperform the BR2 method in terms of the magnitude of the error, the order of accuracy, the size of time steps, and the CPU times required to achieve steady state solutions, indicating that the developed hyperbolic rDG methods provide an attractive and probably an even superior alternative for solving the diffusion equations.