We derive a relatively simple and Lorentz-invariant expression for a Wigner function of a paraxial vortex electron described as a Laguerre-Gaussian wave packet.While relativistic generalizations of a Wigner function [1] of a quantum system have been shown to exist [2-5], these functions may still suffer from a lack of Lorentz invariance [4]. This is all the more so for wave packets with the non-Gaussian spatial profiles, such as a vortex state with an orbital angular momentum (OAM) ℓ [6], an Airy beam [7,8], and their different generalizations. Although there are several works where the Wigner functions of such packets are derived (see, for instance, [9-12]) and even studied experimentally for twisted photons [11,13], the problem of invariance of these functions under the Lorentz transformations for massive wave packets (first of all for vortex electrons, as they have recently been obtained experimentally [6]) remains unsolved.One of the reasons why it has not been done for vortex electrons yet is that the very wave functions describing these states (the so-called Bessel beams and the Laguerre-Gaussian packets [6]) are usually written in terms of the non-invariant quantities and the transformation properties of these functions remain unclear. As noted by , this lack of explicit invariance per se does not violate correct transformation properties of the observables. Nevertheless, aiming at relativistic applications in particle and hadronic physics [14], it is desirable to know how the Wigner functions of the wave packets behave under the Lorentz boosts. To achieve this goal, the packets themselves need to be described in terms of the variables with definite transformation laws, which has been done for Gaussian packets of the massive neutrinos in [15,16] and recently for vortex electrons in [17].Here we derive a Wigner function of a massive particle with a phase vortex (say, the twisted electron) described via the paraxial Laguerre-Gaussian packet. The latter is Lorentz invariant for longitudinal boosts (along a propagation axis) and represents a massive generalization of the corresponding state of a twisted photon. Correct relativistic transformation of all variables was made possible by using an invariant approach to describe wave packets adopted from [15,16]. As we argue in [17], it is convenient to choose the condition of paraxiality for massive particles to be invariant too, in contrast to the twisted photons. In this paper, we imply that the packet is wide in x-space compared to a Compton wavelength, σ ⊥ ≫ λ c = /mc ≡ 1/m, or that it is narrow in p-space, σ ∼ 1/σ ⊥ ≪ m. It is this Lorentz invariant condition of paraxiality that significantly simplifies the calculations and the resultant expression for the Wigner function. Indeed,