HERA data on diffractive DIS show deviations from twist 2 DGLAP predictions below $Q^2sim 5$ GeV$^2$ at low pomeron $xi$, which may reach up to 100%. These deviations are consistent with higher twists effects extracted from the saturation model. It i
s a first direct evidence for the higher twists in DIS. This finding affects determination of the diffractive parton densities that are used for the predictions at the LHC.
