This review delves into the utilization of a sextic oscillator within the β degree of freedom of the Bohr Hamiltonian to elucidate critical-point solutions in nuclei, with a specific emphasis on the critical point associated with the β shape variable, governing transitions from spherical to deformed nuclei. To commence, an overview is presented for critical-point solutions E(5), X(5), X(3), Z(5), and Z(4). These symmetries, encapsulated in simple models, all model the β degree of freedom using an infinite square-well (ISW) potential. They are particularly useful for dissecting phase transitions from spherical to deformed nuclear shapes. The distinguishing factor among these models lies in their treatment of the γ degree of freedom. These models are rooted in a geometrical context, employing the Bohr Hamiltonian. The review then continues with the analysis of the same critical solutions but with the adoption of a sextic potential in place of the ISW potential within the β degree of freedom. The sextic oscillator, being quasi-exactly solvable (QES), allows for the derivation of exact solutions for the lower part of the energy spectrum. The outcomes of this analysis are examined in detail. Additionally, various versions of the sextic potential, while not exactly solvable, can still be tackled numerically, offering a means to establish benchmarks for criticality in the transitional path from spherical to deformed shapes. This review extends its scope to encompass related papers published in the field in the past 20 years, contributing to a comprehensive understanding of critical-point symmetries in nuclear physics. To facilitate this understanding, a map depicting the different regions of the nuclide chart where these models have been applied is provided, serving as a concise summary of their applications and implications in the realm of nuclear structure.