Transitive closure is nothing but $M_{R^*}$ as you know it. And
$M_{R^∗} = M_R ∨ M^{[2]}_R ∨ M^{[3]}_R ∨···∨ M^{[n]}_R .$
So, $n^2(2n-1)(n-1)$ was for computing $M_R^{[i]}$ , $i = 1 \,\,to\,\,n$.
And Now, We need to Join(Union) these matrices to compute $M_{R^*}$ .. So, For that $n^2$ elements in $M_{R^*}$ and each needs $(n-1)$ boolean addition. So, that extra $n^2(n-1)$