We consider steady advection-diffusion about a slender ( 1) body of revolution at arbitrary O(1) Péclet numbers (Pe). The transported scalar attenuates at large distances, and is governed by axisymmetric (either Dirichlet or Neumann) data prescribed at the body boundary. The advecting field is assumed to be an axisymmetric Stokes flow approaching a uniform stream at large distances and satisfying impermeability at the boundary; otherwise, the interfacial distribution of tangential velocity is assumed arbitrary, irrotational and no-slip Stokes flows being particular cases. Employing the method of matched asymptotic expansions, we develop a systematic scheme for calculating the scalar concentration in increasing powers of ln −1 (1/ ). The leading term in the inner expansion coincides with the pure diffusion case, the second term depends nonlinearly on the magnitude of the far-field stream, and higher-order terms depend on the boundary distribution of tangential velocity. In the special case of irrotational flow and Neumann boundary conditions the logarithmic expansion terminates, leaving an algebraic error in . The general formulae developed can be directly applied to numerous physical scenarios. We here consider the classical problem of forced heat convection from an isothermal body, finding a two-term expansion for Nu( , Pe)/Nu( , 0), the ratio of the Nusselt number to its value at Pe = 0. This ratio is insensitive to the particle shape at the asymptotic orders considered; at moderately large Pe ( −1 ) its deviation from unity is O [ln(Pe)/ ln(1/ )], marking the poor effectiveness of advection about slender bodies. The expansion is compared to a numerical computation in the case of a prolate spheroid in both irrotational and no-slip Stokes flows.