In generic curved spacetimes, the unavailability of a natural choice of vacuum state introduces a serious ambiguity in the Fock quantization of fields. In this review, we study the case of fermions described by a Dirac field in non-stationary spacetimes, and present recent results obtained by us and our collaborators about well-motivated criteria capable to ensure the uniqueness in the selection of a vacuum up to unitary transformations, at least in certain situations of interest in cosmology. These criteria are based on two reasonable requirements. First, the invariance of the vacuum under the symmetries of the Dirac equations in the considered spacetime. These symmetries include the spatial isometries. Second, the unitary implementability of the Heisenberg dynamics of the annihilation and creation operators when the curved spacetime is treated as a fixed background. This last requirement not only permits the uniqueness of the Fock quantization but, remarkably, it also allows us to determine an essentially unique splitting between the phase space variables assigned to the background and the fermionic annihilation and creation variables. We first consider Dirac fields in 2 + 1 dimensions and then discuss the more relevant case of 3 + 1 dimensions, particularizing the analysis to cosmological spacetimes with spatial sections of spherical or toroidal topology. We use this analysis to investigate the combined, hybrid quantization of the Dirac field and a flat homogeneous and isotropic background cosmology when the latter is treated as a quantum entity, and the former as a perturbation. Specifically, we focus our study on a background quantization along the lines of loop quantum cosmology. Among the Fock quantizations for the fermionic perturbations admissible according to our criteria, we discuss the possibility of further restricting the choice of a vacuum by the requisite of a finite fermionic backreaction and, moreover, by the diagonalization of the fermionic contribution to the total Hamiltonian in the asymptotic limit of large wave numbers of the Dirac modes. Finally, we argue in support of the uniqueness of the vacuum state selected by the extension of this diagonalization condition beyond the commented asymptotic region, in particular proving that it picks out the standard Poincaré and Bunch–Davies vacua for fixed flat and de Sitter background spacetimes, respectively.