The distribution of fitness effects (DFE) is central to many questions in evolutionary biology.However, little is known about the differences in DFEs between closely related species. We use more than 9,000 coding genes orthologous onetoone across great apes, gibbons, and macaques to assess the stability of the DFE across great apes. We use the unfolded site frequency spectrum of polymorphic mutations (n = 8 haploid chromosomes per population) to estimate the DFE. We find that the shape of the deleterious DFE is strikingly similar across great apes. We confirm that effective population size (N e ) is a strong predictor of the strength of negative selection, consistent with the Nearly Neutral Theory.However, we also find that the strength of negative selection varies more than expected given the differences in N e between species. Across species, mean fitness effects of new deleterious mutations co varies with N e , consistent with positive epistasis among deleterious mutations. We find that the strength of negative selection for the smallest populations: bonobos and western chimpanzees, is higher than expected given their N e . This may result from a more efficient purging of strongly deleterious recessive variants in these populations. Forward simulations confirm that these findings are not artifacts of the way we are inferring N e and DFE parameters. All findings are replicated using only GCconservative mutations, thereby confirming that GCbiased gene conversion is not affecting our conclusions. 2 4 6 8 10 12 14 16All organisms undergo mutation. Studying the effect of mutations on fitness is fundamental to explain the patterns of genetic diversity within and between species and to predict the impact of population size on the probability of survival, or extinction of a species . The fitness effects of all new mutations that can occur in a given genome are described by the distribution of fitness effects (DFE). The allele frequency distributions or the site frequency spectrum (SFS) contains information to infer the DFE of new mutations. Current statistical methods infer the DFE while accounting for demography and other sources of distortion in the SFS (EyreWalker et al. 2006;Keightley and EyreWalker 2007;Boyko et al. 2008;Schneider et al. 2011;Galtier 2016;Kim et al. 2017;Tataru et al. 2017;Barton and Zeng 2018).The DFE of new deleterious amino acid mutations is assumed to be similar in closely related species, while the mean effect of those deleterious mutations (S d = 2N e s d ) is expected to increase as a function of the effective population size (N e ) (Kimura 1983; Ohta 1992). Very little is known about the variability in the DFE of new beneficial mutations across species as well as how this variability affects the estimates of the deleterious DFE. Most studies have described the DFE of new deleterious mutations of individual populations, while the study of the differences in the DFE between populations or species has remained unexplored. To our knowledge, Huber et al. (2017) is the first work that tries to ...