Let $(B_t)_{t \geq 0}$ a Brownian motion on $\mathbb R^n$ and $\ell_t \in \mathbb S$, where $\mathbb S$ is the space of all increasing càdlàg function from $(0,\infty)$ to $(0,\infty)$ with $\lim_{s\downarrow 0} = 0$. Then the quadratic variation of $(B_{\ell_t})_{t\geq 0}$ is given by $$[B_{\ell}]_t = \ell_t - \sum_{0<s\leq t}\Delta \ell_s + \sum_{0<s\leq t} | \Delta B_{\ell_s} |^2.$$
Heuristically, this is very clear: the square bracket can be written as the sum of a continuous and discontinuous part. Since the cts. part $[B_\ell]_t^c = l_t$ we get the first summand + the jumps which is the third summand. Finally, since $\ell_t$ is càdlàg, we have to subtracte the jumps $\Delta \ell_s$.
My thoughts so far: By Definition, the quadratic variation of a semimartingale $X$ is given by $$[X] = X^2 - \int X_- \ \mathrm dX.$$ We have to calculate the stochastic integral for $X_t := B_{\ell_t}$. Let $\Pi = \{0 \leq s_1^k \leq ... \leq s_{n(k)}^k \leq t\}$ be a partition of the interval $[0,T]$. Then $$\sum_i B_{\ell(s_j)} (B_{\ell(s_{j+1})} - B_{\ell(s_j)}) \overset{ucp}\longrightarrow \int_0^t B_{\ell_{t-}} \ \mathrm d B_{\ell_t}.$$
If there are reasonable sets A und B, then using the elementary equality $b^2 - a^2 - (b-a)^2 = 2a (b-a)$, we find \begin{align*} 2 \sum_i B_{\ell(s_j)} (B_{\ell(s_{j+1})} - B_{\ell(s_j)}) &= \sum_i (B_{\ell(s_{j+1})}^2 - B_{\ell(s_j)}^2) - \sum_i (B_{\ell(s_{j+1})} - \sum B_{\ell(s_j)})^2\\ &= \underbrace{ \sum_{i,B} (B_{\ell(s_{j+1})}^2 - B_{\ell(s_j)}^2) }_{= \ B_{\ell_t}} - \sum_{i,A} (B_{\ell(s_{j+1})}^2 - B_{\ell(s_j))^2 }\\ &\qquad - \Bigg( \underbrace{ \sum_{i,B} (B_{\ell(s_{j+1})} - B_{\ell(s_j)})^2 }_{\overset{ucp}{\longrightarrow} \ \ell_t} - \sum_{i,A} (B_{\ell(s_{j+1})} - B_{\ell(s_j)} )^2 \Bigg)\\ &\overset{k\to \infty}\longrightarrow B_{\ell_t}^2 - \sum_{0<s\leq t} (B_{\ell_s} - B_{\ell_{s-}} )^2 - (\ell_t - \sum_{0<s\leq t} (\ell_s - \ell_{s-})) \\ &= B_{\ell_t}^2 - \ell_t - \sum_{0<s\leq t} \Delta(B^2)_s + \sum_{0<s\leq t} \Delta \ell_s. \end{align*}
Is this decomposition of the sums possible? My idea was something like $$J(\epsilon) := \{s \in [0,t] : |\Delta X_s| > \epsilon\}.$$ Since $s \mapsto X_s$ is càdlàg, $J$ is a.s. finite, each $\epsilon > 0$ and $\sum_{s \in J(\epsilon)} \overset{\epsilon \to 0}\longrightarrow \sum_{0<s\leq t}$ a.s.
Is there an easier way? Thanks.