Octonion algebra O has recently been used to study the fundamental physics of the Standard Model, such as its three-generation structure and its possibility of unifying gravity and quantum mechanics. Interestingly, this octonion algebra O has also been related to the SO(8) monopole and, consequently, links to the nine-dimensional MICZ-Kepler problem. This problem has been solved exactly by the variables separation method in three different coordinate systems, spherical, parabolic, and prolate spheroidal. In the present study, we establish a relationship between the variable-separation and the algebraic structure of SO(10) symmetry. Each of the spherical, parabolic, or prolate spheroidal bases is proved to be a set of eigenfunctions of a corresponding nonuplet of algebraically-independent integrals of motion. This finding also helps us to establish connections between the bases by the algebraic method. This connection, in turn, allows calculating a complicated integral of confluent Heun, generalized Laguerre, and generalized Jacobi polynomials, which may be engaging in analytics.