The distribution of fitness effects (DFE) of new mutations plays a fundamental role in evolutionary genetics. However, the extent to which the DFE differs across species has yet to be systematically investigated. Furthermore, the biological mechanisms determining the DFE in natural populations remain unclear. Here, we show that theoretical models emphasizing different biological factors at determining the DFE, such as protein stability, back-mutations, species complexity, and mutational robustness make distinct predictions about how the DFE will differ between species. Analyzing amino acid-changing variants from natural populations in a comparative population genomic framework, we find that humans have a higher proportion of strongly deleterious mutations than Drosophila melanogaster. Furthermore, when comparing the DFE across yeast, Drosophila, mice, and humans, the average selection coefficient becomes more deleterious with increasing species complexity. Last, pleiotropic genes have a DFE that is less variable than that of nonpleiotropic genes. Comparing four categories of theoretical models, only Fisher's geometrical model (FGM) is consistent with our findings. FGM assumes that multiple phenotypes are under stabilizing selection, with the number of phenotypes defining the complexity of the organism. Our results suggest that long-term population size and cost of complexity drive the evolution of the DFE, with many implications for evolutionary and medical genomics. T he distribution of fitness effects (DFE) represents the distribution of selection coefficients, s, of random mutations in the genome. Here, s quantifies the relative change in fitness due to the mutation. The DFE plays a fundamental role in evolutionary genetics because it quantifies the amount of deleterious, neutral, and adaptive mutations entering a population (1). Despite the importance and considerable study of the DFE (1-3), the extent to which the DFE in terms of s differs across species has yet to be systematically quantified. Furthermore, the biological factors determining the DFE in different species remain elusive. Several theoretical models propose different mechanisms for the evolution of the DFE (Fig. 1) (4-8). Although each of these models has a reasonable theoretical basis as well as some support from experimental evolution studies or microbial studies, which model best explains differences in the DFE between species has not yet been determined. Nor have these models been tested with genetic variation data from natural populations in higher organisms. Because experimental evolution studies in laboratory organisms often use a homogeneous environment and genetically homogeneous organisms, they may better satisfy some of the assumptions of these theoretical models. However, natural populations may provide different qualitative results due to increased resolution to measure weakly deleterious mutations and the unnatural selection pressure in the laboratory (2, 9). Importantly, five theoretical models for the evolution of the DFE predic...