Nor Maximum Likelihood Estimators (MLE) neither the Moments Matching Estimators (MME) for the two parameters $\alpha, \beta$ (shape and rate respectively) are unbiased.
Is there a closed formula to compute an unbiased estimator for both the shape and rate of the gamma distribution?