If Mutual Information is $\text{KLD}(P(X, Y) || P(X)P(Y))$, why is $\text{KLD}(P(X)P(Y) || P(X, Y))$ never mentioned/used? Mutual Information is defined as the Kullback-Liebler Divergence between $P(X, Y)$ and $P(X)P(Y)$. Using definitions of KLD, this can be intuitively understood as the cost in average code length if we obtain data from the $P(X, Y)$ distribution but instead encode it optimally assuming the $P(X)P(Y)$ distribution. It can also be thought of as the expected information gain from samples to reject the null hypothesis of $P(X, Y)$ to the alternative $P(X)P(Y)$.
However, why don't we ever discuss/use the expected information gain from samples to reject the null hypothesis of $P(X)P(Y)$ to the alternative of $P(X, Y)$? The KL Divergence is not symmetric, and thus $\text{KLD}(P(X)P(Y) || P(X, Y))$ also defines a divergence in which $X$ and $Y$ are interchangeable. Why is this metric not used as much (never?) in literature?
Also, the MI has a convenient relation to the (conditional) entropy: $H(X) = H(X \mid Y) + I(X;Y)$ as seen in this figure taken from Information Theory and its Relation to Machine Learning on arXiv.org:

Does such a relation exist for the reverse KLD described above?

If the full distribution $P(X, Y)$ is not computationally tractable, then we may choose to work with a simpler distribution $P(X)P(Y)$. In this case,  $\text{KLD}(P(X, Y) || P(X)P(Y))$ will tell us how well the factored distribution approximates the full distribution we are actually interested in. If we are able to work with the full distribution of interest, there's usually no reason to see how well the full version approximates some simpler factored version.
However, the reverse direction does come up in variational Bayesian inference, where we want to learn a distribution over some number of unobserved variables, call them $Y$ and $Z$, but only observe $X$. Variational Bayes is used when the posterior of interest $P(Y, Z | X)$ is computationally intractable. (Typically it's intractable because we have many unobserved variables, rather than just the two $Y$ and $Z$). In order to formulate an optimization problem that is tractable, variational Bayes introduces an approximating distribution over unobserved variables $Q(Y, Z)$ by defining a lower bound on the log marginal probability of the observed data:
\begin{align}
\log P(x) & = \log \sum_{y, z} P(x, y, z) \\
& = \log \sum_{y, z} Q(y, z) \frac{P(x, y, z)}{Q(y, z)}\\
& \geq \sum_{y, z} Q(y, z) \log \left(\frac{P(x, y, z)}{Q(y, z)}\right) \\
& = E_Q\left[P(x, y, z)\right] - E_Q\left[Q(y, z)\right]
\end{align}
Where the inequality is due to Jensen's inequality. We can then learn $Q$ by maximizing this lower bound on the marginal:
\begin{align}
Q^*
& = \underset{Q}{\text{argmax}} E_Q\left[P(x, y, z)\right] - E_Q\left[Q(y, z)\right]
\end{align}
We can see that this is equivalent to minimizing $\text{KLD}(Q || P)$ by using negation to change the maximization problem into a minimization problem and adding $\log P(x)$ (which doesn't affect the optimum since $Q$ doesn't depend on $\log P(x)$):
\begin{align}
Q^* & = \underset{Q}{\text{argmin}}~ -\left(E_Q\left[P(x, y, z)\right] - E_Q\left[Q(y, z)\right]\right)\\
  & = \underset{Q}{\text{argmin}}~
    E_Q\left[Q(y, z)\right] - E_Q\left[P(x, y, z)\right] + \log P(x) \\
  & = \underset{Q}{\text{argmin}}~
    E_Q\left[Q(y, z)\right] - E_Q\left[\log P(y, z | x)\right]\\
  & = \underset{Q}{\text{argmin}}~ \text{KLD}(Q(y, z) || P(y, z | x))
\end{align}
To make the learning problem tractable, $Q$ is typically taken to be from a constrained family of distributions. A popular choice is the "mean-field" approximation that severs all dependencies: $Q(y, z) = Q(y)Q(z)$. In this case, we are then finding $Q$ to minimize $\text{KLD}(Q(y)Q(z) || P(y, z | z))$.
