Since $F_{ij} = \epsilon_{ijk}B^k$ one has $\epsilon^{lij} F_{ij} = \epsilon^{lij}\epsilon_{ijk}B^k = 2\delta_{lk} B^k$, hence $B^l = \frac{1}{2} \epsilon^{lij} F_{ij}$.
The third Maxwell's equation in OP's question can be expressed with the field strength tensor $F_{\mu\nu}$ according to
\begin{align}
0 &= \epsilon^{ijk} \partial_j E_k + \partial_0 B^i \\
&= -\epsilon^{ijk} \partial_j F_{k0} + \frac{1}{2} \epsilon^{ikl} \partial_0 F_{kl}\\
&= -2\epsilon^{ijk} \partial_j F_{k0} + \epsilon^{ijk} \partial_0 F_{jk}\\
&=-\epsilon^{ijk} \partial_0 F_{jk} + \epsilon^{ijk} \partial_j F_{k0} - \epsilon^{ijk} \partial_j F_{0k}\ .
\end{align}
The Bianchi identity, i.e., the relation $\partial_\mu F_{\alpha\beta}+ \partial_\alpha F_{\beta\mu}+ \partial_\beta F_{\mu\alpha} =0$, can be compactly arrange with a four-dimensional totally antisymmetric Levi-Civita symbol $$\epsilon^{\alpha\beta\gamma\delta} \partial_\beta F_{\gamma\delta}=0$$
with the property $\epsilon^{0\beta\gamma\delta}=\epsilon^{ijk}$.
Replacing the three-dimensional Levi-Civita symbol we obtain
\begin{align}
0&=-\epsilon^{ijk} \partial_0 F_{jk} + \epsilon^{ijk} \partial_j F_{k0} - \epsilon^{ijk} \partial_j F_{0k}\\
&=\epsilon^{i0jk} \partial_0 F_{jk} + \epsilon^{ij0k} \partial_j F_{k0} + \epsilon^{ijk0} \partial_j F_{0k}\\
&=\epsilon^{i\beta\gamma\delta} \partial_\beta F_{\gamma\delta}\ .
\end{align}
The still missing case $\alpha=0$ follows immediately from the last Maxwell's equation:
\begin{align}
0&=\partial_i B^i\\
&=\epsilon^{ijk} \partial_i F_{jk}\\
&=\epsilon^{0\beta\gamma\delta}\partial_\beta F_{\gamma\delta}\ ,
\end{align}
which eventually reduces to the Bianchi identity.
PS: note the sign for covariant and contravariant indices, e.g. $F_{0i}=\eta_{00}\eta_{ij}F^{0j}=-F^{0i}$.