For a diffusion process $X(t)$ of drift $μ(x)$ and of diffusion coefficient $D=1/2$, we study the joint distribution of the two local times $A(t)= \int_{0}^{t} dτδ(X(τ)) $ and $B(t)= \int_{0}^{t} dτδ(X(τ)-L) $ at positions $x=0$ and $x=L$, as well as the simpler statistics of their sum $ Σ(t)=A(t)+B(t)$. Their asymptotic statistics for large time $t \to + \infty$ involves two very different cases : (i) when the diffusion process $X(t)$ is transient, the two local times $[A(t);B(t)]$ remain finite random variables $[A^*(\infty),B^*(\infty)]$ and we analyze their limiting joint distribution ; (ii) when the diffusion process $X(t)$ is recurrent, we describe the large deviations properties of the two intensive local times $a = \frac{A(t)}{t}$ and $b = \frac{B(t)}{t}$ and of their intensive sum $σ= \frac{Σ(t)}{t}=a+b$. These properties are then used to construct various conditioned processes $[X^*(t),A^*(t),B^*(t)]$ satisfying certain constraints involving the two local times, thereby generalizing our previous work [arXiv:2205.15818] concerning the conditioning with respect to a single local time $A(t)$. In particular for the infinite time horizon $T \to +\infty$, we consider the conditioning towards the finite asymptotic values $[A^*(\infty),B^*(\infty)]$ or $Σ^*(\infty) $, as well as the conditioning towards the intensive values $[a^*,b^*] $ or $σ^*$, that can be compared with the appropriate 'canonical conditioning' based on the generating function of the local times in the regime of large deviations. This general construction is then applied to the simplest case where the unconditioned diffusion is the Brownian motion of uniform drift $μ$.