We present new (B,I) photometry for the globular cluster systems in eight Brightest Cluster Galaxies (BCGs), obtained with the ACS/WFC camera on the Hubble Space Telescope. In the very rich cluster systems that reside within these giant galaxies, we find that all have strongly bimodal color distributions All the BCGs show population gradients, with much higher relative numbers of red clusters within 5 kpc of their centers, consistent with their having formed at later times than the blue, metal-poor population. A striking new feature of the color distributions emerging from our data is that for the brightest clusters (M_I < -10.5) the color distribution becomes broad and less obviously bimodal. we suggest that it may be a characteristic of many BCGs. Furthermore, the blue (metal-poor) clusters become progressively redder with increasing luminosity, following a mass/metallicity scaling relation Z ~ M^0.55. We argue that these GCS characteristics are consistent with a hierarchical-merging formation picture in which the metal-poor clusters formed in protogalactic clouds or dense starburst complexes with gas masses in the range 10^7 - 10^10 M_Sun, but where the more massive clusters on average formed in bigger clouds with deeper potential wells where more pre-enrichment could occur.

We have used deep I-band (F814W) images from the HST archive to study the globular cluster systems around the brightest cluster galaxies (BCGs) in Abell 262, 3560, 3565 and 3742. Three of these BCGs have inner dust lanes and peculiar structural features which indicate past histories of low-level interaction and accretion. The deep I-band WFPC2 images have photometric limits which, for all four galaxies, reach near or just beyond the GCLF turnover point. Their specific frequencies are 8.24 \pm 1.65, 4.66 \pm 0.93, 2.58 \pm 0.52 and 2.62 \pm 0.52 respectively, all within a factor of two of the normal range for giant ellipticals. We obtain new estimates of the GCLF turnover magnitudes, which are shown to be consistent with an adopted Hubble constant of H_0 ~ 70 km s^{-1} Mpc^{-1} on the ``Hubble diagram'' of GCLF turnover apparent magnitude versus redshift, on a distance scale where the fundamental GCLF calibrator E galaxies (M87 and others) in Virgo are at d=16 Mpc.