This is an addition to the Laurent's answer. First, it should be said that if one derives all the functors, one will get an isomorphim --- $Ri_*M \otimes^L Ri_*N \cong Ri_* (M \otimes^L N)$. Indeed, it is a simple corollary of the projection formula:
$$
Ri_*M \otimes^L Ri_*N \cong
Ri_*(M \otimes^L i^*Ri_*N) \cong
Ri_*(M \otimes^L N)
$$
(the second isomorphism is by the flat base change). What goes wrong with the underived version is that $M$ and $N$ have higher direct images which then have $Tor$'s all of which eventually get canceled. In the particular example of Laurent one has
$$
R^1i_*O_{Y'} = y_1^{-1}y_2^{-1}k[y_1^{-1},y_2^{-1}],
\quad
R^1i_*O_{Z'} = z_1^{-1}z_2^{-1}k[z_1^{-1},z_2^{-1}],
$$
where are $y_1,y_2$ are coordinates on $Y$ and $z_1,z_2$ are coordinates on $Z$.
In addition to $O_Y\otimes O_Z = k$ we have
$$
Tor_2(O_Y,R^1i_*O_{Z'}) = Tor_2(R^1i_*O_{Y'},O_Z) = k,
$$
$$
Tor_4(R^1i_*O_{Y'},R^1i_*O_{Z'}) = k,
$$
and it is easy to see that all this cancels in the spectral sequence calculating $Ri_*O_{Y'} \otimes^L Ri_*O_{Z'}$.