We study $L^p(μ)$ estimates for the commutator $[H,b]$, where the operator $H$ is a dyadic model of the classical Hilbert transform introduced in \cite{arXiv:2012.10201,arXiv:2212.00090} and is adapted to a non-doubling Borel measure $μ$ satisfying a dyadic regularity condition which is necessary for $H$ to be bounded on $L^p(μ)$. We show that $\|[H, b]\|_{L^p(μ) \rightarrow L^p(μ)} \lesssim \|b\|_{\mathrm{BMO}(μ)}$, but to {\it characterize} martingale BMO requires additional commutator information. We prove weighted inequalities for $[H, b]$ together with a version of the John-Nirenberg inequality adapted to appropriate weight classes $\widehat{A}_p$ that we define for our non-homogeneous setting. This requires establishing reverse Hölder inequalities for these new weight classes. Finally, we revisit the appropriate class of nonhomogeneous measures $μ$ for the study of different types of Haar shift operators.