Why must the Bogoliubov transform preserve anticommutation relations?
That's essentially what we mean by a quasiparticle: a degree of freedom that has been otherwise decoupled from the rest of the system and which behaves exactly like a particle as far as quantum mechanics can tell. Since the (anti)commutation relations are the crucial part of the quantum mechanics of the relevant particles, those need to be carried as well.
Moreover, the canonical anticommutation relations are the best you could hope for anyways. This is because the transformation, \begin{align} b_j &= \sum_{k=1}^n \left(U_{jk}c_k + V_{jk}c_k^\dagger\right)\\ b_j^\dagger &= \sum_{k=1}^n \left( U^*_{jk}c_k^\dagger + V_{jk}^*c_k\right) , \end{align} requires the anticommutator to read \begin{align} \{b_i^\dagger,b_j\} &= \left\{ \sum_{k=1}^n \left( U_{ik}^*c_k^\dagger + V_{ik}^*c_k\right), \sum_{l=1}^n \left(U_{jl}c_l + V_{jl}c_l^\dagger\right) \right\} \\&= \sum_{k=1}^n \sum_{l=1}^n \left[ U_{ik}^*U_{jl} \left\{c_k^\dagger,c_l \right\} + V_{ik}^*V_{jl} \left\{c_k,c_l^\dagger\right\} \right] \\&= \sum_{k=1}^n \sum_{l=1}^n \left[ -U_{ik}^*U_{jl} + V_{ik}^*V_{jl} \right] \delta_{kl} \\&= \sum_{k=1}^n\left[ V_{ik}^*V_{jk} - U_{ik}^*U_{jk} \right] \end{align} in the most general case. Thus, in the worst-case scenario, you've bungled up the (anti)commutation relations, and you're going to need a nonzero (anti)commutator any time you need to slide a $b_j$ past a $b_i^\dagger$, even if $i\neq j$. Thus, you require at the very least that the anticommutator be diagonal.
Moreover, since $\{b,b^\dagger\}=bb^\dagger + b^\dagger b$ is the sum of two positive semidefinite hermitian operators, the diagonal anticommutator needs to be real and non-negative, so you've got two possibilities:
- it can be positive, in which case it can be rescaled to one, or
- it can be zero, in which case you're essentially losing information, probably via a degenerate matrix, and the algebra spanned by the $b_i,b_i^\dagger$ will be unable to fully span all of the $c_i,c_i^\dagger$.
More broadly speaking, you're obliged to "conserve the quantumness" of the problem (unless you explicitly want to make mean-field approximations or similar) and this requires you to maintain a full set of generators for the operator algebra and therefore to have nontrivial anticommutators between your new quasiparticle fermionic operators.
Maybe we can think from an "opposite" point of view, here we have a binary Hamiltonian: (a single particle operator in the language of second quantization) $$H=c^{\dagger}Ac$$ the Bogoliubov transformation is essentially a transformation trying to diagonalize the Hamiltonian, or in other words, find a new basis. In quantum mechanics, we usually seek for two normalised bases that are related by a unitary transformation: $$\begin{align} |i\rangle_b&=U_{i,j}|j\rangle_c \\ \rightarrow b^{\dagger}_i &=U_{i,j} c^{\dagger}_j \\ \rightarrow b_i &= c_jU^{\dagger}_{j,i} \\ \end{align}$$ with $U_{i,j}U^{\dagger}_{j,i'} = \delta_{i,i'}$, and this automatedly leads to the fact that: $$\begin{align} [b_i,b^{\dagger}_j]_+ &=[c_{\alpha},c^{\dagger}_{\beta}]_+ U^{\dagger}_{\alpha,i}U_{j,\beta} \\ &=U^{\dagger}_{\alpha,i}U_{j,\alpha} \\ &=\delta_{i,j} \end{align}$$ so we can see that the commutation trlation does not change after the transformation.