We prove L p → L q estimates for local maximal operators associated with dilates of codimension two spheres in Heisenberg groups; these are sharp up to two endpoints. The results can be applied to improve currently known bounds on sparse domination for global maximal operators. We also consider lacunary variants, and extensions to Métivier groups.