Distinct binding pattern of EZH2 and JARID2 on RNAs and DNAs in hepatocellular carcinoma development

Hepatocellular carcinoma (HCC) is one of the most malignant cancers worldwide, with high mortality. However, the molecular regulatory mechanisms of liver cancer, especially transcriptional and post-transcriptional mechanisms, should be further studied. Here we used chromatin and cross-linking immunoprecipitation with high throughput sequencing methods (ChIP-seq and CLIP-seq) to capture the global binding profiles on RNAs and DNAs of Enhancer of zeste homolog 2 (EZH2) and its partner Jumonji And AT-Rich Interaction Domain Containing 2 (JARID2) in liver carcinoma cell lines (HepG2) and normal liver cell line (THLE-2), respectively. We also integrated HCC transcriptome data from the TCGA to analyze the expression pattern of bound genes. We found that EZH2 and JARID2 both showed distinct binding profiles between HepG2 and THLE-2 cells. By binding to the primary RNAs, bound transcripts of EZH2 and JARID2 in HepG2 showed significantly increased transcriptional levels in HCC patients. By performing gene set enrichment analysis (GSEA), the bound transcripts were also highly related to HCC development. We also found EZH2 and JARID2 could specifically bind to several long noncoding RNAs (lncRNAs), including H19. By exploring the DNA binding profile, we detected a dramatically repressed DNA binding ability of EZH2 in HepG2 cells. We also found that the EZH2-bound genes showed slightly increased transcriptional levels in HepG2 cells. Integrating analysis of the RNA and DNA binding profiles suggests EZH2 and JARID2 shift their binding ability from DNA to RNA in HepG2 cells to promote cancer development in HCC. Our study provided a comprehensive and distinct binding profile on RNAs and DNAs of EZH2 and JARID2 in liver cancer cell lines, suggesting their potential novel functional manners to promote HCC development.