We have used comparative genomics to characterize the evolutionary behavior of predicted splicing regulatory motifs. Using base substitution rates in intronic regions as a calibrator for neutral change, we found a strong avoidance of synonymous substitutions that disrupt predicted exonic splicing enhancers or create predicted exonic splicing silencers. These results attest to the functionality of the hexameric motif set used and suggest that they are subject to purifying selection. We also found that synonymous substitutions in constitutive exons tend to create exonic splicing enhancers and to disrupt exonic splicing silencers, implying positive selection for these splicing promoting events. We present evidence that this positive selection is the result of splicing-positive events compensating for splicing-negative events as well as for mutations that weaken splice-site sequences. Such compensatory events include nonsynonymous mutations, synonymous mutations, and mutations at splice sites. Compensation was also seen from the fact that orthologous exons tend to maintain the same number of predicted splicing motifs. Our data fit a splicing compensation model of exon evolution, in which selection for splicing-positive mutations takes place to counter the effect of an ongoing splicing-negative mutational process, with the exon as a whole being conserved as a unit of splicing. In the course of this analysis, we observed that synonymous positions in general are conserved relative to intronic sequences, suggesting that messenger RNA molecules are rich in sequence information for functions beyond protein coding and splicing.[Supplemental material is available online at www.genome.org.]In addition to sequences that specify their cognate polypeptides, the open reading frames of eukaryotic messenger RNAs (mRNAs) are likely to contain additional information governing the functioning of these molecules; these functions include exon splicing and mRNA transport, localization, stability, and translation. If so, the evolution of mRNA sequence would be constrained by purifying selection against the loss of this information. This selection should be evident by examining the rate and quality of base substitutions that do not alter protein coding (synonymous substitutions) but may alter one of the abovementioned functions.Mutation rates in open reading frames have been extensively characterized by two parameters: K a , the rate of mutations that result in amino acid substitutions; and K s , the rate at synonymous sites. The latter has most often been assumed to reflect neutral change and has been used to normalize the mutation rates in a given gene, with the K a /K s ratio serving as an inverse measure of protein sequence conservation (Hurst 2002). More recently, the neutrality of K s has been challenged, with accumulating evidence for selection acting at synonymous sites (for reviews, see Chamary et al. 2006;Xing and Lee 2006). Some of this evidence has been based on comparisons between mutation rates in predicted exonic splicing e...