R. Hirota and K. Kimura discovered integrable discretizations of the Euler and the Lagrange tops, given by birational maps. Their method is a specialization to the integrable context of a general discretization scheme introduced by W. Kahan and applicable to any vector field with a quadratic dependence on phase variables. We report several novel observations regarding integrability of the Kahan-Hirota-Kimura discretization. For several of the most complicated cases for which integrability is known (Clebsch system, Kirchhoff system, and Lagrange top), - we give nice compact formulas for some of the more complicated integrals of motion and for the density of the invariant measure, and - we establish the existence of higher order Wronskian Hirota-Kimura bases, generating the full set of integrals of motion. While the first set of results admits nice algebraic proofs, the second one relies on computer algebra.