Following this question: Can we apply the continuous mapping theorem for the limiting joint distribution of the Tracy-Widom law?.
We know that
$$ \lim_{N\to\infty}P(N^{2/3}(\lambda_N-2)\le s_1,\dotsc, N^{2/3}(\lambda_{N-k+1}-2)\le s_k)=F_{\beta, k}(s_1,\dotsc, s_k). $$
Can we say that for any $\epsilon>0$, there exists a constant $c>0$ such that $$ P(N^{2/3}(\lambda_N-\lambda_{N-1})\le c)\ge 1-\epsilon? $$
Or assume that $(N^{2/3}(\lambda_N-2),\dotsc, N^{2/3}(\lambda_{N-k+1}-2))\to (Y_1,\dotsc, Y_{N-k+1})$ in distribution.
The question becomes that for any $\epsilon>0$, there exists a constant $c>0$ such that $$ P(Y_N-Y_{N-1}\le c)\ge 1-\epsilon? $$
We can say $N^{2/3}(\lambda_N-\lambda_{N-1})=O_p(1)$.