Genetic drift promotes and recombination hinders speciation on holey fitness landscapes.

Dobzhansky and Muller proposed a general mechanism through which microevolution, the substitution of alleles within populations, can cause the evolution of reproductive isolation between populations and, therefore, macroevolution. As allopatric populations diverge, many combinations of alleles differing between them have not been tested by natural selection and may thus be incompatible. Such genetic incompatibilities often cause low fitness in hybrids between species. Furthermore, the number of incompatibilities grows with the genetic distance between diverging populations. However, what determines the rate and pattern of accumulation of incompatibilities remains unclear. We investigate this question by simulating evolution on holey fitness landscapes on which genetic incompatibilities can be identified unambiguously. We find that genetic incompatibilities accumulate more slowly among genetically robust populations and identify two determinants of the accumulation rate: recombination rate and population size. In large populations with abundant genetic variation, recombination selects for increased genetic robustness and, consequently, incompatibilities accumulate more slowly. In small populations, genetic drift interferes with this process and promotes the accumulation of genetic incompatibilities. Our results suggest a novel mechanism by which genetic drift promotes and recombination hinders speciation.

Empirical Evidence That Complexity Limits Horizontal Gene Transfer.

Horizontal gene transfer (HGT) is a major contributor to bacterial genome evolution, generating phenotypic diversity, driving the expansion of protein families, and facilitating the evolution of new phenotypes, new metabolic pathways, and new species. Comparative studies of gene gain in bacteria suggest that the frequency with which individual genes successfully undergo HGT varies considerably and may be associated with the number of protein-protein interactions in which the gene participates, that is, its connectivity. Two nonexclusive hypotheses have emerged to explain why transferability should decrease with connectivity: the complexity hypothesis (Jain R, Rivera MC, Lake JA. 1999. Horizontal gene transfer among genomes: the complexity hypothesis. Proc Natl Acad Sci U S A. 96:3801-3806.) and the balance hypothesis (Papp B, Pál C, Hurst LD. 2003. Dosage sensitivity and the evolution of gene families in yeast. Nature 424:194-197.). These hypotheses predict that the functional costs of HGT arise from a failure of divergent homologs to make normal protein-protein interactions or from gene misexpression, respectively. Here we describe genome-wide assessments of these hypotheses in which we used 74 existing prokaryotic whole genome shotgun libraries to estimate rates of horizontal transfer of genes from taxonomically diverse prokaryotic donors into Escherichia coli. We show that 1) transferability declines as connectivity increases, 2) transferability declines as the divergence between donor and recipient orthologs increases, and that 3) the magnitude of this negative effect of divergence on transferability increases with connectivity. These effects are particularly robust among the translational proteins, which span the widest range of connectivities. Whereas the complexity hypothesis explains all three of these observations, the balance hypothesis explains only the first one.

Recombination drives the evolution of mutational robustness.

Recombination can impose fitness costs as beneficial parental combinations of alleles are broken apart, a phenomenon known as recombination load. Computational models suggest that populations may evolve a reduced recombination load by reducing either the likelihood of recombination events (bring interacting loci in physical proximity) or the strength of interactions between loci (make loci more independent of one another). We review evidence for each of these possibilities and their consequences for the genotype-fitness relationship. In particular, we expect that reducing interaction strengths between loci will lead to genomes that are also robust to mutational perturbations, but reducing recombination rates alone will not. We note that both mechanisms most likely played a role in the evolution of extant populations, and that both can result in the frequently-observed pattern of physical linkage between interacting loci.

Population structure promotes the evolution of costly sex in artificial gene networks.

We build on previous observations that Hill-Robertson interference generates an advantage of sex that, in structured populations, can be large enough to explain the evolutionary maintenance of costly sex. We employed a gene network model that explicitly incorporates interactions between genes. Mutations in the gene networks have variable effects that depend on the genetic background in which they appear. Consequently, our simulations include two costs of sex-recombination and migration loads-that were missing from previous studies of the evolution of costly sex. Our results suggest a critical role for population structure that lies in its ability to align the long- and short-term advantages of sex. We show that the addition of population structure favored the evolution of sex by disproportionately decreasing the equilibrium mean fitness of asexual populations, primarily by increasing the strength of Muller's Ratchet. Population structure also increased the ability of the short-term advantage of sex to counter the primary limit to the evolution of sex in the gene network model-recombination load. On the other hand, highly structured populations experienced migration load in the form of Dobzhansky-Muller incompatibilities, decreasing the effective rate of migration between demes and, consequently, accelerating the accumulation of drift load in the sexual populations.

Evolutionary pathways to NS5A inhibitor resistance in genotype 1 hepatitis C virus.

Direct-acting antivirals (DAAs) targeting NS5A are broadly effective against hepatitis C virus (HCV) infections, but sustained virological response rates are generally lower in patients infected with genotype (gt)-1a than gt-1b viruses. The explanation for this remains uncertain. Here, we adopted a highly accurate, ultra-deep primer ID sequencing approach to intensively study serial changes in the NS5A-coding region of HCV in gt-1a- and gt-1b-infected subjects receiving a short course of monotherapy with the NS5A inhibitor, elbasvir. Low or undetectable levels of viremia precluded on-treatment analysis in gt-1b-infected subjects, but variants with the resistance-associated substitution (RAS) Y93H in NS5A dominated rebounding virus populations following cessation of treatment. These variants persisted until the end of the study, two months later. In contrast, while Y93H emerged in multiple lineages and became dominant in subjects with gt-1a virus, these haplotypes rapidly decreased in frequency off therapy. Substitutions at Q30 and L31 emerged in distinctly independent lineages at later time points, ultimately coming to dominate the virus population off therapy. Consistent with this, cell culture studies with gt-1a and gt-1b reporter viruses and replicons demonstrated that Y93H confers a much greater loss of replicative fitness in gt-1a than gt-1b virus, and that L31M/V both compensates for the loss of fitness associated with Q30R (but not Y93H) and also boosts drug resistance. These observations show how differences in the impact of RASs on drug resistance and replicative fitness influence the evolution of gt-1a and gt-1b viruses during monotherapy with an antiviral targeting NS5A.

Permissivity of Dipeptidyl Peptidase 4 Orthologs to Middle East Respiratory Syndrome Coronavirus Is Governed by Glycosylation and Other Complex Determinants.

Middle East respiratory syndrome coronavirus (MERS-CoV) utilizes dipeptidyl peptidase 4 (DPP4) as an entry receptor. While bat, camel, and human DPP4 support MERS-CoV infection, several DPP4 orthologs, including mouse, ferret, hamster, and guinea pig DPP4, do not. Previous work revealed that glycosylation of mouse DPP4 plays a role in blocking MERS-CoV infection. Here, we tested whether glycosylation also acts as a determinant of permissivity for ferret, hamster, and guinea pig DPP4. We found that, while glycosylation plays an important role in these orthologs, additional sequence and structural determinants impact their ability to act as functional receptors for MERS-CoV. These results provide insight into DPP4 species-specific differences impacting MERS-CoV host range and better inform our understanding of virus-receptor interactions associated with disease emergence and host susceptibility.IMPORTANCE MERS-CoV is a recently emerged zoonotic virus that is still circulating in the human population with an ∼35% mortality rate. With no available vaccines or therapeutics, the study of MERS-CoV pathogenesis is crucial for its control and prevention. However, in vivo studies are limited because MERS-CoV cannot infect wild-type mice due to incompatibilities between the virus spike and the mouse host cell receptor, mouse DPP4 (mDPP4). Specifically, mDPP4 has a nonconserved glycosylation site that acts as a barrier to MERS-CoV infection. Thus, one mouse model strategy has been to modify the mouse genome to remove this glycosylation site. Here, we investigated whether glycosylation acts as a barrier to infection for other nonpermissive small-animal species, namely, ferret, guinea pig, and hamster. Understanding the virus-receptor interactions for these DPP4 orthologs will help in the development of additional animal models while also revealing species-specific differences impacting MERS-CoV host range.

The emergence of performance trade-offs during local adaptation: insights from experimental evolution.

Environmental heterogeneity is considered a general explanation for phenotypic diversification, particularly when heterogeneity causes populations to diverge via local adaptation. Performance trade-offs, such as those stemming from antagonistic pleiotropy, are thought to contribute to the maintenance of diversity in this scenario. Specifically, alleles that promote adaptation in one environment are expected to promote maladaptation in alternative environments. Contrary to this expectation, however, alleles that underlie locally adaptive traits often fail to exhibit fitness costs in alternative environments. Here, we attempt to explain this paradox by reviewing the results of experimental evolution studies, including a new one of our own, that examined the evolution of trade-offs during adaptation to homogeneous versus heterogeneous environments. We propose that when pleiotropic effects vary, whether or not trade-offs emerge among diverging populations will depend critically on ecology. For example, adaptation to a locally homogeneous environment is more likely to occur by alleles that are antagonistically pleiotropic than adaptation to a locally heterogeneous environment, simply because selection is blind to costs associated with environments that are not experienced locally. Our literature review confirmed the resulting prediction that performance trade-offs were more likely to evolve during selection in homogeneous than heterogeneous environments. The nature of the environmental heterogeneity (spatial versus temporal) and the length of the experiment also contributed in predictable ways to the likelihood that performance trade-offs evolved.

An Evolving Genetic Architecture Interacts with Hill-Robertson Interference to Determine the Benefit of Sex.

Sex is ubiquitous in the natural world, but the nature of its benefits remains controversial. Previous studies have suggested that a major advantage of sex is its ability to eliminate interference between selection on linked mutations, a phenomenon known as Hill-Robertson interference. However, those studies may have missed both important advantages and important disadvantages of sexual reproduction because they did not allow the distributions of mutational effects and interactions (i.e., the genetic architecture) to evolve. Here we investigate how Hill-Robertson interference interacts with an evolving genetic architecture to affect the evolutionary origin and maintenance of sex by simulating evolution in populations of artificial gene networks. We observed a long-term advantage of sex-equilibrium mean fitness of sexual populations exceeded that of asexual populations-that did not depend on population size. We also observed a short-term advantage of sex-sexual modifier mutations readily invaded asexual populations-that increased with population size, as was observed in previous studies. We show that the long- and short-term advantages of sex were both determined by differences between sexual and asexual populations in the evolutionary dynamics of two properties of the genetic architecture: the deleterious mutation rate ([Formula: see text]) and recombination load ([Formula: see text]). These differences resulted from a combination of selection to minimize [Formula: see text] which is experienced only by sexuals, and Hill-Robertson interference experienced primarily by asexuals. In contrast to the previous studies, in which Hill-Robertson interference had only a direct impact on the fitness advantages of sex, the impact of Hill-Robertson interference in our simulations was mediated additionally by an indirect impact on the efficiency with which selection acted to reduce [Formula: see text].

Evolutionary rescue and the coexistence of generalist and specialist competitors: an experimental test.

Competition for resources is thought to play a critical role in both the origins and maintenance of biodiversity. Although numerous laboratory evolution experiments have confirmed that competition can be a key driver of adaptive diversification, few have demonstrated its role in the maintenance of the resulting diversity. We investigate the conditions that favour the origin and maintenance of alternative generalist and specialist resource-use phenotypes within the same population. Previously, we confirmed that competition for hosts among φ6 bacteriophage in a mixed novel (non-permissive) and ancestral (permissive) host microcosm triggered the evolution of a generalist phenotype capable of infecting both hosts. However, because the newly evolved generalists tended to competitively exclude the ancestral specialists, coexistence between the two phenotypes was rare. Here, we show that reducing the relative abundance of the novel host slowed the increase in frequency of the generalist phenotype, allowing sufficient time for the specialist to further adapt to the ancestral host. This adaptation resulted in 'evolutionary rescue' of the specialists, preventing their competitive exclusion by the generalists. Thus, our results suggest that competition promotes both the origin and maintenance of biodiversity when it is strong enough to favour a novel resource-use phenotype, but weak enough to allow adaptation of both the novel and ancestral phenotypes to their respective niches.

Coronavirus Host Range Expansion and Middle East Respiratory Syndrome Coronavirus Emergence: Biochemical Mechanisms and Evolutionary Perspectives.

Coronaviruses have frequently expanded their host range in recent history, with two events resulting in severe disease outbreaks in human populations. Severe acute respiratory syndrome coronavirus (SARS-CoV) emerged in 2003 in Southeast Asia and rapidly spread around the world before it was controlled by public health intervention strategies. The 2012 Middle East respiratory syndrome coronavirus (MERS-CoV) outbreak represents another prime example of virus emergence from a zoonotic reservoir. Here, we review the current knowledge of coronavirus cross-species transmission, with particular focus on MERS-CoV. MERS-CoV is still circulating in the human population, and the mechanisms governing its cross-species transmission have been only partially elucidated, highlighting a need for further investigation. We discuss biochemical determinants mediating MERS-CoV host cell permissivity, including virus spike interactions with the MERS-CoV cell surface receptor dipeptidyl peptidase 4 (DPP4), and evolutionary mechanisms that may facilitate host range expansion, including recombination, mutator alleles, and mutational robustness. Understanding these mechanisms can help us better recognize the threat of emergence for currently circulating zoonotic strains.

Dominance effects of deleterious and beneficial mutations in a single gene of the RNA virus ϕ6.

Most of our knowledge of dominance stems from studies of deleterious mutations. From these studies we know that most deleterious mutations are recessive, and that this recessivity arises from a hyperbolic relationship between protein function (i.e., protein concentration or activity) and fitness. Here we investigate whether this knowledge can be used to make predictions about the dominance of beneficial and deleterious mutations in a single gene. We employed a model system--the bacteriophage φ6--that allowed us to generate a collection of mutations in haploid conditions so that it was not biased toward either dominant beneficial or recessive deleterious mutations. Screening for the ability to infect a bacterial host that does not permit infection by the wildtype φ6, we generated a collection of mutations in P3, a gene involved in attachment to the host and in phage particle assembly. The resulting collection contained mutations with both deleterious and beneficial effects on fitness. The deleterious mutations in our collection had additive effects on fitness and the beneficial mutations were recessive. Neither of these observations were predicted from previous studies of dominance. This pattern is not consistent with the hyperbolic (diminishing returns) relationship between protein function and fitness that is characteristic of enzymatic genes, but could have resulted from a curve of increasing returns.

Comparison of SIV and HIV-1 genomic RNA structures reveals impact of sequence evolution on conserved and non-conserved structural motifs.

RNA secondary structure plays a central role in the replication and metabolism of all RNA viruses, including retroviruses like HIV-1. However, structures with known function represent only a fraction of the secondary structure reported for HIV-1(NL4-3). One tool to assess the importance of RNA structures is to examine their conservation over evolutionary time. To this end, we used SHAPE to model the secondary structure of a second primate lentiviral genome, SIVmac239, which shares only 50% sequence identity at the nucleotide level with HIV-1NL4-3. Only about half of the paired nucleotides are paired in both genomic RNAs and, across the genome, just 71 base pairs form with the same pairing partner in both genomes. On average the RNA secondary structure is thus evolving at a much faster rate than the sequence. Structure at the Gag-Pro-Pol frameshift site is maintained but in a significantly altered form, while the impact of selection for maintaining a protein binding interaction can be seen in the conservation of pairing partners in the small RRE stems where Rev binds. Structures that are conserved between SIVmac239 and HIV-1(NL4-3) also occur at the 5' polyadenylation sequence, in the plus strand primer sites, PPT and cPPT, and in the stem-loop structure that includes the first splice acceptor site. The two genomes are adenosine-rich and cytidine-poor. The structured regions are enriched in guanosines, while unpaired regions are enriched in adenosines, and functionaly important structures have stronger base pairing than nonconserved structures. We conclude that much of the secondary structure is the result of fortuitous pairing in a metastable state that reforms during sequence evolution. However, secondary structure elements with important function are stabilized by higher guanosine content that allows regions of structure to persist as sequence evolution proceeds, and, within the confines of selective pressure, allows structures to evolve.

Competition and the origins of novelty: experimental evolution of niche-width expansion in a virus.

Competition for resources has long been viewed as a key agent of divergent selection. Theory holds that populations facing severe intraspecific competition will tend to use a wider range of resources, possibly even using entirely novel resources that are less in demand. Yet, there have been few experimental tests of these ideas. Using the bacterial virus (bacteriophage) 6 as a model system, we examined whether competition for host resources promotes the evolution of novel resource use. In the laboratory, 6 exhibits a narrow host range but readily produces mutants capable of infecting novel bacterial hosts. Here, we show that when 6 populations were subjected to intense intraspecific competition for their standard laboratory host, they rapidly evolved new generalist morphs that infect novel hosts. Our results therefore suggest that competition for host resources may drive the evolution of host range expansion in viruses. More generally, our findings demonstrate that intraspecific resource competition can indeed promote the evolution of novel resource-use phenotypes.

Genetic architecture and the evolution of sex.

Theoretical investigations of the advantages of sex have tended to treat the genetic architecture of organisms as static and have not considered that genetic architecture might coevolve with reproductive mode. As a result, some potential advantages of sex may have been missed. Using a gene network model, we recently showed that recombination imposes selection for robustness to mutation and that negative epistasis can evolve as a by-product of this selection. These results motivated a detailed exploration of the mutational deterministic hypothesis, a hypothesis in which the advantage of sex depends critically on epistasis. We found that sexual populations do evolve higher mean fitness and lower genetic load than asexual populations at equilibrium, and, under moderate stabilizing selection and large population size, these equilibrium sexual populations resist invasion by asexuals. However, we found no evidence that these long- and short-term advantages to sex were explained by the negative epistasis that evolved in our experiments. The long-term advantage of sex was that sexual populations evolved a lower deleterious mutation rate, but this property was not sufficient to account for the ability of sexual populations to resist invasion by asexuals. The ability to resist asexual invasion was acquired simultaneously with an increase in recombinational robustness that minimized the cost of sex. These observations provide the first direct evidence that sexual reproduction does indeed select for conditions that favor its own maintenance. Furthermore, our results highlight the importance of considering a dynamic view of the genetic architecture to understand the evolution of sex and recombination.

Architecture and secondary structure of an entire HIV-1 RNA genome.

Single-stranded RNA viruses encompass broad classes of infectious agents and cause the common cold, cancer, AIDS and other serious health threats. Viral replication is regulated at many levels, including the use of conserved genomic RNA structures. Most potential regulatory elements in viral RNA genomes are uncharacterized. Here we report the structure of an entire HIV-1 genome at single nucleotide resolution using SHAPE, a high-throughput RNA analysis technology. The genome encodes protein structure at two levels. In addition to the correspondence between RNA and protein primary sequences, a correlation exists between high levels of RNA structure and sequences that encode inter-domain loops in HIV proteins. This correlation suggests that RNA structure modulates ribosome elongation to promote native protein folding. Some simple genome elements previously shown to be important, including the ribosomal gag-pol frameshift stem-loop, are components of larger RNA motifs. We also identify organizational principles for unstructured RNA regions, including splice site acceptors and hypervariable regions. These results emphasize that the HIV-1 genome and, potentially, many coding RNAs are punctuated by previously unrecognized regulatory motifs and that extensive RNA structure constitutes an important component of the genetic code.

Cross-sectional characterization of HIV-1 env compartmentalization in cerebrospinal fluid over the full disease course.

OBJECTIVES: To characterize HIV-1 env compartmentalization between cerebrospinal fluid (CSF) and peripheral blood plasma over all stages of the HIV-1 disease course, and to determine the relationship between the extent of CSF HIV-1 env compartmentalization and clinical neurologic disease status. DESIGN: Paired blood plasma and CSF specimens were collected from 66 HIV-infected patients cross-sectionally representing all major clinical stages relating to HIV-associated neurologic disease, including primary infection, asymptomatic chronic infection, chronic infection with minor global impairment, and immune deficiency with HIV-associated dementia. METHODS: Heteroduplex tracking assays and bulk sequence analysis targeting the V1/V2, C2-V3, and V4/V5 regions of env were performed to characterize the genetic makeup of complex HIV-1 populations in the cross-sectional blood plasma and CSF specimens. The levels of blood plasma/CSF env compartmentalization were quantified and compared across the different clinical stages of HIV-1 neurologic disease. RESULTS: Blood plasma/CSF env compartmentalization levels varied considerably by disease stage and were generally consistent across all three regions of env characterized. Little or no compartmentalization was observed in non-impaired individuals with primary HIV-1 infection. Compartmentalization levels were elevated in chronically infected patients, but were not significantly different between mildly impaired and non-impaired patients. Patients with HIV-associated dementia showed significantly greater blood plasma/CSF env compartmentalization relative to other groups. CONCLUSION: : Increased CSF compartmentalization of the HIV-1 env gene, which may reflect independent HIV-1 replication and evolution within the central nervous system, is specifically associated with HIV-associated dementia and not the less severe forms of HIV-1 neurologic disease.

Hotter is better and broader: thermal sensitivity of fitness in a population of bacteriophages.

Hotter is better is a hypothesis of thermal adaptation that posits that the rate-depressing effects of low temperature on biochemical reactions cannot be overcome by physiological plasticity or genetic adaptation. If so, then genotypes or populations adapted to warmer temperatures will have higher maximum growth rates than those adapted to low temperatures. Here we test hotter is better by measuring thermal reaction norms for intrinsic rate of population growth among an intraspecific collection of bacteriophages recently isolated from nature. Consistent with hotter is better, we find that phage genotypes with higher optimal temperatures have higher maximum growth rates. Unexpectedly, we also found that hotter is broader, meaning that the phages with the highest optimal temperatures also have the greatest temperature ranges. We found that the temperature sensitivity of fitness for phages is similar to that for insects.

Major coexisting human immunodeficiency virus type 1 env gene subpopulations in the peripheral blood are produced by cells with similar turnover rates and show little evidence of genetic compartmentalization.

A distinctive feature of chronic human immunodeficiency virus type 1 (HIV-1) infection is the presence of multiple coexisting genetic variants, or subpopulations, that comprise the HIV-1 population detected in the peripheral blood. Analysis of HIV-1 RNA decay dynamics during the initiation of highly active antiretroviral therapy (HAART) has been a valuable tool for modeling the life span of infected cells that produce the bulk HIV-1 population. However, different HIV-1 target cells may have different turnover rates, and it is not clear whether the bulk HIV-1 RNA decay rate actually represents a composite of the decay rates of viral subpopulations compartmentalized in different cellular subsets with different life spans. Using heteroduplex tracking assays targeting the highly variable V3 or V4-V5 regions of the HIV-1 env gene in eight subjects, we found that all detectable coexisting HIV-1 variants in the peripheral blood generally decayed at similar rates during the initiation of HAART, suggesting that all of the variants were produced by cells with similar life spans. Furthermore, single genome amplification and coreceptor phenotyping revealed that in two subjects coexisting HIV-1 variants with distinct CXCR4 or CCR5 coreceptor phenotypes decayed with similar rates. Also, in nine additional subjects, recombination and a lack of genetic compartmentalization between X4 and R5 variants were observed, suggesting an overlap in host cell range. Our results suggest that the HIV-1 env subpopulations detectable in the peripheral blood are produced by cells with similar life spans and are not genetically isolated within particular cell types.

Coinfection rates in Φ6 bacteriophage are enhanced by virus-induced changes in host cells.

Two or more viruses infecting the same host cell can interact in ways that profoundly affect disease dynamics and control, yet the factors determining coinfection rates are incompletely understood. Previous studies have focused on the mechanisms that viruses use to suppress coinfection, but recently the phenomenon of enhanced coinfection has also been documented. In the experiments described here, we explore the hypothesis that enhanced coinfection rates in the bacteriophage Φ6 are achieved by virus-induced upregulation of the Φ6 receptor, which is the bacterial pilus. First, we confirmed that coinfection enhancement in Φ6 is virus-mediated by showing that Φ6 attaches significantly faster to infected cells than to uninfected cells. Second, we explored the hypothesis that coinfection enhancement in Φ6 depends upon changes in the expression of an inducible receptor. Consistent with this hypothesis, the closely related phage, Φ12, that uses constitutively expressed lipopolysaccharide as its receptor, attaches to infected and uninfected cells at the same rate. Our results, along with the previous finding that coinfection in Φ6 is limited to two virions, suggest that viruses may closely regulate rates of coinfection through mechanisms for both coinfection enhancement and exclusion.

Beneficial fitness effects are not exponential for two viruses.

The distribution of fitness effects for beneficial mutations is of paramount importance in determining the outcome of adaptation. It is generally assumed that fitness effects of beneficial mutations follow an exponential distribution, for example, in theoretical treatments of quantitative genetics, clonal interference, experimental evolution, and the adaptation of DNA sequences. This assumption has been justified by the statistical theory of extreme values, because the fitnesses conferred by beneficial mutations should represent samples from the extreme right tail of the fitness distribution. Yet in extreme value theory, there are three different limiting forms for right tails of distributions, and the exponential describes only those of distributions in the Gumbel domain of attraction. Using beneficial mutations from two viruses, we show for the first time that the Gumbel domain can be rejected in favor of a distribution with a right-truncated tail, thus providing evidence for an upper bound on fitness effects. Our data also violate the common assumption that small-effect beneficial mutations greatly outnumber those of large effect, as they are consistent with a uniform distribution of beneficial effects.