For a semilinear PDE, we usually have this FBSDE representation:
$\mathcal{X}_t=\mathcal{X}_0+\int^t_0 \mu (s,\mathcal{X}_s)\, ds\, +\int_0^t \sigma (s,\mathcal{X}_s)dW_s,\quad 0\leq t\leq T, \\ Y_t = g(\mathcal{X}_T)+ \int_t^T f(s, \mathcal{X}_s,Y_s,Z_s)ds - \int_t^TZ_sdW_s, \quad 0\leq t\leq T,$
where $W$ is a $d$-dimensional standard Brownian motion on some probability space $(\Omega,\mathbb{F},\mathbb{P})$ equipped with a filtration $\mathbb{F}=(\mathcal{F}_t)_t$, and $\mathcal{X}_0$ is an $\mathcal{F}_0$-measurable random variable valued in $\mathbb{R}^d$.
Now in this paper: https://arxiv.org/pdf/2006.01496.pdf on page 7, equation (2.11) they use this representation: $Y_t = g(\mathcal{X}_T) - \int_t^T f(s, \mathcal{X}_s,Y_s,Z_s)ds - \int_t^TZ_sdW_s, \quad 0\leq t\leq T,$
My question: How we get the first minus? I didn't find any literature/reference who stated this result. Can someone link me some reference or explain it? Any help is appreciated. Thanks.