Photon-jet azimuthal correlations in proton-nucleus collisions are a promising tool for gaining information on the gluon distribution of the nucleus in the regime of nonlinear color fields. We compute such correlations from the process g → qqγ in the rapidity regime where both the projectile and target lightcone momentum fractions are small. By integrating over the phase space of the quark which emits the photon, subject to the restriction that the photon picks up most of the transverse momentum (to pass an isolation cut), we effectively obtain a g þ A → qγ process. For nearly back-to-back photon-jet configurations we find that it dominates over the leading-order process q þ A → qγ by two less powers of Q ⊥ =Q S , where Q ⊥ and Q S denote the net photon-jet pair momentum and the saturation scale of the nucleus, respectively. We determine the transverse-momentum-dependent gluon distributions involved in g þ A → qγ and the scale where they are evaluated. Finally, we provide analytic expressions for hcos nϕi moments, where ϕ is the angle between Q ⊥ and the average photon-jet transverse momentumP ⊥ , and first qualitative estimates of their transverse momentum dependence.