We have performed a spectroscopic study of globular cluster (GC) system associated with the Virgo cD galaxy M87 using the Subaru/FOCAS MOS mode. We derive ages, metallicities and abundance ratios from the GC spectra using Simple Stellar Population (SSP) models. The metallicity distribution function (MDF) obtained empirically based on Milky Way GCs is consistent with the MDF derived from SSP models. A comparison with a meta-analysis using literature data sample of 15 other GC systems shows good agreement with our results. The properties of GCs acquired from the spectra will be used to test the recent theoretical prediction of a significant inflection along the colour-metallicity relations (Yoon et al. 2006). If confirmed, the non-linearity of the relations would shed new light on the interpretation of the GC colour bimodality. The robustness of our results is being tested against the choice of a SSP model, measurement errors and sample selection towards the goal of better understanding the formation history of GCs and host galaxy.