Comparative studies of gene regulation suggest an important role for natural selection in shaping gene expression patterns within and between species. Most of these studies, however, estimated gene expression levels using microarray probes designed to hybridize to only a small proportion of each gene. Here, we used recently developed RNA sequencing protocols, which sidestep this limitation, to assess intra-and interspecies variation in gene regulatory processes in considerably more detail than was previously possible. Specifically, we used RNA-seq to study transcript levels in humans, chimpanzees, and rhesus macaques, using liver RNA samples from three males and three females from each species. Our approach allowed us to identify a large number of genes whose expression levels likely evolve under natural selection in primates. These include a subset of genes with conserved sexually dimorphic expression patterns across the three species, which we found to be enriched for genes involved in lipid metabolism. Our data also suggest that while alternative splicing is tightly regulated within and between species, sex-specific and lineage-specific changes in the expression of different splice forms are also frequent. Intriguingly, among genes in which a change in exon usage occurred exclusively in the human lineage, we found an enrichment of genes involved in anatomical structure and morphogenesis, raising the possibility that differences in the regulation of alternative splicing have been an important force in human evolution.[Supplemental material is available online at http://www.genome.org. The RNA-seq data have been submitted to the NCBI Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/) under series accession no. GSE17274.]Changes in gene regulation are thought to play an important role in adaptive evolution and speciation (Britten and Davidson 1971;King and Wilson 1975;Jin et al. 2001;Carroll 2003Carroll , 2008Abzhanov et al. 2004;Iftikhar et al. 2004;Shapiro et al. 2004;Taron et al. 2004;Wray 2007). In support of this notion, comparative genome-wide studies of gene regulation within and between populations and species have revealed evidence consistent with the action of both stabilizing as well as directional selection on gene expression levels (Oleksiak et al. 2002;Lemos et al. 2005;Rifkin et al. 2005;Gilad et al. 2006;Whitehead and Crawford 2006). Most of these studies, however, focused on estimates of overall gene expression levels, probably because prior to the development of next-generation sequencing, it was very challenging to characterize expression level variation of individual exons on a genome-wide scale.Indeed, previous studies of alternative splicing patterns in mammalian species focused on relatively small numbers of exons and genes. For example, Su et al. (2008) studied variation in exon usage and alternative splicing in liver samples from a number of mouse strains from both sexes, by using a custom microarray designed to probe the expression levels of 25,760 exons and exonexon junctions f...