Genomic selection is a genetic evaluation of animals from their DNA, based on a huge number of markers covering the whole genome. It requires advanced approaches and in particular feature selection methods. Feature selection is a combinatorial problem that may be addressed by combinatorial optimization methods. We propose to combine an iterated local search (ILS) with a statistical evaluation of a multivariate regression and we compared three criteria in order to analyse their impact on the performance of the local search.