[This](https://www.rdocumentation.org/packages/mgcv/versions/1.9-1/topics/anova.gam) seems like the appropriate function to test both for smooth effects (e.g., age) and for effects of regressors.