  We investigate the properties of the hot gas in four fossil galaxy systems
detected at high significance in the Planck Sunyaev-Zeldovich (SZ) survey.
XMM-Newton observations reveal overall temperatures of kT ~ 5-6 keV and yield
hydrostatic masses M500,HE > 3.5 x 10e14 Msun, confirming their nature as bona
fide massive clusters. We measure the thermodynamic properties of the hot gas
in X-rays (out to beyond R500 in three cases) and derive their individual
pressure profiles out to R ~ 2.5 R500 with the SZ data. We combine the X-ray
and SZ data to measure hydrostatic mass profiles and to examine the hot gas
content and its radial distribution. The average Navarro-Frenk-White (NFW)
concentration parameter, c500 = 3.2 +/- 0.4, is the same as that of relaxed
`normal' clusters. The gas mass fraction profiles exhibit striking variation in
the inner regions, but converge to approximately the cosmic baryon fraction
(corrected for depletion) at R500. Beyond R500 the gas mass fraction profiles
again diverge, which we interpret as being due to a difference in gas clumping
and/or a breakdown of hydrostatic equilibrium in the external regions. Overall
our observations point to considerable radial variation in the hot gas content
and in the gas clumping and/or hydrostatic equilibrium properties in these
fossil clusters, at odds with the interpretation of their being old, evolved
and undisturbed. At least some fossil objects appear to be dynamically young.
