Setup: Let $(\mathcal{X},d_{\mathcal{X}})$ and $(\mathcal{Y},d_{\mathcal{Y}})$ be two separable metric spaces. Let $M^1(\mathcal{X})$ be the space of Borel probability measures on $\mathcal{X}$ with finite first moment, i.e. a Borel probability measure $\mu$ on $\mathcal{X}$ is in $M^1(\mathcal{X})$ if $\int d_{\mathcal{X}}(x,o) d\mu(x)<\infty$ for any $o\in\mathcal{X}$. The space $M^1(\mathcal{Y})$ is defined in similar fashion.
Fix $\mu\in M^1(\mathcal{X})$ and $\nu\in M^1(\mathcal{Y})$ and define $$ d_\mu(x_1,x_2)= d_{\mathcal{X}}(x_1,x_2) -\int d_{\mathcal{X}}(x_1,x)\, d\mu(x) - \int d_{\mathcal{X}}(x_2,x)\, d\mu(x) + \int d_{\mathcal{X}}(x,x')\, d\mu^2(x,x'), $$ and a similar definition of $d_\nu:\mathcal{Y}\times \mathcal{Y}\to\mathbb{R}$.
Now let $S:L^2(\mathcal{X}\times \mathcal{Y},\mathcal{B}(\mathcal{X})\otimes \mathcal{B}(\mathcal{Y}),\mu\times \nu) \to L^2(\mathcal{X}\times \mathcal{Y},\mathcal{B}(\mathcal{X})\otimes \mathcal{B}(\mathcal{Y}),\mu\times \nu),$ be a Hilbert-Schmidt operator given by $$ S(f)(x,y) = \int d_\mu(x,x')d_\nu(y,y') f(x',y') d\mu\times \nu(x',y'). $$ and let $\{\lambda_i\}_{i\geq 1}$ denote the non-zero eigenvalues of $S$ repeated according to multiplicity.
Question: How do I prove the following identity: $$\sum_{i=1}^\infty\lambda_i=\int d_\mu(x,x)d_\nu(y,y) \, d\mu\times \nu(x,y).$$ I tried but failed to show that $S$ is of trace class, since they under certain conditions (which I also can't verify in this setup) satisfy that $$ Trace(S)=\int d_\mu(x,x)d_\nu(y,y) \, d\mu\times \nu(x,y), $$ which yields the result since $Trace(S)=\sum_{i=1}^\infty\lambda_i$ (if it were of trace class).
The identity $\sum_{i=1}^\infty\lambda_i=\int d_\mu(x,x)d_\nu(y,y) \, d\mu\times \nu(x,y)$ is stated in Distance covariance in metric spaces by Russell Lyons (2013) Theorem 2.7, without proof, so my approach with traces is only an idea.
If another way of proving the identity appears or if there is a counterexample, that would more than satisfy my needs. In the case of a counterexample i would very much appreciate stronger initial conditions rendering the identity true.
Please bear in mind that i am a novice in the theory of operators and trace class operators (only went down this road to explain the equality above), so references would be much appriciated.
Update: A counterexample to the operator being of trace-class is presented by Russell Lyons in the errata to the mentioned paper. Furthermore a proof that the formula holds and that the operator is of traceclass, whenever the marginal spaces posses additional nice proporties (isometric embeddability into hilbert spaces), is also presented in this errata.