We study magnetohydrodynamical (MHD) turbulence excited by magnetorotational instability in differentially rotating flows with cylindrical shearing box simulations. Compared with the simulation in a local Cartesian shearing box, the cylindrical simulation exhibits considerably larger intermittency, although the time-averaged field strengths are comparable in both cases. The large temporal variation is manly due to the compressible amplification of magnetic fields. The radial gradient of epicyclic frequency, \(\kappa \) , which cannot be handled in the normal Cartesian model, is the cause of the compressible amplification. The spatial variation of \(\kappa \) also causes continuous and ubiquitous formation of localized high-density and weak-field zonal-flow structures.