Within a generalized parton model approach, with inclusion of spin and intrinsic transverse momentum effects, we show how the latest, highly precise, midrapidity data on the transverse single spin asymmetry measured in $pptopi^0, X$ by the PHENIX Col
laboration at RHIC [1], can be used to get a first estimate on the still poorly known gluon Sivers distribution. To this end we also adopt the present information on the quark Sivers functions, as extracted from semi-inclusive deeply inelastic scattering data. This analysis updates a previous study by some of us where a first bound on this distribution was obtained [2].
