Single-locus equilibrium frequencies of a partially recessive deleterious mutation under the mutation–selection balance model are derived for partially selfing autotetraploid populations. Assuming multiplicative fitness interactions among loci, approximate solutions for the mean fitness and inbreeding depression values are also derived for the multiple locus case and compared with expectations for the diploid model. As in diploids, purging of deleterious mutations through consanguineous matings occurs in autotetraploid populations, i.e. the equilibrium mutation load is a decreasing function of the selfing rate. However, the variation of inbreeding depression with the selfing rate depends strongly on the dominance coefficients associated with the three heterozygous genotypes. Inbreeding depression can either increase or decrease with the selfing rate, and does not always vary monotonically. Expected issues for the evolution of the selfing rate consequently differ depending on the dominance coefficients. In some cases, expectations for the evolution of the selfing rate resemble expectations in diploids; but particular sets of dominance coefficients can be found that lead to either complete selfing or intermediate selfing rates as unique evolutionary stable state.