Bogoliubov - de Gennes equation reads, $$\left( \begin{array}{cc} H_{0} - E_{F} & -i\sigma_{y}\Delta \\ i\sigma_{y}\Delta^{*} & E_{F} - H_{0}^{*} \end{array}\right) \left( \begin{array}{c} \psi_{e}^{\uparrow} \\ \psi_{e}^{\downarrow} \\ \psi_{h}^{\uparrow} \\ \psi_{h}^{\downarrow} \end{array}\right) = \mathcal{E} \left( \begin{array}{c} \psi_{e}^{\uparrow} \\ \psi_{e}^{\downarrow} \\ \psi_{h}^{\uparrow} \\ \psi_{h}^{\downarrow} \end{array}\right) $$
The above Hamiltonian which is $4$x$4$ obeys particle-hole symmetry, $$H = -\mathcal{C}H\mathcal{C}^{-1}$$ where $\mathcal{C} = \tau_{x}\mathcal{K}$, with $\tau_{x} = \sigma_{x} \oplus \sigma_{x}$ and complex conjugation operator $\mathcal{K}$. This symmetry implies that if You have solution $\Psi$ with energy $\mathcal{E}$, then You also have solution $\mathcal{C}\Psi$ with energy $-\mathcal{E}$. I always thought that if Hamiltonian posseses some symmetry then You have solution with the same energy. I am confused, cause I also read that holes are time-reversed electrons, $$\psi_{h} = \mathcal{T}\psi_{e}$$ with $\mathcal{T} = i\sigma_{y}\mathcal{K}$, but I can't show this from above BdG equation.