Protein-protein interactions may impose constraints on both structural and regulatory evolution. Here we show that protein-protein interactions are negatively associated with evolutionary variation in gene expression. Moreover, interacting proteins have similar levels of variation in expression, and their expression levels are positively correlated across strains. Our results suggest that interacting proteins undergo similar evolutionary dynamics, and that their expression levels are evolutionarily coupled. These patterns hold for organisms as diverse as budding yeast and fruit flies.Recent studies on variation in gene expression in natural populations indicate that a number of factors influence regulatory evolution, including membership in a particular functional class or biological process 1 and the pattern of sex-biased expression 2,3 . Protein-protein interactions may also be relevant to regulatory evolution if they mean that the interacting partners of a given protein impose stricter stoichiometric requirements. For proteins that participate in complexes, such requirements may result in increased selection pressure to maintain evolutionarily stable expression levels and might be manifested in a negative association between evolutionary variation in gene expression and the number of protein interactions.We analyzed variation in expression of 5,978 yeast (Saccharomyces cerevisiae) genes among four natural isolates grown under identical laboratory conditions 4 , of 4,439 fly (Drosophila melanogaster) genes in adult males from eight wild-type strains of geographically diverse origin 2 and of a similar number of genes assayed in both sexes of two species of Drosophila 3 . We matched these data to protein-protein interactions that could be assigned with high confidence (confidence score 4 0.55): 8,728 in yeast 5 and 3,964 in flies 6 . Choosing interactions with more stringent selection criteria produced similar results (Supplementary Note online).The breadth of variation of expression for a given gene in both flies and yeast (gene expression polymorphism) was negatively correlated with the number of protein-protein interactions found for the product of that gene (yeast: r ¼ -0.11, P o 0.0001, n ¼ 3,536; flies: r ¼ -0.15, P ¼ 0.022, n ¼ 225; Fig. 1). Gene expression divergence between two Drosophila species was also negatively associated with the number of protein-protein interactions (males: r ¼ -0.18, P ¼ 0.012, n ¼ 205; females: r ¼ -0.19, P ¼ 0.007, n ¼ 205). Furthermore, the connectivity of a protein was an excellent predictor of the mean variation in expression in classes of genes with a similar number of interactions ( Supplementary Fig. 1 online). Taken together, these results suggest that protein interactions might be a constraint that reduces evolutionary variation in gene expression both within and between species.Among yeast proteins, the rate of accumulation of replacement substitutions decreases as the number of protein-protein interactions increases 7,8 . The importance of this conclusion has been d...