Sorry for the delay, I have been away for a few days. This calculation is horrible and I don't guarantee that this is the best way to do this. I just did it in coordinates and it seems to work.
We want to show that
$$f_*[X,Y]=[f_*(X),f_*(Y)]$$
We will show this by showing that both sides act the same way on functions on $N$ and so must be the same vector field.
Let $h:N\to\mathbb R$, we have by definition
$$f_*([X,Y])(h)=[X,Y](h\circ f)$$
We use coordinates $x^i$ for $M$ and $y^m$ for $N$. We have
\begin{equation}
\begin{aligned}f_*([X,Y])(h)&=[X,Y](h\circ f)(x)\\
&=X^i\frac{\partial}{\partial x^i}\left(Y^j\frac{\partial}{\partial x^j}(h\circ f)\right)(x)-Y^i\frac{\partial}{\partial x^i}\left(X^j\frac{\partial}{\partial x^j}(h\circ f)\right)(x)\\
&=X^i\frac{\partial}{\partial x^i}\left(Y^j\frac{\partial h}{\partial y^m}(f(x))\frac{\partial y^m}{\partial x^i}(x)\right)-Y^i\frac{\partial}{\partial x^i}\left(X^j\frac{\partial h}{\partial y^m}(f(x))\frac{\partial y^m}{\partial x^i}(x)\right)\\
&=X^i\frac{\partial Y^j}{\partial x^i}\frac{\partial h}{\partial y^m}(f(x))\frac{\partial y^m}{\partial x^i}(x)-Y^i\frac{\partial X^j}{\partial x^i}\frac{\partial h}{\partial y^m}(f(x))\frac{\partial y^m}{\partial x^i}(x)
\end{aligned}
\end{equation}
While going from the other end
\begin{equation}
\begin{aligned}
\phantom a[f_*(X),f_*(Y)](h)&=f_*(X)(f_*(Y)(h))-f_*(X)(f_*(Y)(h))\\
&=f_*(X)\left(Y^j\frac{\partial y^m}{\partial x^j}(x)\frac{\partial h}{\partial y^m}(f(x))\right)-f_*(Y)\left(X^j\frac{\partial y^m}{\partial x^j}(x)\frac{\partial h}{\partial y^m}(f(x))\right)\\
&=X^iY^j\frac{\partial y^m}{\partial x^j}(x)\frac{\partial y^n}{\partial x^i}(x)\frac{\partial^2 h}{\partial y^m\partial y^n}(f(x))+X^i\frac{\partial Y^j}{\partial x^m}\frac{\partial y^m}{\partial x^j}(x)\frac{\partial x^m}{\partial y^n}(f(x))\frac{\partial y^n}{\partial x^i}(x)\frac{\partial h}{\partial y^m}(f(x))\\
&\phantom=-Y^iX^j\frac{\partial y^m}{\partial x^j}(x)\frac{\partial y^n}{\partial x^i}(x)\frac{\partial^2 h}{\partial y^m\partial y^n}(f(x))-Y^i\frac{\partial X^j}{\partial x^m}\frac{\partial y^m}{\partial x^j}(x)\frac{\partial x^m}{\partial y^n}(f(x))\frac{\partial y^n}{\partial x^i}(x)\frac{\partial h}{\partial y^m}(f(x))\\
&=X^i\frac{\partial Y^j}{\partial x^i}\frac{\partial y^m}{\partial x^j}(x)\frac{\partial h}{\partial y^m}(f(x))-Y^i\frac{\partial X^j}{\partial x^i}\frac{\partial y^m}{\partial x^j}(x)\frac{\partial h}{\partial y^m}(f(x))\\
\end{aligned}
\end{equation}
We have now shown that both sides are the same. In this derivation I have used the chain rule as well as the Jacobian, some cancellations and the fact that $\frac{\partial y^m}{\partial x^j}(f(x))\frac{\partial x^k}{\partial y^m}(x)=\delta^k_j$.