We present results on the isovector momentum fraction, hxi u−d , helicity moment, hxi Δu−Δd , and the transversity moment, hxi δu−δd , of the nucleon obtained using nine ensembles of gauge configurations generated by the MILC Collaboration using 2 þ 1 þ 1-flavors of dynamical highly improved staggered quarks. The correlation functions are calculated using the Wilson-Clover action, and the renormalization of the three operators is carried out nonperturbatively on the lattice in the RI 0-MOM scheme. The data have been collected at lattice spacings a ≈ 0.15, 0.12, 0.09, and 0.06 fm and M π ≈ 310, 220, and 135 MeV, which are used to obtain the physical values using a simultaneous chiral-continuum-finite-volume fit. The final results, in the MS scheme at 2 GeV, are hxi u−d ¼ 0.173ð14Þð07Þ, hxi Δu−Δd ¼ 0.213ð15Þð22Þ, and hxi δu−δd ¼ 0.208ð19Þð24Þ, where the first error is the overall analysis uncertainty and the second is an additional systematic uncertainty due to possible residual excited-state contributions. These results are consistent with other recent lattice calculations and phenomenological global fit values.