A Higgs boson produced in association with one or more bottom quarks is of great theoretical and experimental interest to the high-energy community. A precise prediction of its total and differential cross-section can have a great impact on the discovery of a Higgs boson with large bottom-quark Yukawa coupling, like the scalar (h^0 and H^0) and pseudoscalar (A^0) Higgs bosons of the Minimal Supersymmetric Standard Model (MSSM) in the region of large tanbeta. In this paper we apply the threshold resummation formalism to determine both differential and total cross-sections for b g to bPhi (where Phi = h^0, H^0), including up to next-to-next-to-next-to-leading order (NNNLO) soft plus virtual QCD corrections at next-to-leading logarithmic (NLL) accuracy. We present results for both the Fermilab Tevatron and the CERN Large Hadron Collider (LHC).