The Taylor-Couette flow between a stationary outer cylinder and rotating inner cylinder undergoes a supercritical transition. After becoming linearly unstable, the flow becomes progressively more complex: as the inner cylinder rotation Reynolds number Rein increases, the flow state changes to the Taylor vortex flow (TVF)- wavy Taylor vortex flow (WVF)- modulated wavy Taylor vortex flow (MWV). In contrast, annular Poiseuille flow, driven by an axial pressure gradient in concentric cylinders, undergoes a subcritical transition. Its subcritical turbulent flow features helical-shaped localized turbulence (HLT). The Taylor-Couette-Poiseuille flow, which is a combined shear flow of cylinder-rotation-driven flow and axial pressure-driven flow, is the subject of this study. We investigated the flow state transition processes for a high radius ratio of 0.883 at three different Rein values, using direct numerical simulations. We demonstrated that in the TVF and WVFbased cases, the pressure-driven axial flow stabilized into the Taylor-vortex-free flow field, with the WVF state transitioning to the TVF state before laminarization. A further increase in the axial pressure gradient led to intermittent turbulence, similar to HLT. These facts indicate that the switch from supercritical to subcritical transitions occurs across laminarization. In the MWV-based case, at a higher Rein, the flow does not exhibit laminarization but becomes fully turbulent, unlike in the lower Rein cases. However, the waviness of the Taylor vortex disappeared, and the pre-multiplied energy spectra confirmed partial stabilization before the transition to turbulence. From the perspective of Lumley's anisotropic invariant map, the TVF- and WVF-based cases have one- or two-component anisotropy under all conditions. However, the MWV-based case becomes continuously similar to the anisotropic map of typical turbulent channel flow as F(P) increases.