Turbulent transport in a high ion temperature discharge of the Large Helical Device (LHD) is investigated by means of electromagnetic gyrokinetic simulations, which include kinetic electrons, magnetic perturbations, and full geometrical effects. Including kinetic electrons enables us to firstly evaluate the particle and the electron heat fluxes caused by turbulence in LHD plasmas. It is found that the electron energy transport reproduces the experimental result, and that the particle flux is negative. The contribution of magnetic perturbation to the transport is small because of very low beta. The turbulence is driven by the ion temperature gradient instability, and the effect of kinetic electrons enhances the growth rate larger than that from the adiabatic electron calculation. The ion energy flux is larger than that observed in the experiment, while the flux is close to the experimental observation when the temperature gradient is reduced 20% in the simulation. This significant sensitivity of the energy flux implies that the profile in the experiment is close to the critical temperature gradient. The critical gradient for turbulent energy flux is similar to that for the linear instability, i.e., the Dimits shift is small. This is because the zonal flow in the LHD is weaker than that in tokamaks.