Let $M$ and $N$ be smooth manifolds, and let $F: M \to N$ be a smooth map. Let $X$ and $Y$ be vector fields on $M$, and let $\tilde{X}$ and $\tilde{Y}$ be vector fields on $N$. We say that $X$ and $\tilde{X}$ are $F$-related if $dF_p(X_p) = \tilde{X}_{F(p)}$ for every $p \in M$. There is a theorem (naturality of the Lie bracket):
If $X$ and $\tilde{X}$ are $F$-related and $Y$ and $\tilde{Y}$ are $F$-related, then $[X,Y]$ and $[\tilde{X}, \tilde{Y}]$ are $F$-related.
A very concise proof is as follows:
Let $f \in C^\infty(N)$. We only need to show that $[X,Y](f \circ F) = [\tilde{X}, \tilde{Y}]f$, which is $$(XY - YX)(f \circ F) = (\tilde{X}\tilde{Y} - \tilde{Y}\tilde{X})f.$$ This equality follows from $XY(f \circ F) = X(\tilde{Y}f \circ F) = \tilde{X}\tilde{Y}f$.
My question is, how can we prove this theorem without introducing the smooth function $f$?
When $F$ is a local diffeomorphism, we can prove the theorem directly using coordinates. Let the matrix $(b^i_j) = (\partial_jF^i)$ with inverse $(a^i_j)$. We compute: \begin{align*} \tilde{x}^j \tilde{\partial}_j \tilde{y}^i - \tilde{y}^j \tilde{\partial}_j \tilde{x}^i &= b^j_k x^k a^l_j \partial_l(b^i_my^m) - b^j_k y^k a^l_j \partial_l(b^i_mx^m)\\ &= x^l\partial_l(b^i_my^m) - y^l\partial_l(b^i_mx^m)\\ &= b^i_m(x^l \partial_ly^m - y^l \partial_lx^m) + \partial_l\partial_mF^i(x^ly^m - x^my^l)\\ &= b^i_m(x^l \partial_ly^m - y^l \partial_lx^m). \end{align*}
This proof is also in the accepted answer of this MSE question: Pushforward of Lie Bracket.
However, when $F$ is not locally invertible, then the matrix $(b^i_j)$ won't be invertible either. The matrix $(a^i_j)$ in the proof above would be undefined! This situation seems very strange. It appears that the only way to prove this theorem is to use the abstract definition of $[X,Y]f = (XY - YX)f$. The coordinate expression $[X,Y] = (x^j \partial_j y^i - y^j \partial_j x^i)$ seems inadequate.
Another abstract definition of $[X,Y]$ is the Lie derivative $\mathcal{L}_XY$, which is computed using the flow map of $X$. We can also discuss the naturality of the Lie derivative: are $\mathcal{L}_XY$ and $\mathcal{L}_{\tilde{X}}\tilde{Y}$ $F$-related? If $F$ is not locally invertible, then it seems we can't prove this only using the abstract notion of Lie derivative either. We have to use the fact that $\mathcal{L}_XY = [X,Y]$ and return to the Lie bracket case. This also feels strange because the definition $[X,Y]f = (XY - YX)f$ is so special here.
I also discussed the question in this article on Banana Space (a Chinese website dedicated to sharing math): Intuitive interpretation of Lie bracket. However, it is written in Chinese.
Update: I seem to find a way to prove $\mathcal{L}_XY$ and $\mathcal{L}_{\tilde{X}}\tilde{Y}$ are $F$-related only using the abstract notion of Lie derivative. I have written an anwser.