Axial turbine stages for gas generator cycle type rocket engines typically employ highly-loaded supersonic stator vanes. Consequently, the flow pattern downstream of the vanes is characterized by shock waves which induce high-frequency excitation on the subsequent rotor. For this reason, the optimal design of the stator is crucial in the context of development of the next generation of high-performance rocket engines, where reusability is a principal design criterion. A thorough comprehension of the loss mechanisms combined with the adoption of automated optimisation techniques can therefore enable new stator designs that may provide large benefits in terms of overall turbine performance and lifespan. The scope of this study stems from these considerations and its objective is twofold, namely i) the shape optimisation of a su- personic stator for rocket engines and ii) the investigation of the loss mechanisms in supersonic axial turbine stator vanes at on- and off-design conditions. The investigation is performed on stator vanes that are under development for the first turbine stage of a gas generator cycle type rocket engine. The stator vanes are therefore optimised in order to reduce the profile losses by exploiting a novel adjoint optimisation framework for turbomachinery implemented in the open-source code SU2. The effect on the resulting flow field and loss sources is finally investigated. Results show that entropy based loss coefficient gains of 6% can be achieved via shape optimisation and that the fluid- dynamic performance of these vanes is less sensitive to changes in pressure ratio compared to the performance provided by the baseline configuration. Eventually, shock-waves remain the primary loss source.