© 2017 Cambridge University Press. The statistical properties are presented for the direct numerical simulation of a self-similar adverse pressure gradient (APG) turbulent boundary layer (TBL) at the verge of separation. The APG TBL has a momentum thickness-based Reynolds number range from to 13Â 800, with a self-similar region from to 12Â 300. Within this domain the average non-dimensional pressure gradient parameter , where for a unit density , with the displacement thickness, the mean shear stress at the wall and the far-field pressure gradient. This flow is compared with previous zero pressure gradient and mild APG TBL results of similar Reynolds number. All flows are generated via the direct numerical simulation of a TBL on a flat surface with far-field boundary conditions tailored to apply the desired pressure gradient. The conditions for self-similarity, and the appropriate length and velocity scales, are derived. The mean and Reynolds stress profiles are shown to collapse when non-dimensionalised on the basis of these length and velocity scales. As the pressure gradient increases, the extent of the wake region in the mean streamwise velocity profiles increases, whilst the extent of the log-layer and viscous sublayer decreases. The Reynolds stress, production and dissipation profiles of the APG TBL cases exhibit a second outer peak, which becomes more pronounced and more spatially localised with increasing pressure gradient. This outer peak is located at the point of inflection of the mean velocity profiles, and is suggestive of the presence of a shear flow instability. The maximum streamwise velocity variance is located at a wall normal position of of spanwise wavelength of . In summary as the pressure gradient increases the flow has properties less like a zero pressure gradient TBL and more akin to a free shear layer.