Suppose we have nonindependent random variables $X \sim P$ and $Y \sim Q$, where $P$ and $Q$ denote their marginal distributions. We are interested in upper bounding $$ |\mathbf{E}_{X, Y\sim P \otimes Q } f(X, Y) - \mathbf{E}_{X, Y\sim P, Q } f(X, Y)| $$ for a given function $f$. Here $P \otimes Q$ denotes the product distribution of $P$ and $Q$, while $P, Q$ denotes their joint distribution (or equivalently a coupling).
Question: Can we upper bound such a mean difference using the Wasserstein distance between $P$ and $Q$, or the Wasserstein distance between $P \otimes Q$ and $P, Q$?
If needed, we are allowed to impose regularity conditions, such as Lipschitz continuity, on $f$.