Abstract
Long noncoding RNAs (lncRNAs) have emerged as vital regulators of gene expression during embryonic stem cell (ESC) self-renewal and differentiation. Here, we systemically analyzed the differentially regulated lncRNAs during ESC-derived cardiomyocyte (CM) differentiation. We established a perspicuous profile of lncRNA expression at four critical developmental stages and found that the differentially expressed lncRNAs were grouped into six distinct clusters. The cluster with specific expression in ESC enriches the largest number of lncRNAs. Investigation of lncRNA-protein interaction network revealed that they are not only controlled by classic key transcription factors, but also modulated by epigenetic and epitranscriptomic factors including N(6)-methyladenosine (m(6)A) effector machineries. A detailed inspection revealed that 28 out of 385 lncRNAs were modified by methylation as well as directly recruited by the nuclear m(6)A reader protein Ythdc1. Unlike other 27 non-coding transcripts, the ESC-specific lncRNA Gm2379, located in both nucleus and cytoplasm, becomes dramatically upregulated in response to the depletion of m(6)A or Ythdc1. Consistent with the role of m(6)A in cell fate regulation, depletion of Gm2379 results in dysregulated expressions of pluripotent genes and crucial genes required for the formation of three germ layers. Collectively, our study provides a foundation for understanding the dynamic regulation of lncRNA transcriptomes during ESC differentiation and identifies the interplay between epitranscriptomic modification and key lncRNAs in the regulation of cell fate decision.