  We compute the complete next-to-leading order (NLO) SUSY-QCD corrections for
the associated production of a pseudoscalar Higgs boson with a bottom quark via
bottom-gluon fusion at the CERN Large Hadron Collider (LHC) and the Fermilab
Tevatron. We find that the NLO QCD correction in the MSSM reaches $40%\sim50%$
at the LHC and $45%\sim80%$ at the Tevatron in our chosen parameter space.
