We investigate the structural relaxation of a soft-sphere liquid quenched isochorically (φ = 0.7) and instantaneously to different temperatures Tf above and below the glass transition. For this, we combine extensive Brownian dynamics simulations and theoretical calculations based on the non-equilibrium self-consistent generalized Langevin equation (NE-SCGLE) theory. The response of the liquid to a quench generally consists of a sub-linear increase of the α-relaxation time with system's age. Approaching the ideal glass-transition temperature from above (Tf > Ta ) sub-aging appears as a transient process describing a broad equilibration crossover for quenches to nearly arrested states. This allows us to empirically determine an equilibration timescale t eq(Tf) that becomes increasingly longer as Tf approaches Ta . For quenches inside the glass (Tf less than or equal to Ta) the growth rate of the structural relaxation time becomes progressively larger as Tf decreases and, unlike the equilibration scenario, τα remains evolving within the whole observation time-window. These features are consistently found in theory and simulations with remarkable semi-quantitative agreement, and coincide with those revealed in a previous and complementary study [Phys. Rev. 96, 022608 (2017)] that considered a sequence of quenches with fixed final temperature Tf = 0 but increasing φ towards the hard-sphere dynamical arrest volume fraction φ a HS = 0.582. The NESCGLE analysis, however, unveils various fundamental aspects of the glass transition, involving the abrupt passage from the ordinary equilibration scenario to the persistent aging effects that are characteristic of glass-forming liquids.