Let's forget about normalization of the two sides of the inequality as if $\operatorname{rank}(A) \leq \mathrm{rank}(B)$ then $\mathrm{rank}(c A) \leq \mathrm{rank}(d B)$ for any $c,d \neq 0$.
Now let's write $I_n = \rho + (I-\rho)$. As $\rho \leq I_n$ because $\rho$ is a quantum state, we have $(I-\rho) \geq 0$. Then
$$
\begin{aligned}
\Phi(I_n) &= \Phi(\rho + (I_n-\rho)) \\
&= \Phi(\rho) + \Phi(I_n-\rho)
\end{aligned}
$$
Thus $\mathrm{rank}(\Phi(I_n)) = \mathrm{rank}(\Phi(\rho) + \Phi(I_n-\rho))$. But as $\Phi(\rho)$ and $\Phi(I_n - \rho)$ are both positive-semidefinite we have $\mathrm{rank}(\Phi(\rho) + \Phi(I_n-\rho)) \geq \mathrm{rank}(\Phi(\rho))$ (see this answer for a proof of that).
Putting it all together we get
$$
\mathrm{rank}(\Phi(I_n)) \geq \mathrm{rank}(\Phi(\rho) )
$$
which implies the idenity in the question.
Regarding the title of the question
You should not interpret this result as saying that the maximally mixed state is always mapped to a state of maximal rank in the output system. For example, you can take take a channel $\Phi(\rho) = \mathrm{Tr}(\rho) |0\rangle \langle 0 |$ whose output will always have rank $1$ regardless of the dimension of the output sysytem.