Functional divergences of coding genes can be caused by divergences in their coding sequences and expression. However, whether and how expression divergences and coding sequence divergences coevolve is not clear. Gene expression divergences in differentiated cells and tissues recapitulate developmental models within a species, while gene expression divergences between analogous cells and tissues resemble traditional phylogenies in different species, suggesting that gene expression divergences are molecular traits that can be used for evolutionary studies. Using transcriptomes and evolutionary proxies to study gene expression divergences among differentiated cells and tissues in Arabidopsis, expression divergences of coding genes are shown to be strongly anti-correlated with phylostrata (gene ages), indicators of selective constraint Ka/Ks (nonsynonymous replacement rate/synonymous substitution rate) and indicators of positive selection (frequency of loci with Ka/Ks > 1), but only weakly or not correlated with indicators of neutral selection (Ks). Our results thus suggest that expression divergences largely coevolve with coding sequence divergences, suggesting that expression divergences of coding genes are selectively fixed by natural selection but not neutral selection, which provides a molecular framework for trait diversification, functional adaptation and speciation. Our findings therefore support that positive selection rather than negative or neutral selection is a major driver for the origin and evolution of Arabidopsis genes, supporting the Darwinian theory at molecular levels.