Modelling stellar populations in galaxies is a key approach to gain knowledge on the still elusive process of galaxy formation as a function of cosmic time. In this review, after a summary of the state-of-art, I discuss three aspects of the modelling, that are particularly relevant to massive galaxies, the focus of this symposium, at low and high-redshift. These are the treatment of the Thermally-Pulsating Asymptotic Giant Branch phase, evidences of an unusual Initial Mass Function, and the effect of modern stellar libraries on the model spectral energy distribution.