In this Letter, we present mean optical+NIR color gradient estimates for 5080 early-type galaxies (ETGs) in the grizYJHK wavebands of the Sloan Digital Sky Survey (SDSS) plus UKIRT Infrared Deep Sky Survey (UKIDSS). The color gradient is estimated as the logarithmic slope of the radial color profile in ETGs. With such a large sample size, we study the variation of the mean color gradient as a function of waveband with unprecedented accuracy. We find that (i) color gradients are mainly due, on average, to a metallicity variation of about -0.4dex per decade in galaxy radius; and (ii) a small, but significant, positive age gradient is present, on average, in ETGs, with the inner stellar population being slightly younger, by ~0.1dex per radial decade, than the outer one. Also, we show that the presence of a positive mean age gradient in ETGs, as found in the present study, implies their effective radius to be smaller at high z, consistent with observations.