I think the easiest way to see this is through the fact that the trace is independent of the basis in which we take it. That is, for any operator $A$ and any orthonormal basis $\{|i\rangle\}_i$ we can write
$$
\mathrm{Tr}[A] = \sum_i \langle i | A | i \rangle.
$$
So if $|\psi_j\rangle$ is normalised we can just choose any orthonormal basis which includes $|\psi_j\rangle$, let's write such a basis as $\{|\psi_j\rangle, |\psi_{j,1}^{\perp}\rangle, |\psi_{j,2}^{\perp}\rangle, \dots\}$. Then we have
$$
\begin{aligned}
\mathrm{Tr}[|\psi_j\rangle \langle\psi_j|M] &= \langle\psi_j|\psi_j\rangle\langle\psi_j|M|\psi_j\rangle + \sum_{i} \langle\psi_{j,i}^{\perp}|\psi_j\rangle\langle\psi_j|M|\psi_{j,i}^{\perp}\rangle \\
&= \langle\psi_j|\psi_j\rangle\langle\psi_j|M|\psi_j\rangle \\
&= \langle\psi_j|M|\psi_j\rangle,
\end{aligned}
$$
where on the second line we used the fact that $\langle\psi_{j,i}^{\perp}|\psi_j\rangle = 0$ for all $i$ as they form an orthonormal basis and on the final line we used normalization $\langle\psi_j|\psi_j\rangle = 1$.