To understand the irradiation resistance and hydrogen (H) behavior in tungsten borides (WxBy) in a burning plasma fusion environment, the energetics of intrinsic point defects and H in six stable ground state W borides, including W2B, WB, WB2, W2B5, WB3 and WB4, have been investigated using first-principles density functional theory calculations. The results show that the formation energies of interstitials and vacancies do not directly depend on the W and B content in W borides. However, the interaction between vacancies of a stoichiometric vacancy (SV) cluster in WxBy is related to the atomic ratio of B to W (y/x). The vacancies of a minimum-size SV cluster in WxBy are energetically repulsive for y/x ⩽ 1, while the vacancies energetically bind together for y/x larger than 1. The formation energy of B Frenkel pairs in each W boride is lower than that of W Frenkel pairs. Among the six evaluated WxBy compositions, WB has the highest and the lowest formation energy of H interstitials and H-vacancy complexes, respectively; however, these two energies in WB2 are in reverse order. The average H binding energies to single vacancies in WB and WB4 are comparable with that in W, while this binding energy in WB3 and WB2 is obviously higher or lower than in W, respectively. The diffusion activation energy of H in W borides is anisotropic. One dimensional (1D) diffusion of H in W2B, as well as 1D/2D diffusion in WB, W2B5, WB3 and WB4 are preferred at relatively low temperatures; however, three-dimensional diffusion of H is predicted in WB2. The diffusion activation energy of H generally increases with B content in W borides due to the increasing local charge deficit caused by strong B–B covalent bonds. This study is useful for evaluating the performance of W borides in a fusion environment.