$\newcommand{\M}{M}$ $\newcommand{\N}{N}$ $\newcommand{\TM}{TM}$ $\newcommand{\TN}{TN}$ $\newcommand{\TstarM}{T^*M}$ $\newcommand{\Ga}{\Gamma}$
Let $\M,\N$ be smooth manifolds, $\phi:\M \to \N$ be a smooth map. Let $\nabla$ be a symmetric connection on $\TN$, and let $X,Y \in \Ga(\TM)$.
Then the following holds: $$ \nabla^{\phi^*\TN}_X d\phi(Y) - \nabla^{\phi^*\TN}_Y d\phi(X) = d\phi([X,Y]) $$
where $\nabla^{\phi^*\TN}$ is the pullback connection induced on $\phi^*\TN$ by $\phi,\nabla$.
It is not hard to prove this by choosing local coordinates on $M,N$.
(This result comes up naturally in many different scenarios, for instance when calculating the first variation of the Dirichlet's integral)
Question:
Is there an "invariant" (coordinate-independent) proof of this result?
I do not think the "coordinates-proof" is simple enough. Perhaps there is an elegant way to show this result is indeed obvious at a glance?
The problem seems to be that the pullback-connection is defined by a characterising property which is local in nature (action on pullback sections), but I wonder if there is a way to bypass it somehow.