We review the present status of the QCD corrected cross sections and kinematic distributions for the production of a Higgs boson in association with bottom quarks at the Fermilab Tevatron and CERN Large Hadron Collider. Results are presented for the Minimal Supersymmetric Standard Model where, for large tan beta, these production modes can be greatly enhanced compared to the Standard Model case. The next-to-leading order QCD results are much less sensitive to the renormalization and factorization scales than the lowest order results, but have a significant dependence on the choice of the renormalization scheme for the bottom quark Yukawa coupling. We also investigate the uncertainties coming from the Parton Distribution Functions and find that these uncertainties can be comparable to the uncertainties from the remaining scale dependence of the next-to-leading order results. We present results separately for the different final states depending on the number of bottom quarks identified.