I give a short review of the theory of twisted symmetries of differential equations, emphasizing geometrical aspects. Some open problems are also mentioned.Dedicated to Juergen Scheurle on the occasion of his retirement 2 G. GAETA (a = 1, ..., m); partial derivatives will be denoted by u a J , where J is a multi-index J = {j 1 , ..., j n } of order |J| = j 1 + ... + j n and(here and somewhere in the following we moved the vector index of the x for typographical convenience). We denote by u (k) the set of all partial derivatives of order k, and by u [n] the set of all partial derivatives of order k ≤ n.2.1. Geometrical description of differential equations and solutions. The x are local coordinates in a manifold B, while u are local coordinates in a manifold U ; we consider the phase manifold M = B × U , which has a natural structure of bundle (M, π 0 , B) over B with fiber U . As well known [1,13,37,57,58,70] we can associate to (M, π 0 , B) its Jet bundles J k M of any order; these have a structure of fiber bundle (J k M, π k , B) over B (with projection π k ) but also of fiber bundle over those of lower order (with projectionWe also recall that the Jet bundle is equipped with a contact structure Ω [4,33,58,68,71]; this can be encoded in the contact forms 2 ω a J := du a J − u a J,i dx i . (k) f in (J k M, π k , B), which is just the set of points (x, u [k] ) with u a = f a (x) and u a J = f a J (x). Now u = f (x) is a solution to ∆ if and only if γ(k) f ⊂ S ∆ . 4 G. GAETAstart with a complete set of DIs of order zero and one, we can generate in this way a complete set of DIs of any order, basically because if U is q-dimensional, we have k · q DIs of order k (these includes those of lower orders), as follows at once from J k M being of dimension d k = (k + 1)q + 1. The situation is different for PDEs, as the dimension of J k M grows combinatorially; see e.g. the discussion in [57].2.3. Lie-point symmetries. A (Lie-point 3 ) symmetry, or more precisely a Liepoint symmetry generator, is a vector field X on M such that its prolongation X (k) to J k M is tangent to S ∆ . For a given ∆ in the form (5), this condition is expressed in local coordinates asIn these equations, also called determining equations, the F are given and one looks for ξ, ϕ (i.e. components of the vector field X) satisfying them. As the components ψ a J of X (k) along u J are given in terms of ξ, ϕ and their derivatives, all dependencies of u a J (with |J| = 0) are fully explicit, and hence (12) decouple into a system of simpler equations, one for each monomial in the u a J ; each of these is a linear PDE for the ξ and ϕ, and they can be solved algorithmically -usually with the help of a symbolic manipulation program, as the dimension of the system can be quite large. This fails in the case of first order equations.