We study the large order behaviour in perturbation theory of the Bjorken, Ellis-Jaffe and Gross-Llewellyn-Smith sum rules. In particular, we consider their first infrared renormalons, for which we obtain their analytic structure with logarithmic accuracy and also an approximate determination of their normalization constant. Estimates of higher order terms of the perturbative series are given. The Renormalon subtracted scheme is worked out for these observables and compared with experimental data. Overall, good agreement with experiment is found. This allows us to obtain {\hat a}_0 and some higher-twist non-perturbative constants from experiment: {\hat a}_0=0.141\pm 0.089; f_{3,RS}(1 GeV)=-0.124^{+0.137}_{-0.142} GeV^2.