Density functional theory has been used to study lithium intercalation into TiO2(B) at low to moderate concentrations (0 < x(Li) ≤ 0.25) with a range of density functionals: LDA, GGA (PW91, PBE, PBEsol), and GGA+U (PBE+U, PBEsol+U), with the GGA+U calculations employing a Hubbard +U correction to the Ti d states. LDA and GGA functionals give the same general behaviour, whereas qualitatively different behaviour is predicted by GGA+U for electronic structure and the order of stability of occupied intercalation sites. LDA/GGA functionals predict LixTiO2(B) to be metallic, with the excess charge distributed over all the Ti sites. In contrast, GGA+U predicts defect states in the band gap corresponding to charge strongly localised at specific Ti sites. All the considered functionals predict A1 and/or A2 site occupation at x(Li)=0.25, which challenges the interpretation of previous neutron data that at this composition the C site is preferentially occupied.