Another way
\begin{equation}
X_t^2=X_0^2-2\alpha \int_{0}^{t}X^2_s ds+2\sigma\int_0^t X_s\sqrt{X_s} dW_s+\sigma^2\int_{0}^{t} X_s ds
\end{equation}
thus
$$\mathbb{E}[X_t^2]=X_0^2-2\alpha\int_{0}^{t}\mathbb{E}[X^2_s] ds+\sigma^2\int_{0}^{t} \mathbb{E}[X_s] ds$$
As you saw in the first answer,
$\mathbb{E}[X_s]=X_0^2e^{-\alpha s}$,
thus
$$\mathbb{E}[X_t^2]=X_0^2-2\alpha\int_{0}^{t}\mathbb{E}[X^2_s] ds-\frac{\sigma^2}{\alpha}X_0^2\left(e^{-\alpha t}-1\right) $$
Set $m(t)=\mathbb{E}[X_t^2]$, we have
$$m(t)=X_0^2-2\alpha \int_0^t m(s)ds-\frac{\sigma^2}{\alpha}X_0\left(e^{-\alpha t}-1\right)$$
Take differentiate with respect to time,
$$m'(t)=-2\alpha\,m(t)+\sigma^2 X_0\,e^{-\alpha t}$$
In other words
$$m'(t)+2\alpha\,m(t)=\sigma^2 X_0\,e^{-\alpha t}$$
This ODE is a First-order equation ,
$$\mu(t)=e^{\int {2\alpha } dt}=e^{2\alpha t}$$
and
$$m(t)=\frac{1}{e^ {2\alpha t}}\left(\sigma^2 X_0 \int e^ {2\alpha t}e^ {-\alpha t}dt+c\right)$$
where $c$ is a constant, we have
$$m(t)=e^ {-2\alpha t}\left(\frac{\sigma^2 X_0}{\alpha}e^{\alpha t} +c\right)$$
on the other hand $m(0)=X_0^2$ and $m(0)=\frac{\sigma^2 X_0}{\alpha}$, therefor $c=X_0^2-\frac{\sigma^2 X_0}{\alpha}$. Finall we have
$$m(t)=\frac{\sigma^2}{\alpha}X_0e^{-\alpha t}+\left(X^2_0-\frac{\sigma^2}{\alpha}X_0\right)e^{-2\alpha t}$$