One has the ordering of the three $\lambda$'s, i.e.
$$
\lambda_{\text{convex}} \leq \lambda_{\text{LSI}} \leq \lambda_{\text{SG}},
$$
where $\lambda_{\text{convex}}$ is the one from convexity, $\lambda_{\text{LSI}}$ is the one from the inverse Log-Sobolev constant (both as in your question) and $\lambda_{\text{SG}}$ the spectral gap, characterized via the smallest constant in the weighted Poincaré inequality
$$
\Vert f -1 \Vert_{L^2(\pi)}^2 \leq (2\lambda)^{-1} \int |\nabla f |^2 d\pi , \qquad \forall f : \int f \, d\pi = 1.
$$
here $f$ plays the role of $d\rho/ d\pi$. The first inequality is the HWI-inequality after Otto-Villani and the second is linearization of the LSI.
The equivalence $\lambda_{\text{convex}} = \lambda_{\text{LSI}} = \lambda_{\text{SG}}$ holds for $V$ being a positive quadratic form, i.e. $V(x) = x \cdot H x$ for a fixed symmetric positive matrix $H$. Then, the eigenvalues and eigenvectors of the according Ornstein-Uhlenbeck process are explicit (products of Hermite-polynomials) and obtained in dependence of the eigenvalues of $H$. The smallest non-negative eigenvalue is then the smallest eigenvalue of $H$, in accordance to the $\lambda$-convexity. The LSI is sandwiched inbetween anyways.
Except for this particular case, I expect that the equivalence breaks down for generic $V$ (non-quadratic).
This is easiest observed in the non-convex case. So let, $V$ be a double-well, with local max in $0$ and two minima in $~\pm 1$ and have convex quadratic growth outside a bounded region. Then, the Fokker-Planck evolution is $\lambda$-convex for a $\lambda<0$ being the lowest bound on the Hessian again, determined by the non-convexity of $V$ around the local maximum. However, the Log-Sobolev and spectral gap constants are still positive finite and can be obtained by combining the Bakry-Emery criterion with the Holley-Stroock perturbation principle, since $V$ is a bounded perturbation of a convex potential.
This becomes even more apparent, by considering the vanishing diffusion limit, i.e. consider for $\varepsilon>0$ the Fokker-Planck equation
$$
\partial_t \rho = \varepsilon \Delta \rho + \nabla \cdot(\rho \nabla V)
$$
In this case, it becomes clear that the time-scales captured by convexity and Log-Sobolev constants or spectral gaps are rather different.
From the comment of @Tobsn follows that, convexity measures local stability at every point in the space of probability measures and the setting of double-well potential it is readily checked that, one actually gets the opposite comparison of the type
$$
W(\rho_t, \hat\rho_t) \geq c e^{|\lambda| t} W(\delta_{-\eta},\delta_{\eta}) \quad\text{for } t \in [0, t_0]
$$
where $\eta>0$ is small and $t_0$ is also not too large. I write $|\lambda|=-\lambda$ to make clear, that the trajectories expand and do not converge.
This follows, because $\rho$ and $\hat\rho$ will follow mainly the deterministic ODE $\dot X_t = - \nabla V(X_t)$, which expands at rate $|\lambda|$ close to the local maximum. In particular, this result also holds for $\varepsilon=0$.
In comparison, the log-Sobolev constant and the inverse spectral gap measure a global time-scale quantifying the ergodicity in time, showing that the diffusion converges to the measure $\pi$ in entropy or the relative density $f_t = \rho_t / \pi$ in the $L^2(\pi)$-sense. Those are global averaged quantities and hence behave in general better, as the non-convex setting shows. However, note that for the $\varepsilon$-dependent case, both $\lambda_{LSI}$ and $\lambda_{SG}$ will degenerate to $0$ exponentially like $e^{-C/\varepsilon}$, which is a statement about metastability being present in this setting (the limiting ODE dynamic is non-ergodic, since it has actually three stationary states).
There is a theory of variable Ricci bounds, which can catch a bit this different local stability and might improve the gap between those constants. One, can probably think of the Bracamp-Lieb inequality as an instance of variable curvature, since it shows that $\text{Hess}\, V$ behaves like the local Ricci-tensor for the diffusion, that is
$$
\Vert f - 1 \Vert_\pi^2 \leq \int \left\langle \nabla f , \text{Hess} V\, \nabla f \right\rangle d\pi .
$$
Other ways, to make the different behaviour of the constants precise, at least in the scaling regime with vanishing diffusion $\varepsilon\ll 1$, are in 1d upper and lower estimates on the Log-Sobolev constant and spectral gap with the help of the Bobkov-Götze or Muckenhoupt criterion.
I can provide more details or references on the individual mentioned observations, but I'm not aware of a general result showing that for generic V (not perfectly quadratic potentials), there are strict inequalities between all of the three $\lambda$'s, i.e.
$$
\lambda_{\text{convex}} < \lambda_{\text{LSI}} < \lambda_{\text{SG}}.
$$
Best Answer
Yes this is true, formally this follows by the envelope theorem. In an abstract and very smooth setting, the envelope theorem says that for an objective functional depending on a parameter $t$ $$ F(t)=\max\limits_z f(t,z), $$ then the derivative of the optimal value can be computed as $$ \frac{dF}{dt}(t)=\partial_t f(t,z_t) \qquad \mbox{for any smooth selection of a maximizer }z_t \mbox{ of }F(t). $$ This can be seen easily: for any such choice of a maximizer, just apply a chain rule and use the optimality condition of $z_t$ in the maximization problem for fixed $t$: $$ F'(t)=\frac d{dt}f(t,z_t)=\partial_t f(t,z_t)+\underbrace{\partial_zf(t,z_t)}_{=0}\frac {dz_t}{dt}. $$ This means, roughly speaking, that one can simply forget that the minimizer varies, only the variation of the functional matter.
In your specific context, you are trying to differentiate (w.r.t $\rho$) the optimal value of the optimization problem given by the Kantorovich dual formulation $$ W^2(\rho,\eta) =F_\eta(\rho) =\max\limits_\phi f_\eta(\rho,\phi) =\max\limits_\phi \left\{\int \rho\phi+\int\eta\phi^c\right\} $$ (here $\eta$ is fixed once and for all, I'm mimicking my $F,f$ notations above to give some perspective and I hope the notation is sufficiently self-explanatory). Although the Kantorovich potential $\phi$ from $\rho$ to $\eta$ (the optimizer) varies when $\rho$ varies, the envelope theorem strongly suggests that you can actually argue as it did not vary at leading order (same for its $c$-transform $\phi^c$), and one can simply differentiate the functional w.r.t. the varying "parameter" $\rho$. Since the Kantorovich functional is linear in $\rho$, the conclusion is indeed that the first variation is given by $ \frac{\partial f_\eta}{\partial_\rho}(\rho,\phi)=\phi$. Of course various subtle problems may arise owing essentially to the infinite-dimensional setting and functional-analytic details, but this is the rough idea.
For a completely rigorous statement and proof I can recommend Filippo Santambrogio's book [1], in particular chapter 7 and Proposition 7.17
[1] Santambrogio, Filippo. "Optimal transport for applied mathematicians." Birkäuser, NY 55.58-63 (2015): 94.