lm.gls() is to fit Linear Models by Generalized Least Squares estimation method

GLM() is for Generalized linear models. GLM can use different link functions to fit different models including linear models