  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.
