Search for Neutral Supersymmetric Higgs Bosons in $bbh$

Many extensions of the Standard Model predict greatly enhanced production rates of neutral Higgs bosons in association with bottom quarks. The Tevatron Run II data collected by the D\O\ detector are used to search for Higgs bosons in $bh(\rightarrow b\bar{b})$ and $b\bar{b}h(\rightarrow b\bar{b})$ channels that result to $3b$- jet and $4b$-jet final states, respectively. The status of the current data analysis, and expectations for Run II are presented.