Determination of heat loads is a key issue in the design of gas turbines. In order to optimize the cooling, an exact knowledge of the heat flux and temperature distributions on the airfoils surface is necessary. Heat transfer is influenced by various factors, like pressure distribution, wakes, surface curvature, secondary flow effects, surface roughness, free stream turbulence, and separation. Each of these phenomenons is a challenge for numerical simulations. Among numerical methods, large eddy simulations (LES) offers new design paths to diminish development costs of turbines through important reductions of the number of experimental tests. In this study, LES is coupled with a thermal solver in order to investigate the flow field and heat transfer around a highly loaded low pressure water-cooled turbine vane at moderate Reynolds number (150,000). The meshing strategy (hybrid grid with layers of prisms at the wall and tetrahedra elsewhere) combined with a high fidelity LES solver gives accurate predictions of the wall heat transfer coefficient for isothermal computations. Mesh convergence underlines the known result that wall-resolved LES requires discretizations for which y+ is of the order of one. The analysis of the flow field gives a comprehensive view of the main flow features responsible for heat transfer, mainly the separation bubble on the suction side that triggers transition to a turbulent boundary layer and the massive separation region on the pressure side. Conjugate heat transfer computation gives access to the temperature distribution in the blade, which is in good agreement with experimental measurements. Finally, given the uncertainty on the coolant water temperature provided by experimentalists, uncertainty quantification allows apprehension of the effect of this parameter on the temperature distribution.