Abstract
Although epidermal keratinocyte development and differentiation proceeds in similar fashion between humans and mice, evolutionary pressures have also wrought significant species-specific physiological differences. These differences between species could arise in part, by the rewiring of regulatory network due to changes in the global targets of lineage-specific transcriptional master regulators such as p63. Here we have performed a systematic and comparative analysis of the p63 target gene network within the integrated framework of the transcriptomic and epigenomic landscape of mouse and human keratinocytes. We determined that there exists a core set of ∼1600 genomic regions distributed among enhancers and super-enhancers, which are conserved and occupied by p63 in keratinocytes from both species. Notably, these DNA segments are typified by consensus p63 binding motifs under purifying selection and are associated with genes involved in key keratinocyte and skin-centric biological processes. However, the majority of the p63-bound mouse target regions consist of either murine-specific DNA elements that are not alignable to the human genome or exhibit no p63 binding in the orthologous syntenic regions, typifying an occupancy lost subset. Our results suggest that these evolutionarily divergent regions have undergone significant turnover of p63 binding sites and are associated with an underlying inactive and inaccessible chromatin state, indicative of their selective functional activity in the transcriptional regulatory network in mouse but not human. Furthermore, we demonstrate that this selective targeting of genes by p63 correlates with subtle, but measurable transcriptional differences in mouse and human keratinocytes that converges on major metabolic processes, which often exhibit species-specific trends. Collectively our study offers possible molecular explanation for the observable phenotypic differences between the mouse and human skin and broadly informs on the prevailing principles that govern the tug-of-war between evolutionary forces of rigidity and plasticity over transcriptional regulatory programs.