Mean fitness and inbreeding depression values in multi-locus models of the control of fitness were studied, using both a model of mutation to deleterious alleles, and a model of heterozygote advantage. Synergistic fitness interactions between loci were assumed, to find out if this more biologically plausible model altered the conclusions we obtained previously using a model of multiplicative interactions. Systems of unlinked loci were assumed. We used deterministic computer calculations, and approximations based on normal or Poisson theory. These approximations gave good agreement with the exact results for some regions of the parameter space. In the mutational model, we found that the effect of synergism was to lower the number of mutant alleles per individual, and thus to increase the mean fitness, compared with the multiplicative case. Inbreeding depression, however, was increased. Similar effects on mean fitness and inbreeding depression were found for the case of heterozygote advantage. For that model, the results were qualitatively similar to those previously obtained assuming multiplicativity. With the mutational load model, however, the mean fitness sometimes decreased, and the inbreeding depression increased, at high selfing rates, after declining as the selfing rate increased from zero. We also studied the behaviour of modifier alleles that changed the selfing rate, introduced into equilibrium populations. In general, the results were similar to those with the multiplicative model, but in some cases an ESS selfing rate, with selfing slightly below one, existed. Finally, we derive an approximate expression for the inbreeding depression in completely selfing populations. This depends only on the mutation rate and the dominance coefficient and can therefore be used to obtain estimates of the mutation rate to mildly deleterious alleles for plant species.