Direct numerical simulation (DNS) of a fully developed turbulent channel flow for various Reynolds numbers has been carried out to investigate the Reynolds number dependence. The Reynolds number is set to be Reτ=180, 395, and 640, where Reτ is the Reynolds number based on the friction velocity and the channel half width. The computation has been executed with the use of the finite difference method. Various turbulence statistics such as turbulence intensities, vorticity fluctuations, Reynolds stresses, their budget terms, two-point correlation coefficients, and energy spectra are obtained and discussed. The present results are compared with the ones of the DNSs for the turbulent boundary layer and the plane turbulent Poiseuille flow and the experiments for the channel flow. The closure models are also tested using the present results for the dissipation rate of the Reynolds normal stresses. In addition, the instantaneous flow field is visualized in order to examine the Reynolds number dependence for the quasi-coherent structures such as the vortices and streaks.