It is generally assumed that stabilizing selection promoting a phenotypic optimum acts to shape variation in quantitative traits across individuals and species. Although gene expression represents an intensively studied molecular phenotype, the extent to which stabilizing selection limits divergence in gene expression remains contentious. In this study, we present a theoretical framework for the study of stabilizing and directional selection using data from between-species divergence of continuous traits. This framework, based upon Brownian motion, is analytically tractable and can be used in maximum-likelihood or Bayesian parameter estimation. We apply this model to gene-expression levels in 7 species of Drosophila, and find that gene-expression divergence is substantially curtailed by stabilizing selection. However, we estimate the selective effect, s, of gene-expression change to be very small, approximately equal to Ns for a change of one standard deviation, where N is the effective population size. These findings highlight the power of natural selection to shape phenotype, even when the fitness effects of mutations are in the nearly neutral range.evolution Í nearly neutral Í Ornstein-Uhlenbeck Í phenotypic optima A bundant evidence indicates that natural selection is remarkably powerful in shaping nucleotide sequences (1, 2). Many tests of natural selection rely on a comparison between nonsynonymous sites, in which mutations affect protein sequence, and synonymous sites, in which mutations do not. Synonymous sites serve as a proxy for neutral sites, enabling the effects of selection to be distinguished from background mutational and demographic patterns. Although changes in gene expression are hypothesized to play a major role in adaptation (3, 4), changes in expression cannot be so easily partitioned into neutral and selected categories. Thus, methods derived to analyze selection in coding sequences cannot be readily applied to gene-expression data. In part because of this ambiguity, general forces acting on gene-expression divergence have remained unclear. At this point, there exists considerable debate over the relative importance of selection and random drift in shaping gene-expression levels (5-8).The benefits of optimal gene regulation seem in many ways obvious. In the simple case of metabolic enzymes, underexpression may slow metabolic flux, while over-expression may expose the cell to additional toxic misfolded proteins (9). At the morphological level, gene regulation can be tightly coupled to phenotype (10, 11). Genetic mutations whose effects cascade into morphological differences are expected to have especially large fitness impacts, and as such will be heavily influenced by natural selection. A straightforward example of selection on gene-expression level can be seen in ribosomal proteins, which contrary to the neutral prediction are found to be highly expressed across a variety of organisms (12).In this article, we present a model of gene-expression divergence that explicitly distinguishes betw...