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 Re-delta 2 = 570 to 13 800, with a self-similar region from Re-delta 2 = 10 000 to 12 300. Within this domain the average non-dimensional pressure gradient parameter beta = 39, where for a unit density beta = delta P-1'(e)/tau(w), with delta(1) the displacement thickness, tau(w) the mean shear stress at the wall and P'(e) the far-field pressure gradient. This flow is compared with previous zero pressure gradient and mild APG TBL (beta = 1) 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 delta(1) of spanwise wavelength of 2 delta(1). 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.