Clearly all eigenvalues apart from the eigenvalue $\lambda(\epsilon)$ with $\lambda(0) = 1$ stay in the open unit disk for $\epsilon$ sufficiently small. To see what happens to the last eigenvalue, use eigenvalue first-order perturbation theory, for instance Theorem 1 in [Greenbaum, Li, and Overton - First-order Perturbation Theory for Eigenvalues and Eigenvectors](https://arxiv.org/abs/1903.00785): if $v$ and $z^*$ are the right and left eigenvalue associated to $\lambda(0)$, then the right and left eigenvectors of the $2\times 2$ matrix in the text are $\begin{bmatrix}v\\0\end{bmatrix}$ and $\begin{bmatrix}z^* & z^*(I-B)^{-1}\end{bmatrix}$ and one gets
$$
\frac{d\lambda(\epsilon)}{d\epsilon} = \frac{1}{z^*v} \begin{bmatrix}z^* & z^*(I-B)^{-1}\end{bmatrix}\begin{bmatrix}-I & 0\\-C & 0\end{bmatrix} \begin{bmatrix}v\\0\end{bmatrix} = -1 - \frac{1}{z^*v} z^*(I-B)^{-1}Cv.
$$
If this derivative is negative, for a sufficiently small $\epsilon > 0$ we have that $\lambda(\epsilon)$ is in the unit disk.

However, it should be possible to choose $z,v,B$ such that $z^*(I-B)^{-1}v < 0$, and $C = \alpha I$; then for sufficiently large $\alpha > 0$ the derivative is larger than $0$ and I don't think that the result holds.

For instance, $A = \begin{bmatrix}1 & 8 \\ 0 & 0.9\end{bmatrix}$, $B = \begin{bmatrix}1/4 & -1/4\\ -1/4 & 1/4\end{bmatrix}$, $C=I$ numerically gives me matrices $X$ with eigenvalues larger than $1$ for small $\epsilon$.