“…Finding similarity matrices for each isospectral matrix L in Propositions 2, 4, 5, 6, 7, 8 would provide a tool for construction of additional algebraic identities satisfied by the zeros of the appropriate polynomials. A method that compares the spectral and the pseudospectral matrix representations of linear differential operators has been employed in [4,7,16] to construct new and re-obtain known algebraic identities satisfied by classical, Krall and Sonin-Markov orthogonal polynomials. It would be interesting to adapt the last method to the cases of the generalized hypergeometric, generalized basic hypergeometric, Wilson and Racah as well as Askey-Wilson and q-Racah polynomials and to compare the results with those reviewed in this chapter as well as to utilize the new developments to study properties of the solvable nonlinear systems (45), (68), (87), (100), (117), (127).…”