In this study, we investigate the behaviour of decaying compressible anisotropic turbulence using Direct Numerical Simulations (DNS). High-order schemes are employed to accurately solve the full three-dimensional Navier-Stokes equations. The study of compressibility effects is based on decay simulations of anisotropic turbulence with two parametric studies, one concerning the ratio between the compressible part and the total kinetic energy, and the other one concerning the turbulent Mach number. We focuse our attention on the return to isotropy phenomenon. We analyse our results in respect with previous studies