Directional selection at gene expression level contributes to the speciation of Asian rice cultivars

Differences in expression levels play important roles in phenotypic variation across species, especially those closely related species with limited genomic differences. Therefore, studying gene evolution at expression level is important for illustrating phenotypic differentiation between species, such as the two Asian rice cultivars, Oryza sativa L. ssp. indica and Oryza sativa L. ssp. japonica. In this study, we evaluated the gene expression variation at inter-subspecies and intra-subspecies level using transcriptome data from seedlings of three indica and japonica rice and defined four groups of genes under different natural selections. We found a substantial of genes (about 79%) that are under stabilizing selection at the expression level in both subspecies, while about 16% of genes are under directional selection. Genes under directional selection have higher expression level and lower expression variation than those under stabilizing selection, which suggest a potential explanation to subspecies adaptation to different environments and interspecific phenotypic differences. Subsequent functional enrichment analysis of genes under directional selection shows that indica rice have experienced the adaptation to environmental stresses, and also show differences in biosynthesis and metabolism pathways. Our study provide an avenue of investigating indica-japonica differentiation through gene expression variation, which may guide to rice breeding and yield improvement.