Particle-hole symmetry is only antilinear in the one-particle space. It is linear and unitary when acting on the many-particle Fock space. See footnote after eq 4 in S.Ryu, A.Schnyder, A.Furusaki, A.Ludwig, Topological insulators and superconductors: ten-fold way and dimensional hierarchy New J. Phys. 12, 065010 (2010). ArXiv:0912.2157.
Supose that the one particle Hamiltonian has the property that
$$
C H^* C^{-1}= -H
$$
for some unitary matrix $C$. Then
$$
Hu_n=\lambda_n u_n\quad\Rightarrow \quad HCu_n^* = -\lambda_n Cu^*_n,
$$
so, when $\lambda$ is non zero, the single-particle eigenfunctions come in opposite-eigenvalue pairs. In the absence of zero energy states the ground state $|0\rangle$ has all negative-energy states occupied and is non-degenerate.
We define the action of a unitary particle-hole operator ${\mathsf C}$ on the many-body Fock space by
$$
{\mathsf C}\Psi_\beta {\mathsf C}^{-1}= \Psi_\alpha^\dagger C_{\alpha\beta}, \quad {\mathsf C}\Psi^\dagger_\beta {\mathsf C}^{-1}= C^\dagger_{\beta\alpha}\Psi_\alpha.
$$
When
${\mathsf C}$ acts on the Hamiltonian we have
$$
{\mathsf C}\hat H {\mathsf C}^{-1}=
{\mathsf C}\Psi^\dagger_\alpha H_{\alpha\beta} \Psi_\beta {\mathsf C}^{-1}\\
={\mathsf C}\Psi^\dagger_\alpha {\mathsf C}^{-1}{\mathsf C}H_{\alpha\beta}{\mathsf C}^{-1}{\mathsf C} \Psi_\beta {\mathsf C}^{-1}\\
= {\mathsf C}\Psi^\dagger_\alpha {\mathsf C}^{-1}H_{\alpha\beta}{\mathsf C} \Psi_\beta {\mathsf C}^{-1}\\
= C^\dagger_{\alpha\rho}\Psi_{\rho} H_{\alpha\beta} \Psi^\dagger_\sigma C_{\sigma\beta}\\
=- \Psi^\dagger_\sigma C_{\sigma\beta} H_{\alpha\beta}C^\dagger_{\alpha\rho}\Psi_{\rho}\\
=- \Psi^\dagger_\sigma C_{\sigma\beta} H^T_{\beta \alpha}C^\dagger_{\alpha\rho}\Psi_{\rho}\nonumber\\
=- \Psi^\dagger_\sigma C_{\sigma\beta} H^*_{\beta \alpha}C^\dagger_{\alpha\rho}\Psi_{\rho}\nonumber\\
=+ \Psi^\dagger_\sigma H_{\sigma\rho} \Psi_{\rho}.\\
= \hat H.
$$
So the one-particle transformation on $H$ leaves the many-particle hamiltonian invariant.
Note: In expressions like $ \Psi^\dagger_\sigma C_{\sigma\beta} H_{\alpha\beta}C^\dagger_{\alpha\rho}\Psi_{\rho}$, the quantities $\Psi^\dagger_\sigma$ and $\Psi_{\rho}$ are second quantized operators, whereas $C_{\sigma\beta}$ or $H_{\alpha\beta}$ are regular complex numbers (they are the components of the relevant matrices).
Note that we have used $C^{*\dagger} = C^T$ and the tracelessness (line 5 $\to$ 6) and hermiticity of $H$ in the above manipulations.
More importantly, and despite the appearance of ``$*$" in the action on $H$, the many-body operator ${\mathsf C}$ must act on the Fock space linearly:
$$
{\mathsf C}(\lambda |\psi_1\rangle+\mu |\psi_2\rangle)= \lambda {\mathsf C}|\psi_1\rangle+\mu {\mathsf C}|\psi_2\rangle.
$$
The linearity is required in the step
$$
{\mathsf C}H_{\alpha\beta}{\mathsf C}^{-1}= H_{\alpha\beta}.
$$