Massive early-type galaxies typically have two subpopulations of globular clusters (GCs) which often reveal radial colour (metallicity) gradients. Collating gradients from the literature, we show that the gradients in the metal-rich and metal-poor GC subpopulations are the same, within measurement uncertainties, in a given galaxy. Furthermore, these GC gradients are similar in strength to the stellar metallicity gradient of the host galaxy. At the very largest radii (e.g. greater than 8 galaxy effective radii), there is some evidence that the GC gradients become flat with near constant mean metallicity.Using stellarmetallicity gradients as a proxy,we probe the assembly histories of massive early-type galaxies with hydrodynamical simulations from the Magneticum suite of models. In particular, we measure the stellar metallicity gradient for the in situ and accreted components over a similar radial range as those observed for GC subpopulations. We find that the in situ and accreted stellar metallicity gradients are similar but have a larger scatter than the metal-rich and metal-poor GC subpopulations' gradients in a given galaxy. We conclude that although metal-rich GCs are predominately formed during the in situ phase and metal-poor GCs during the accretion phase of massive galaxy formation, they do not have a strict one-to-one connection.