How to prove Bianchi identity? \begin{align*} \varepsilon^{\mu\nu\rho\sigma}D_{\nu}F_{\rho\sigma}=0 \end{align*} using Jacobi identity; \begin{align*} \epsilon^{\mu\nu\rho\sigma}[D_{\mu},[D_{\rho},D_{\sigma}]]=0 \end{align*} where covarient derivative is given as \begin{align*} D_{\mu}=\partial_{\mu}-igA_{\mu} \end{align*}
I know the same question was asked on this site before; Bianchi identity of a non-Abelian gauge theory? In this answer, he used the fact that covariant derivation satisfies the Leibniz rule. So, I would like to know why this fact holds.