The efficacy of selection may increase or decrease with selfing depending upon the recombination environment.

Much theory has focused on how a population's selfing rate affects the ability of natural selection to remove deleterious mutations from a population. However, most such theory has focused on mutations of a given dominance and fitness effect in isolation. It remains unclear how selfing affects the purging of deleterious mutations in a genome-wide context where mutations with different selection and dominance coefficients co-segregate. Here, we use individual-based forward simulations and analytical models to investigate how mutation, selection and recombination interact with selfing rate to shape genome-wide patterns of mutation accumulation and fitness. In addition to recovering previously described results for how selfing affects the efficacy of selection against mutations of a given dominance class, we find that the interaction of purifying selection against mutations of different dominance classes changes with selfing and recombination rates. In particular, when recombination is low and recessive deleterious mutations are common, outcrossing populations transition from purifying selection to pseudo-overdominance, dramatically reducing the efficacy of selection. At these parameter combinations, the efficacy of selection remains low until populations hit a threshold selfing rate, above which it increases. In contrast, selection is more effective in outcrossing than (partial) selfing populations when recombination rates are moderate to high and recessive deleterious mutations are rare.