Scaling and intermittency in compressible isotropic turbulence at the turbulent Mach number Mt ranging from 0.5 to 1.0 are studied by using numerical simulations with solenoidal forcing. Linear relations between the structure functions of the compressible velocity component and those of thermodynamic variables are modeled based on the shock jump conditions and are verified by numerical simulations. At a turbulent Mach number around 1.0, the relative scaling exponent of the structure functions saturates with an increase of the order. After proper normalization, the tails of the probability density functions (PDFs) of the increments of the compressible velocity component and thermodynamic variables overlap one another for different separations. Moreover, we study the conditional PDFs of the increments with respect to the shocklet. Linear relations between the tails of unconditional PDFs and conditional PDFs are established. The shocklet plays an important role in the determination of the PDF tails. The compressible velocity increment is decomposed into a negative component and a positive component. The negative component of the compressible velocity increment exhibits a scaling behavior with the saturation of the scaling exponent at high orders, which is similar to the Burgers turbulence, while the positive component of the compressible velocity increment exhibits a power-law scaling behavior, which is similar to the incompressible turbulence.