Two aims are pursued in this article. The first one is methodological and consists of demonstrating the symmetry calculation techniques based on the commutation relation on the example of the quasi-linear heat equation. The second one consists of investigating the following open question: is the algebra of higher symmetries of the linear heat equation exhausted by those which may be obtained by means of the classical symmetries and well-known recursion operators? A positive solution of this question is given.