Gianluca,
This works:
ParametricPlot3D[{x, y, x}, Element[{x, y}, Disk[]],
RegionFunction -> Function[{x, y, z}, x^2 + y^2 < 1]]
or better yet,
ParametricPlot3D[{x, y, x}, Element[{x, y}, Disk[]],
RegionFunction -> Function[{x, y, z}, Element[{x, y}, Disk[]]]]
but you should report that the element syntax no longer automatically limits the region.
Regards
Neil