SMARCAD1 Contributes to the Regulation of Naive Pluripotency by Interacting with Histone Citrullination.

Histone citrullination regulates diverse cellular processes. Here, we report that SMARCAD1 preferentially associates with H3 arginine 26 citrullination (H3R26Cit) peptides present on arrays composed of 384 histone peptides harboring distinct post-transcriptional modifications. Among ten histone modifications assayed by ChIP-seq, H3R26Cit exhibited the most extensive genomewide co-localization with SMARCAD1 binding. Increased Smarcad1 expression correlated with naive pluripotency in pre-implantation embryos. In the presence of LIF, Smarcad1 knockdown (KD) embryonic stem cells lost naive state phenotypes but remained pluripotent, as suggested by morphology, gene expression, histone modifications, alkaline phosphatase activity, energy metabolism, embryoid bodies, teratoma, and chimeras. The majority of H3R26Cit ChIP-seq peaks occupied by SMARCAD1 were associated with increased levels of H3K9me3 in Smarcad1 KD cells. Inhibition of H3Cit induced H3K9me3 at the overlapping regions of H3R26Cit peaks and SMARCAD1 peaks. These data suggest a model in which SMARCAD1 regulates naive pluripotency by interacting with H3R26Cit and suppressing heterochromatin formation.

Spatiotemporal clustering of the epigenome reveals rules of dynamic gene regulation.

Spatial organization of different epigenomic marks was used to infer functions of the epigenome. It remains unclear what can be learned from the temporal changes of the epigenome. Here, we developed a probabilistic model to cluster genomic sequences based on the similarity of temporal changes of multiple epigenomic marks during a cellular differentiation process. We differentiated mouse embryonic stem (ES) cells into mesendoderm cells. At three time points during this differentiation process, we used high-throughput sequencing to measure seven histone modifications and variants--H3K4me1/2/3, H3K27ac, H3K27me3, H3K36me3, and H2A.Z; two DNA modifications--5-mC and 5-hmC; and transcribed mRNAs and noncoding RNAs (ncRNAs). Genomic sequences were clustered based on the spatiotemporal epigenomic information. These clusters not only clearly distinguished gene bodies, promoters, and enhancers, but also were predictive of bidirectional promoters, miRNA promoters, and piRNAs. This suggests specific epigenomic patterns exist on piRNA genes much earlier than germ cell development. Temporal changes of H3K4me2, unmethylated CpG, and H2A.Z were predictive of 5-hmC changes, suggesting unmethylated CpG and H3K4me2 as potential upstream signals guiding TETs to specific sequences. Several rules on combinatorial epigenomic changes and their effects on mRNA expression and ncRNA expression were derived, including a simple rule governing the relationship between 5-hmC and gene expression levels. A Sox17 enhancer containing a FOXA2 binding site and a Foxa2 enhancer containing a SOX17 binding site were identified, suggesting a positive feedback loop between the two mesendoderm transcription factors. These data illustrate the power of using epigenome dynamics to investigate regulatory functions.